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ABSTRACT 


The main purpose of this thesis was to obtain the symbolic solution of a 
multilayer (four fluid media) ocean waveguide problem. The waveguide was 
assumed to have depth-dependent ambient density and sound-speed 
profiles in all fluid media, and arbitrarily shaped boundaries between all 
fluid media. A system of 28 equations in 17 unknowns was generated by 
satisfying all of the boundary conditions (including the boundary condition 
at the source) in cylindrical coordinates. The problem was set up as a 
weighted least squares problem for symbolic solution by the computer 
program Afathematica. Due to software and hardware constraints, a symbolic 
solution for the most general case was not obtained. However, by making all 
of the boundaries plane, parallel! boundaries, two cases were successfully 
program med, yielding symbolic solutions which were verified by comparison 


to previously known results. 
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I. A MORE GENERAL SOLUTION TO THE LINEAR, THREE- 
DIMENSIONAL, LOSSLESS, HOMOGENEOUS WAVE EQUATION 


The primary objective of this section is to derive the solution of the 
linear, three-dimensional, lossless, homogeneous wave equation (1) in the 
cylindrical coordinate system defined in Figure | for an arbitrary sound 


speed profile (a function of the depth coordinate, y, only). 


1 a@gltr) 
Vegttit) = 0 i, 
ote) Cy) eat ( 
The Laplacian expressed in the defined cylindrical coordinate system is 
given by (2). 


+ + + 
ar? tar rtate ay? 





,@ iat @ 


The remaining quantities of (1) are defined as follows: 
* g(t.r) is the velocity potential at time t at a position r = (roy) expressed 
in units of square meters per second, and 
¢c(y) is the arbitrary, depth-dependent speed of sound expressed in units 


of meters per second. 


The first assumption we will make in this derivation is that the source 
of acoustic energy in the waveguide has a time-harmonic dependence. This 
assumption may be justified by the facts that any arbitrary time dependence 


can be expressed as a summation of time-harmonic terms using Fourier 


analysis. In addition, one of the basic tenets of linear acoustics allows the 
velocity potential to be expressed as a linear combination (using the 


superposition principle) of such terms. 





Figure 1. The Cylindrical Coordinate System 


Therefore, the resulting velocity potential will also have a time- 


harmonic dependence given by 


(tr) = or) ebaatt (3) 
where, f is the source frequency in Hertz. 


Substituting (3) into (1) gives 


2 jaxft 
Vlglr) el2att] - _1 alge) el?) 0 


c*(y) at? @) 


Taking the time derivatives first yields 


oe San ae a » j2nf g(r) el2all 


alg(r) elt] a (alg) sal) alj2xf ir) ei!) 
at> st at : at 


= (j2nf)* g(r) el2att 


a*[o(r) el2aft] \ 


rv - wo? gr) el@att | (5) 


where w= 2nf is the angular source frequency in radians per second. 


Substituting (5) into (4) and observing that the Laplacian does not 
operate on the time-harmonic term e/**!t yields 
me 


c*(y) 





el2ait V-qir) ay gir) el2alt =O. 


Dividing out the complex exponential terms reveals the time- 


independent lossless Helmholtz equation 


Voir) + ky) gir) = 0, (6) 


where, 


k(y)- a SE (7) 
¥™ ely) 7 ely)” Aly) 


is the wave number expressed in units of radians per meter and A is the 


wavelength expressed in meters. 


The next step is to find the solution to the lossless Helmholtz equation 
(6). This will be accomplished using the method of separation of variables. 


We will assume that the solution for g(r) = g(r.o,y) has the form 


g(r) = pr.gy) = R(r) (6) Y(y). (8) 


Substituting (8) into (6) yields 


V-[R(r) (6) Y(y)] + k2(y) [Rir) (6) Yly)] = 0. 


Performing the Laplacian operation reveals 


a[R(r) (6) Y(y)] ; 


alR(r) (6) Yl a-[Rir) (6) Y(y)] 
ar? 


i 1 
r ar r? ap? 


: a[R(r) @(6) Y(y)] 


ap + k*(y) [R(r) (6) Y(y)] - 0. (9) 


The partial derivatives of (9) may be replaced with total derivatives 
since the functions R, 6, and Y are each simply functions of the single 


variables r, $, and y respectively. Continuing, (9) may be written as 








d2R(r) 1dR(r) 1 d7o(4) 
#19) 10") SB + 4) YO) FE ROO) YO) 
2 
+ R(r) o(@) 3 ) , £2(y) {Ric} @(6) Y(y)] = 0. (10) 


Dividing (10) by the product [R(r) (4) Y(y)] yields 


1 d@Rir) 1 dR(r) 1 d7@($) 1 d?¥(y) 


Rir) dr2 “rR(r) dr “r2e(g) dg?" Y(y) dy? 











+k*(y) =0. 


Separating out the depth dependence reveals 

















1 ay) | ee 1 d’R(r) 1 dR(r) 1 76/9) 


¥(y) dy? Rr) dr? eR) dr ery) age EY 


Since the left-hand side of (11) is a function of a single variable (y) and 


the right-hand side is a function of two variables (r, $), the equality can only 


be true if each side is equal to the same constant, say k?. Thus, the right- 
hand side of (11) becomes 


i d’R(r) 1 Rr) 1 d*olo) 2 





R(t) dr? —orR(r) dr r2@(g) dg? ~k,. 











Multiplying through by r* and rearranging reveals 





~ &(o) do? 





cr? d*R(r) or dR(r) 2 > 1 d*(4) 
r 


ee r 


R(r) dr? Rir) dr 


Since the left-hand side (LHS) of (12) is a function of a single variable 
(r) and the right-hand side (RHS) is a function of a different single variable 
(>), the equality can hold only if each side is equal to the same constant, say 


n*. Thus, the right-hand side of (12) becomes 


1 d’o(o) 


(>) do? 





(13) 


do(9) 


Equation (14) is a second order ordinary differential equation with the 


following exact solution: 


(9) = A, cos(ng) + B, sin(n¢) , 


(Gs) 


where Ay and By are in general complex constants whose values are 


determined by satisfying boundary conditions. 


Now, the left-hand side of (12) must be evaluated. Since the left-hand 


side of (12) must also be equal to n’, we have 


fa d°R(r) me cles) ) 





Ric) ar Rip) ode kere n?, 


Multiplying (16) Oars Ets en rearranging reveals 


d*Rir) 
dr? 








dR(r) F [2 


ral 
r dr 


Let 
R(r) = g(k, r). 


Then 


dR(r) dig(k, c)]  dig(k, r)] d(k, r) 
dr 8 dre Midi da 


or, 


(17) 


(18) 


dR(r) — die(k, rl 
de dik tr) oy) 





Additionally, 


d?R(r) 4 [dR(e)] il dtelk, 
dre mor’ dp - dev * “dtk f) 





{g(k, £)} 
i 


hdl 


or, 
d?R(r) 2 dle(k, rl ee 
dr? dic 
Substituting (18) through (20) into (17) yields 
d*[g(k dig(k 
2 Migtky tk Ke dates ol . Le?- Bate) r) (21) 


od(k, 2 of d(k, 1) 


Dividing (21) by k* yields 


dlgtk, ri) 1 


d(k, 1)? aE ts s ote ls (kpr) = 0 (22) 


Equation (22) is known as Bessel's differential equation, which has the 
following exact solution: 

g(k.r) = A, J,(k,r) + B, N,(k,r), (23) 

where J, and N, are the Bessel functions of the first and second kind 


respectively, and A, and B, are in general complex constants whose values 


are determined by satisfying boundary conditions. 


Since we let R(r) - g(k, r), the exact solution for R(r) can be written as 


follows: 


R(r) = Ay Jn (Ky) + Be Na(k,P). (24) 


Now, the left-hand side of (11) must be set equal to ke, revealing 


| dey) ky) = kK? 


Multiplying (25) by Y(y) and rearranging yields 
2 
sak [k2¢y) - k2] ly) = 0. (26) 
Let 
kity) = ky) - Ki (27) 
Using (27) to rewrite (26) yields 
2 
OO) By) vey) <0. (28) 


dy* 


Since the coefficient Ky) in (28) is an arbitrary function of depth, y, an 


exact solution to this differential equation cannot be determined. 


In order to continue this derivation of a generalized solution to the 
wave equation, we will assume that the solution to (28) can be determined 


by other means, and is given simply by 
Y(y) = Y(y). (29) 
In certain cases, the sound speed profile, c(y), may be a function such 
that an exact solution to (28) may be found (say by consulting tables). In 
other specific cases, approximations can be made in order to put (28) in a 


form whereby a known solution may be found. However, in the most general 


case, (28) will have to be evaluated numerically. 


Recall that (8) specifies the velocity potential as 


pir.¢y) - R(r) O(o) Yy). (30) 


Substituting (15), (24), and (29) into (30) yields the following general 


solution for the velocity potential: 
g{r.oy) = [A, Jo(k. r) +B. Na(k, r)| lA, cos(n) + B, sin(ng)| ¥y). Gt) 


Finally, recall that (3) described the velocity potential as a function of 


time and position as 


ptr) = or) elett (32) 


Therefore, the complete general solution for the velocity potential is 


given by 


g(tir.oy) = [A, Ja(kp er) + B, Na(k, r)| la, cos(n¢) + By sin(ng)| Y(y) el2aft (33) 


II. DESCRIPTION OF THE GENERALIZED OCEAN WAVEGUIDE MODEL 


The next logical step in this analysis is to apply the solution to the 
linear, three-dimensional, lossless, homogeneous wave equation developed in 
the previous section to the generalized ocean waveguide model. Before 
continuing with the mathematical derivation, a short description of the 


generalized ocean waveguide model chosen for analysis will be presented. 


As shown in Figure 2, space has been separated into four distinct 
media. Medium I (which may represent the air) is completely characterized 
by its density (p,(y)) and sound speed profile (c,(y)), which may both be, in 
general, arbitrary functions of depth coordinate, y, only. Medium II (which 
may represent the ocean water) is completely characterized by its density 
(p2(y)) and sound speed profile (c.(y)), which may both be, in general, 
arbitrary functions of depth coordinate, y, only. Medium III (which may 
represent the upper layer (sediment) of the ocean bottom) is completely 
characterized by its density (p3(y)) and sound speed profile (c3(y)), which 
may both be, in general, arbitrary functions of depth coordinate, y, only. 
Medium IV (which may represent a second fluid layer of the ocean bottom) 
is completely characterized by its density (p4(y)) and sound speed profile 
(c4(y)), which may both be, in general, arbitrary functions of depth 


coordinate, y, only. 


Media I and II are separated by a boundary (which may represent the 


ocean surface) at a depth of ys(r.6) meters, where y, may in general be an 


arbitrary function of horizontal range, r, and azimuthal angle $. Media I] and 
III] are separated by a boundary (which may represent the ocean bottom) at 


a depth of ¥p, (ro) meters, where Vp, may in general be an arbitrary function 


of horizontal range, r, and azimuthal! angle 4. Media II] and IV are separated 


by a boundary (which may represent the interface between two different 


layers in the bottom) at a depth of Y(t.) meters, where Yp, may in general 
be an arbitrary function of horizontal range, r, and azimuthal angle 9. 


Medium I] requires some additional examination. First, medium I] 
contains the omnidirectional, time-harmonic point source located at a fixed 
depth of y, meters and at a fixed horizontal range of zero meters. Medium I] 
is also further separated into two distinct subregions, labelled medium Ila 
and medium IIb, by an artificial boundary located at the depth of the source. 
This additional boundary is required to satisfy the source boundary 
conditions. The density and sound speed profiles for these subregions are 
completely defined by the generic medium II descriptions previously 


discussed. 
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Figure 2. Generalized Ocean Waveguide Model 
(Based on Ziomek (1991, Figure 3.9-1)) 
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III. DEVELOPMENT OF THE GENERAL WAVEGUIDE MODEL SOLUTION 


The primary objective of this section is to derive general expressions 
for the velocity potentials in media Ila and I]b based on the wave equation 
solution previously derived and the appropriate general boundary 


conditions. 


The next assumption required in this derivation is that the source of 
acoustic energy in the waveguide is an omnidirectional point source which is 
surrounded by the medium it is to excite acoustically. This assumption may 
be justified by the fact that the velocity potential field generated by an 
arbitrary acoustic source array may be expressed as the summation of the 
velocity potential fields generated by the individual omnidirectional sources 
which make up the array. Assuming that the source is surrounded by the 
medium it is to drive acoustically is reasonable for many practical acoustic 


systems. 


Recall that the velocity potential is given by 


troy) = [A, Gen) BaNe(ke r)| la, cos (n¢) + B, sin (ng) Y(y) ef@!" (3.4) 


The omnidirectional point source assumption implies that the source 
radiates acoustic energy equally well in all directions. However, since we 
have assumed an arbitrarily varying surface and an arbitrarily varying 
bottom for this waveguide, this assumption goes not imply that the velocity 


potential field lacks an angular dependence (that is, in general, the velocity 


potential at a given range and depth will depend on the azimuthal angle 
because of the differences in the interactions of the sound rays with the 
arbitrary surface and bottom). Assuming that the source is surrounded by 
the medium (and since the source is located on the y-axis, this assumption 
implies that the y-axis is surrounded by the medium) implies that the 
arbitrary constant n must be an integer in order to ensure that the velocity 
potential is single-valued for azimuthal angles $ in excess of 2x radians (ie., 


the velocity potential must be a periodic function with period 2x radians). 


Since the velocity potential must be evaluated at the source (ie, at r = 
0), the arbitrary constant B, must be set equal to zero in order to eliminate 
the Neumann function solution (since all of the Neumann functions tend to 
infinity as r tends to zero) (see Boas (1983, p. 513 and p. 525) for further 


proof). Thus, the general solution is reduced to 
gtrdy) - A, Jak, 6) la, cos (nd) + B, sin (ng)| Y(y) el2att (35) 
A. VELOCITY POTENTIAL IN MEDIUM I 
Based on (35), the velocity potential in medium | can be expressed as 


gp, (troy) = Ay Jalk r) lA, cos (nd) + B,, sin (ng) | Y,(y) el2att (36) 


where Y,(y) represents the solution to (28) for the traveling wave in the 


negative y direction. 


Carrying out the indicated multiplications, and defining the following 


new constants: 


reveals the general form for the velocity potential in medium | 


gy (troy) = [B, cos (no) + Ay sin (ng)| Jalk,, eh yiypelert, (37) 


where A, and B, are arbitrary complex constants whose values will be 


determined by satisfying boundary conditions. 


Additionally, 


2nf)? 
kily) = ke + ky (9) a 3 (38) 





where k,(y) is the wave number in medium I, Kr is the constant radial 
component of the propagation vector in medium I, ky ly) is the depth 


component of the propagation vector in medium I, and c,(y) is the depth- 


dependent sound speed in medium I. 


It is worth noting that there is not a solution of the form Y,(y) in (37). 
This is due to the fact that energy is considered to propagate in the negative 


y direction out to negative infinity without reflection. 
B. VELOCITY POTENTIAL IN MEDIUM II 


It is clear from the configuration that the velocity potentials in media 
Ila and IIb will be combinations of “incident” waves traveling toward the 
respective boundaries and "reflected" waves traveling away from these 
boundaries. With this in mind, the velocity potential in medium I]a may be 


written as (see equation (35)) 
Mraltr.oy) = a Jatk,., r) (a, cos nd + Bios sin al 


* Ayo, Yoaly) “ Baa Y;34(y)] alent : (39) 


where Y.,(y) represents the solution to (28) for the traveling wave in the 


Negative y direction (incident on the boundary at y = ys) and Y,,(y) 


represents the solution to (28) for the traveling wave in the positive y 


direction (reflected from the boundary at y = ys). 


Carrying out the indicated multiplications and defining the following 


new constants: 


A A 


ron Mba, vag Boa 
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Tos Sp Bey 


A C 


roa Pb2a Aya ™ S20 


Aco, Boe By oa i Dog 


reveals the general form for the velocity potential in medium Ila 
vraltiry) - (A, 005 n6 Yq(y) + B,, cos nd Yo,(y) 


+C,, sin noY>,(y) + D,, sin ng Yoav) Jat Kava) clea (40) 


where A, B,. C. and D.. are arbitrary complex constants whose values 
will be determined by satisfying boundary conditions. 
Similarly, the velocity potential in medium IIb is given by 
Pop t.r.dy) = A, cos no Y5,(y) + B.,, COs nb Y>,(y) 


+C,, sin no Yop(y) + D,, sin ng Y>,(y)) In{k,, #) eal (41) 


where Y,,(y) represents the solution to (28) for the traveling wave in the 


positive y direction (incident on the boundary at y - Yp, ) Yop(y) represents 


the solution to (28) for the traveling wave in the negative y direction 


(reflected from the boundary at y - Yp,) and Any: B C,. and Do. are 


2b’ 


arbitrary complex constants whose values will be determined by satisfying 
boundary conditions. 


Additionally, 


(2nf)? 


cy) 





ky) = key + kyg(¥) =e (42) 


where ko(y) is the wave number in medium I], ky, is the constant radial 
component of the propagation vector in medium II, ky,(y) is the depth 
component of the propagation vector in medium IJ, and co{y) is the depth- 
dependent sound speed in medium II. 


C. VELOCITY POTENTIAL IN MEDIUM III 


It is again clear from the configuration of our waveguide that the 
velocity potential in medium III will be a combination of a traveling wave in 
the positive y direction and a traveling wave in the negative y direction. 
Thus, the velocity potential in medium III can be derived using the 
techniques presented previously (for medium II). Performing this analysis 


yields 
p3(t0,>y) = {A, cos nd Y,(y) +B, cos no Y3(y) 


+C, sin ng Y,(y) + D, sin ng Yx(y)} In(k,, r) eb2att (43) 
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where Y3(y) represents the solution to (28) for the traveling wave in the 
positive y direction (incident on the boundary at y = Ypo)s Y,(y) represents 


the solution to (28) for the traveling wave in the negative y direction 


(reflected from the boundary at y = Ypo and A,, B,, C,, and D, are 


arbitrary complex constants whose values will be determined by satisfying 
boundary conditions. 
Additionally, 


2 2 2 (2xf)? 
k,(y) = kr, 2 y,(¥) oe 





c3(y) | 


where k3(y) is the wave number in medium II], kr, is the constant radial 
component of the propagation vector in medium ITI, ky,(y) is the depth 


component of the propagation vector in medium III, and c;(y) is the depth- 


dependent sound speed in medium III. 
D. VELOCITY POTENTIAL IN MEDIUM IV 
Finally, in medium IV, the velocity potential is given by (see (35)) 


g4(t.r.¢y) = A, Jafk,, r) A,, COS Nd + B,, sin no} Viv) cate (45) 


where Y,(y) represents the solution to (28) for the traveling wave in the 


positive y direction. 


Carrying out the indicated multiplications and defining the following 


new constants: 


A. Ay.® we As 


A, By, = Ba. 


reveals the general form for the velocity potential in medium IV 


pg(tr.oy) = (A, cos nd + B, sin no} Jatk,, AGA (46) 


where A, and Bs, are arbitrary complex constants whose values will be 


determined by satisfying boundary conditions. 


Additionally, 


2 ea 
ki (y) = K ae (y) = , (47) 





where k.(y) is the wave number in medium IV, k, ; is the constant radial 


component of the propagation vector in medium IV, ky Aly) is the depth 


20 


component of the propagation vector in medium IV, and c,(y) is the depth- 


dependent sound speed in medium IV. 


It is worth noting that there is not a solution of the form Y,(y) in (46). 


This is due to the fact that energy is considered to propagate in the positive 


y direction out to positive infinity without reflection. 
E. BOUNDARY CONDITIONS 


There are three different types of boundary conditions that must be 
applied to the solution of our problem. The first is the condition of 
“continuity of acoustic pressure’ across a boundary. This condition requires 
that the acoustic pressure evaluated at a particular spatial location and time 
on one side of the boundary be identically equal to the pressure evaluated at 
the same spatial location and time on the other side of the boundary. Kinsler 
(1982, p. 125) describes this condition as meaning that there is no net force 


acting on the massless boundary separating the two media. 


The second type of boundary condition is that of continuity of the 
normal component of the acoustic particle velocity across a boundary. This 
condition requires that the normal component of the acoustic particle 
velocity evaluated at a particular spatial location and time on one side of the 
boundary be identically equal to the normal component of the acoustic 
particle velocity evaluated at the same spatial location and time on the other 
side of the boundary. Kinsler (1982, p. 126) describes this condition as 


meaning that the media remain in contact with each other. 


Z| 


The final type of boundary condition is that of discontinuity of the 
normal component of the acoustic particle velocity across the boundary at 
the depth of the source. Officer (1958, p. 124 and following) and Ziomek 
(1991, discussion following equation (3.9-37)) describe this condition as 
being required to ensure that the solution to the wave equation reduces to 
that of an omnidirectional point source when the boundaries at y - ys, 


y Sa, and y = Yp, are removed. 


In developing the required boundary condition expressions, we will 
consider the various types of conditions in the order previously discussed. 
Thus, the first boundary condition to be applied to this problem is that of 


continuity of acoustic pressure across the boundary at y = ys. This implies: 


Py(t.r.b.ys) = Poalt.r.oys) . 
In general (from the fluid dynamics derivation), the acoustic pressure 


can be related to the velocity potential by the following: 


Se nc ee ae (48) 


where po(r.o,y), representing the ambient (equilibrium) density of the 


medium, is, in general, an arbitrary function of the spatial variables r, 4, y, 


but is not a function of time. 
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Since we have assumed that density is simply a function of depth 
coordinate y, and assumed a time-harmonic source, the acoustic pressure 


may be written as 


p(tr.dy) = -jwpoly) otroy). (49) 


Thus, the acoustic pressures in the various media may be expressed as 


p(troy) = -jwp;ly) o(troy), (50) 
Paaltir.gy) = -j opaly) gaaltr.oy) , (51) 
Pap(tr.oy) = - jo paly) gap(trey), (52) 
paltir.oy) = - j wp3ly) w(trgy), (53) 
and, 
p4(t.r.oy) = -jopgly) gglt.r.oy) , (54) 


where p,(y), poly), p3(y), and ply) are the ambient densities in media I, I], 


III, and IV, respectively. 


Returning to the specific boundary condition being examined, we set p, 
(given by (50)) equal to po, (given by (51)) at y = ys, and divide out the 


common factor of -j w revealing 
p1l¥s) gi{tr.b.Ys) = palys) oalt.r.g.ys) . (55) 


Substituting the derived velocity potentials for media I (37) and Ila 
(40) into (55) yields 
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p1(Ys) {B, cos n+ A, sin no} Ja(Kr, £) ¥y(ys) el2att 
= plYs) {A,, cos nd Y>a(¥s) + B,, cos no Yoal¥s) 
+C,, sin ng Y5_(¥s) + D,, sin ng Y>_(Vs)) Jake, ©) el2alt (56) 
The time dependence is eliminated by dividing (56) through by the 


complex exponential term el2alt, Carrying out the indicated multiplications 


and factoring yields the following: 
[oolys) Jakr, r) Y5,(Ys) A,, i £2(Ys) Jalkr, r) Ya(¥s) B,,| cos no 
+ [ostys) Inte, ©) ¥3q(¥s) C,, + pals) Ja(ke, ©) Yog(¥s) D,,] sin ng 
= pil¥s) Inlke, £) ¥yCys) By cos no + pylYs) In(kr, 1) ¥y(ys) Ay sin ng. (57) 
Setting the coefficients of cos nd on the left-hand side of (57) equal to 
the coefficients of cos nd on the right-hand side of (57), setting the 
coefficients of sin nd on the left-hand side of (57) equal to the coefficients of 
sin nd on the right-hand side of (57), and rearranging the equations 


generated by this process reveals the pair of equations representing the first 


boundary condition (BC #1) 
p2(Ys) Jn(kr, r) Yoal¥s) Xs 2 PolYs) Ja{kr, r) Yo4(¥s) Bs 


- p1(¥s) Jake, £) ¥i(ys) By = 0. (58) 
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and 
p2lys) Jntkr, ie), Y3,(ys) on : p2lys) Jatkr, r) Y,(¥s) ye 
zs p1(¥s) In{kr, r) Y (ys) Ay =Q. (59) 


It should be noted here that (58) and (59) are valid only if the 
associated trigonometric function is not identically zero for all values of 
azimuthal angle 9 (for instance, if n = 0, then sin nd is identically zero for all 


values of ¢ implying that (59) is no longer a valid boundary condition). 


The second boundary condition to be applied to this problem is that of 


continuity of acoustic pressure across the boundary at y = y,. This implies 


Paalt.r.b.Yo) = Paplt.r.d.yo) - 


Setting po, (given by (51)) equal to pop (given by (52)) at y - y, and 


dividing out the common terms reveals 


Pralt.l.o.Yo) = Papltr.gyYo) - (60) 


Substituting the derived velocity potentials for media Ila (40) and Ilb 
(41) into (60) yields 


Jatkr, r) ee cos np Y5,(Vo) + Bag cos nd Y>,(Yo) + Cog sin no Y>,yo) 
+ Dog Sin no Y>V¥o)| leat = Jalke, 1) [Asp cos n¢ VERE) 


+ Boy COS NO Yop (Vo) + Cop Sin no Yop(¥o) + Dap sin no Y34(¥o)! elf (61) 


2) 


Dividing out the common terms eliminates both the horizontal range 


and time dependences, revealing, after factoring, the following: 
Ly35(¥0) Ao, * YaalVo) B,,| COS Nd + Ly3,(v¢) C,. * Yoal¥o) D,,| sin nd 
= Ly5 (V0) A,, * Yonl¥o) Boy COS Nd + Ly3,(¥0) Con * YonlVo) Dap! sin np. (62) 


Setting the respective coefficients of cos nd and sin nd equal and 


rearranging yields the following pair of equations representing the second 


boundary condition (BC #2): 
Y2—(¥o) Age ‘ Yoal¥o) Bo, . Yop(¥o) A, 7 Yop(Vo) B., = UF (63) 
and 


You) Coe Yo alo) Dy oon pom onW aD aren (64) 


Again, (63) and (64) are valid only if the associated trigonometric 


function is not identically zero for all values of azimuthal angle 9. 


The third boundary condition to be applied to this problem is that of 


continuity of acoustic pressure across the boundary at y = YBr This implies 


Panltr.dyg,) - Palt.reyg)- 
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Setting po, (given by (52)) equal to pz (given by (53)) at y - Yay and 


dividing out the common factor of -j w reveals 


p2l¥p,) g20lte.d,Ye,) - pal¥p,) galLCoYp, } ; (65) 


Substituting the derived velocity potentials for media IIb (41) and III 
(43) into (65) yields 


p2l¥g,) Jn(ke, 1) [An cos nd You(¥p,) * Bop cos np Yo, (Yp,) 
+ Cop sin nd You(¥p,) + Dap sin nd You(¥p,)! el2xtt 
= p3lY¥p,) Jntkr, r) [A; cos nd Ya(¥p,) + Bz cos nd Y3(¥p,) 
+ C3 sin no Y3(¥p,) + D; sin no ¥3(¥5,)! glare (66) 
The time dependence is eliminated by dividing (66) through by the 


complex exponential term e!@*!" Carrying out the indicated multiplications 


and factoring reveals the following: 
lpalyg,) In{kr, £) Yon(¥p,) nat polyp.) Jake, 0) Yoo(¥p,) B,, | cos nd 


: lpolyp, ) Jake, r) Y2o(¥p,) C.. ; pal¥p,) Intkr, r) Yoo(Yp,) D,»! sin nd 
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~ lpsl¥p,) Jatkegt) Vglvp,) Ag * al¥p,) Jn(Ke, t) Ya(vp,) By | cos np 


+ loal¥p,) Jnlkr, r) Y3(¥p,) C, + pal¥p,) Jnlkr, r) Y3(¥p,) D, sin Nd. (67) 


Setting the respective coefficients of cos nd and sin nd equal and 
rearranging yields the following pair of equations representing the third 


boundary condition (BC #3): 
p2(¥g,) In(Kr, 0) Yop(¥p,) Any * O2l¥p,) In(Key 0) Yong.) Bop 
- £3(¥p,) In(ke, 1) Y3(¥p,) Ag - 03(¥p,) Jnlkr, 0) Ya(vp,)Bg=0, (68) 
and 
p2l¥p,) Ialkr, F) Yool¥p ) Coy * Pal¥p,) In(kry ) Yap(¥g,) Day 


- p3(¥g,) In(Kr, 0) Ya¥g.) Cz - o3(¥g,) In(ke,t) Yglyg.)D3-0. (69) 


Again, (68) and (69) are valid only if the associated trigonometric 


function is not identically zero for all values of azimuthal angle 9. 


The fourth boundary condition to be applied to this problem is that of 


continuity of acoustic pressure across the boundary at y = YB, This implies 


Pa(tr.oyp,) = Paltrdyp.) - 
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Setting pz (given by (53)) equal to p, (given by (54)) at y = Vp, and 


dividing out the common factor of -j w reveals 


pal py) aaltirbyp,) ~ Pal¥p,) altirbyg,)- (70) 


Substituting the derived velocity potentials for media II] (43) and IV 
(46) into (70) yields 


p3(¥p,) Jatkr, r) lA, Cos nd Y3(¥p,) + Bz cos nd Y3(¥p,) 
; + ; < Oy 

+ Cz sin nd Y3(¥p,) + Dz sin no ¥3(¥p,) el2xit 
- p4l¥p,) Ja(kr, 1) LA. cos no Y4(¥p,) + B, sin no Yi(¥5,) elésft = (74) 
The time dependence is eliminated by dividing (71) through by the 


complex exponential term e/2*!" Carrying out the indicated multiplications 


and factoring reveals the following: 
loalvp,) Jake, r) Y3(¥p,) A,+ p3l¥p,) Jntkr, r) Y3(yz,) B, | cos nd 
+ loalvp,) Jn(kr, t) Ya(yp,) C,+ p3l¥p,) Jnlkr, r)Ya(yp,) D, | sin nd 


- lpa(vp,) Intke, r) Y4(¥p,) A, | cos n¢ + lpalvs,) Intkrg r)Ya(¥p,) B, | sin nd. 


(72) 


29 


Setting the respective coefficients of cos n@ and sin nd equal and 
rearranging yields the following pair of equations representing the fourth 


boundary condition (BC #4): 
03(¥p,) In(kr, 1) Y3(vp,) Ag + 3(¥g,) In(kr, 1) Y3(¥p,) Bg 
- p4l¥g,) In(kryt) Yap.) Ay =O. (73) 
and 
p3(¥p,) Jar, F) Y3(¥p,) C, + 3(¥g,) Jake, £) Ya(yp,) D, 
- p4(¥g,) Jake, £) Y4l¥p,) B, = 0. (74) 
Again, (73) and (74) are valid only if the associated trigonometric 
function is not identically zero for all values of azimuthal angle 9. 


The fifth boundary condition is that of continuity of the normal 


component of acoustic particle velocity across the boundary at y = ys. This 


implies 
Uni(t.r.¢.ys) = Unoalt.r.d.ys) , 
where: 
Uni(trpy) = Uy(trgy) © ng(r.by), (75) 
Unaaltit.py) - Urltr.oy) * ag(roy), (76) . 
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A 
and n.(r,¢,y) represents the unit vector normal to the boundary at y = ys. 


In order to continue with the evaluation of this boundary condition. 
expressions for the velocities and the unit normal vector must be developed. 
Recalling that the acoustic particle velocity is simply the gradient of the 


velocity potential implies 


Ui(t.r.oy) = V gy (t.roy). (77) 


where g,(t.r.¢,y) is given by (37), and the gradient of g, expressed in 


cylindrical coordinates, is given by 


dg(t.r.4,y) 


ag(tr.g,y) a  agfltsr.d,y) a 
ar a 


V gltr.dy) = ret % re (78) 


Performing the indicated partial differentiations results in the 
following: 


Gtk t) 
ag (tr.4,y) ; = 
rail = Ge cos np + A, sin ng] erielse dr 


OJalke, r) d(Kky, r) 





- {B, cos no+ Ay sin no} ¥ (yet 


d(Ke r) ws 
di (car) 
Ege =k, {B, cos ng+ A, sin no} YGjea a (79) 
d(kr, r) 


oul 


ary - {- n B, sin ng+n A, cos ng) Jalke 0) ¥ j(y) e!*1* , geeealiaia) 


and 
dY,(y) 





agi (tt my) - {B, cos no + A, sin no} In(Kr r) elaxft 


ay (81) 


Substituting (79) through (81) into (77) (and using (78)) yields 


dJalke, r) , 
U,(t.r.¢y) = lke, {B, cos np + A, sin no} Y,(y) = ae 
d(kr, r) 
~ > LAY cos nd - B, sin no} Jalkr, r) Y,(y)¢ 
dY,(y) 
dy 





+ {B, cos no+ A, sin ng} Ja(kr, £) y\ elaatt (82) 


Conducting a similar analysis on the appropriate expressions reveals 


the following set of equations for the velocities in the remaining media: 
U2, (troy) = lk, { Avs cos nd Y,,(y) + Bog cos no Y>,(y) 


d]qkr, r) ° 
+ Co, sin nd Y>,(y) + Do, sin nd Y>a(v)d —— at 
ry r) 


n : ; 
or Ge cos nd Yo,(y) + Dz, cos nd Y,,(y) 
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- Ao sin nd ene. - Bo, sin nd Y3,(y)) In(kr, re 


as ) dY,,(y) 
+ Bo, cos nd 








+ (Aos cos Nd dy 


dY,,(y) dY,,{y) 
* Coq Sin 06 Ga * Dag sin 2 GS dy ——} }, In{Kr, 0) ¥ y Jett (83) 


Ua, (troy) - ike Lise COS Nd Y>(¥) + Bo, cos nd Y>,ly) 


WpqlKr, r) , 
+ Cop sin no Y>,(y) + Da, sin nd Y>,(y)J ————r 


d(kr, r) 


n ' 2 
an on cos No Yo,(y) + Dap cos no Y,(y) 
= Aap sin nd nn s Bop sin nb Y>,(y)3 Infkr, r) } 


wi dY>,(y) 
+ Bop Cos nd 








+ {Aop cos n¢ 


dY5,(y) dY>,(y) , 
“Co cinh = wee Dop sin nd ay ——} Jake, ¥ | eal (84) 


dy 


U3(tr.oy) = Ike, (A, cos nd Y,;(y) + Bs cos nd Y;(y) 


AylKr, c 
+ C3 sin no Y3(y) + D3 sin né Y(y)} 
d(kr, r) 


oe) 


+ = {Cy cos no Y;ly) + D3 cos ng Y3(y) 


- Az sin ng Y,(y) - Bz sin ng Y3(y)} Jake, r)¢ 


dY;(y) dY;(y) 
+ Bz cos nd dy 





+ | Ag cos no a 


dY;(y) dY,(y) a 
+ D3 sin nd re. Intkr, r)y el2att (85) 





+ Cz sin nd dy 


and 
| ; Walker, r) F 
U,(tr.oy) = lke, hy cos nd + By sin no} VG Soa oa 
d(ky ; r) 
+ = {B, cos nd - A, sin no} Y,,y) In{Kr, re 
dYz(y) 
dy 





+ ON cos no + B, sin no} Jn(Kr, 1) y\ el2att (86) 

Now, we must turn our attention to deriving an expression for the unit 
normal vector. A review of texts covering calculus and analytic geometry 
(for instance, Berkey (1988, p. 830), and Leithold (1972, p. 934)) remind us 
that if the surface can be expressed as a constant function of all three spatial 
variables (ie., f(r.¢.y) = a constant), then the normal vector to this surface at 
any point is simply the gradient of the function describing the surface 


evaluated at that point. Recalling that both the boundaries ys, YBy and YBo 


were defined to be arbitrary functions of the other spatial variables, let 
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y~-fs(ro), (87) at the ocean surface (i.e., y = ys), 


y- fp (r.@) , (88) at the ocean bottom (ie. y = Vp, and 


y= fp,{t.¢) , (89) at the bottom layer interface (ie., y = Yp,): 


Concentrating on the ocean surface for the moment, subtracting fs(r,9) 


from both sides of (87) yields 


y -f.(r.o) = 0. 


(90) 


Letting a new function, t.{r,¢,y), equal the left-hand side of (90), we 


have an equation of the form 


t(r.py) - y - fs(r.g) = 0. 


(91) 


Taking the gradient of t(r.o,y) using the cylindrical coordinate system 


gradient operator (78) yields 








atc(r.o.y) atc(r.o,V) dtc(r by) 
yadrey) GLEE ee Ren es 
ar rab oy 
afs(r.o) A ls af(r.o) A A 
Vv t(r.,y) 7 or cr 7 r ad % = y 
Therefore, the normal vector is given by 

f(r.o)a 1afslrgha a 
oe je ot ab oy 
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(92) 


(93) 


Since we need a unit normal vector, the magnitude of Nc (for which we 
use the symbol Nc) must be evaluated. Recalling that the magnitude of a 
vector is the square root of the sum of the squares of the individual 


components, then Nc can be written directly as 





(2uica)}2 ; 4 (Saar a (94) 


Therefore, the desired unit normal vector is given by 


A 1f oOfelr.o)a 1 afs(r.d) a 
Ds = ees oy (95) 


Similar analysis reveals the following expressions for the unit normal 


vectors at the boundaries y = YB, (n, > and y = YB, (ng, ): 

















Sees i A A 
eae fe Apert age ab 4 fh (96) 
By 
where 
af, (ro) af, (ro) 
1 2 1 1 7 
Ne, ar } + r2 { ab } + es (97) 
and 
’ 1{ ap (r4) 1a ) a 
oe ins | cal ad e+yJ. j 
Bo 
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where 






af, (0.9) af, (r.4) 


ap 





re 


=: Weenie (99) 


Returning to the evaluation of the fifth boundary condition, substituting 
(82) and (95) into (75), substituting (83) and (95) into (76), and performing 
the indicated dot products yields 


dqlkr r) 
Uns (troy) = i {B, cos nd+ A, sin no} Y,(y) ae 9 
dik, r) 


+ - {A, cos nd - B, sin no} In{kr r) Y,(y)¢ 


dY (y) , 
. {B, cos no + A, sin no) Intkr, #) ane Ye el2xtt 


Lf aflroha 1oafslrg)a a 
Ng oar ee ot ) i) 
Walker, r) 
afs(r, 
Uni (t.r.by) = [- ky {B, cos no+ A, sin nb! Y,(y) ants IN 
d(ky, r) : 
a afs(r.o) 


Nn ; 
2 (A, cos ng - B, sin ng) Jn{kr 6) ¥y(y) .:s 


ass ei2att 
dy ° Ng ’ 





+ {B, cos no+ A, sin no} In{Kr, r) (100) 


oF 


Unralt.e.o.y) = Ike, ene cos nd Y>,y) + Bag cos nd Y>,(y) 


dJalkr, r) ; 
+ Coq sin no Yo,(y) + Dg, sin no Y,,(y)J ————— r 
ry 0) 


Nn + a 
7 {Coa cos nd Yo,(y) + Do, cos nd Y,,(y) 


- Ao, sin no Y,,(y) - By, sin nd Y>,(v)4 Jn(Kr, r) ¢ 


dY,,(y) dY,,(y) 
+ Bog COS Nd dy 





+ (ae cos nd 


dY,,(y) dY,,(y) 
+ Dog Sin nd dy ——} }, Jnl kr, ¢) y ye 





+ Co, sin nd Ieee 





(r.o)a 1 afs(r.d) 
x ssa Lashes 49) 


Unaalt.r.oy) = [- kr, {Avs COs Nd Ya(y) + Bog COS Nb Y>,(Y) 


Gnlkr, r) afc(ro) 


+ Coq sin no Y>,(y) + Dag sin no Y>(¥)) 5 
r) 


n ' : . ' 
3 er cos np Yo,(y) + Dz, cos no Y>,(y) - Aza sin no Y2,(y) 


af(r.o) 


a Boa sin Nd Y>,(v)) In{Kr, r) ao 
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lie ) = dY,,(y) 








“LAB cos nd + Bog COS Nd + Co, Sin Nd dy 
dY.,(y) cia 
+ Dog sin nd ——} }, KO (101) 


Setting U,j(tr.oys) (100) equal to Unog(tr.oys) (101), and eliminating 
el2ait 
the common term Nc yields 





dake, r) 
z f(r, 
- kr {B, cos n>+ A; sin no} Y ys) —- ie 
f 
I 
s felr, 
r z {A, cos nd - B, sin no} In{kr r) Y,(ys) mate. = ) 
dY (Ys) 


+ {B, cos no + A, sin no} Jn(Kr, r) dy 


PK is cos no Y5,(Ys) + By cos nd Y>,(¥s) 
dnlKr, g 


+ Co, sin ng Y5,(ys) + Dog sin nd Y5,(ys) 
d{kr, 1) 


) ater.) 
or 


n . 9 
ee ice cos no Y»,(ys) + Dog cos nd Yo, (ys) 


afc(r.o) 


- Aza Sin ng Yoa(¥s) ~ Baa sin nb Yo4(¥s)} Jn(kr, 1) ; 


o2 


ae “s dY3,(Ys) 


























+ (Abe COs Nd + Bo, COS nd + Cy, sin nd dy 
a 
+ Dog sin nd = Jnlkr, ) . (102) 
Factoring (102) yields 
dake, 9) atte) afsir) , 
= Kr, me 1) = a Be Y4 (ys) rr - Ay 
(kr, 1) 
dY, (ys) 
d)alke, 1) afs(r.o) afs(r.o) 
+ [-ke, ASG) wrrmel ty: alls Yi) Ss, 
d(ke, r) : 
dY, (ys) 
+ Jalke, £) a A,| sin no 
dIn(ke, 1) Fear dIa(ke) ar rg) 
= 17 ke) ————~ Yaal¥s) ap 2a Keg 2al Ys) ar ee 
d(ky, rc rv) r) 
+e Jalkr, t) Yoalys) a Oa 2 citi, r) Y3(¥s) = oa Da 
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dY,(¥s) dY,,(y5) 


+ InlKr, Fr) dy Ada * Jq(Kr, 0) oy ie Boal cos nd 


dJalkr, r) 


+ 


2al¥s) a Coa kr, 2alYs) 5, 





_ E Kr, 
d{kr, r) 


) afc(r.g) 


afs(r.o) 
Boa 
ag : 


D a 
Aaa +79 Jalkrg£) Yoel) og 


Nn + 
si oA Jalkr, £) Yalys 
dY5,(¥s) dY2,(¥s) 


‘ Jalke, r) C5, # Intke, r) dy Daal SiN Nd. (103) 


dy 


Setting the respective coefficients of cos nd and sin nd equal and 


rearranging yields the following: 


dY {(ys) dJn(Kr. 1) 
. flr, 


ni : afs(r.o) 
= 7-2 Jalke, r) Y, (ys) Ay <a 


dY5,(¥s) dY,, (ys) 


es JnlKr, r) A2at In{kr, r) ay Boa 


dy 


4] 


da(Kr, r) dale, r) 


s fe ie Y>,(Ys) Ada ai kr, ay Y>,(¥s) Boa 
) d(kr, 1) 


| ice? 


- 7 [dale 1) r) Y3a(¥s) Co, + Ja(Kr, f r) Yaa Ys) ) Daal eae (104) 


and 


dY (ys) An(kr, 1) 


A a oe '¥s) Ay, 


afs(r.o) 


Nn s 
+ 0 Jar, £) Y(ys) By 8, 


dY,,(ys) dY,,(ys) 


= Ink. r) Co + Jntkr, r) Daal 


dy dy 


OJnlkr, r) Anlkr, r) 


- ks, ——— Ygqlrs) canes — eee oi) De 
) d Kr, r 


| af<(r.o) 


; : afs(r.o) 
: 2 [Ja(key#) Ypa(¥s) Aza * Ja(Ke, £) Yoq(¥s) Baal ras . (105) 


Again, (104) and (105) are valid only if the associated trigonometric 


function is not identically zero for all values of azimuthal angle $. It should 


be noted here that both (104) and (105) are of the form 
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(LHS term 1) + (LHS term 2) meg) (LHS term 3) a 
f(r, fo(r.) 
= (RHS term 1) + (RHS term 2) aise + (RHS term 3) ae 


Therefore, we may simplify (104) and (105) by setting LHS term 1 
equal to RHS term 1, LHS term 2 equal to RHS term 2, and LHS term 3 equal 
to RHS term 3. Performing this analysis and rearranging the resulting 
expressions yields the set of six equations representing the fifth boundary 


condition (BC #5) 


dY,,(y5) dY,,(ys) dY (ys) 
Jake, r) ian Ara + Jnlke, r) Ta - Jake, c) Loe: = 0, (106) 
dJn(ke, 1) dJn(ke, 1) 
Kr, Yoal¥s) Ara+ kr, Yoa(¥s) Boa 
d{ky, r) d(kr, r) 
dJq(Kr, r) 
- kp, ——— ¥ (ys) B, = 0, (107) 
dike 1) 


Jntkr, r) Y>,(¥s) Coase Jntkr, r) Y>,(¥s) Dog - Jntke r) Y ;(ys) A,;-90, (108) 


dY,,(¥s) dY,,(Ys) dY, (ys) 
Jalkr, r) dy Co, 3 Jntke, r) dy Dog a Jntke, r) ay Ay a 0 ' (109) 
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kr Vo .(y6) Coa + Kr, Yoal¥s) Dea 
d(k- r ky r 
2 2 
dJalke, r) 
= Ke enn Says ice Oe (110) 
d(k, r 


and 


Ja(ke, t) Yoal¥s) Aza * Ja(ke, t) Yoa(¥s) Boa Jake, 1) Yy(ys)By=0. (111) 


afc(r.) 
or 





Here, (107) and (110) are valid only if is not identically zero for 


all values of range, r, and azimuthal angle, $. Similarly, (108) and (111) are 
afc(r.o) 
valid only if ae is not identically zero for all values of range, r, and 


azimuthal angle, 9. 


The sixth boundary condition is that of continuity of the normal 


component of the acoustic particle velocity at the boundary y = Yer This 


implies 


Unaplt.r.oyg) = Unalt.r.oyp ) , 


where 


Unaltroy) = Usltirgy) ¢ Dy (r6y), (113) 
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and Dp (r.6,.y) represents the unit vector normal to the boundary at y = Ye, 


Substituting (84) and (96) into (112), substituting (85) and (96) into 


(113), and performing the indicated dot products yields 
Unap(t.r.oy) = lke, { Aop cos nd Yopty) + Bop cos no Y>,(y) 
dake, an 


+ Cop sin nd Y>,(y) + Dap sin nd Y>,(y)J — r 
kK; 
2 





r) 
+ : te cos Np Y>,(y) + Dap cos nd Y>,(y) 
- Any sin ng Y>,(y) - Bop sin no Y>n(v) Jatkr, r¢ 


dY apy ) dY>,(y) 
mn cosi¢ 4 dy + Bop COS nd dy 





sic dY>,(y) 
+ Dap sin ng aera kn nhy y | iat 





+ Coy sin nd 


fy (r.9) ae i .) 
r-— +y 
cr 8h 








or 


Unap(t.r.oy) = E kr, eae cos Nd YY) + Boy cos no Y>,(y) 


dIn(kr, 1) ap (r.¢) 





+ Cop sin no Y>,(y) + Dap sin nd Y>,(y) ae 


d(ky, r) 
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n ' : 
“2 (Co, cos nd Yo(¥) + Dap cos nd Yo,(y) 


of, (re) 
= Adp sin nd Yon (y) = Bop sin nd Y>,(v) Intkr, r) ma 
xy) ey) tal 


ee —— Bop COS Nd —__— dy + Cop Sin nd —{—— dy 


dY>,(y) ei2att 
Bs gigi 





(114) 
By 


and 
Unalt.r.by) = Ike, (A, cos nd Y3(y) + Bz cos ng Y,(y) 


dJa(kr, r) ; 
+ C3 sin nd Y,(y) + D3 sin nd Y;(y)} ee 
d(kr, ft) 


+ 7 C3 cos n¢ Y;(y) + D3 cos né Y,(y) 


- Az sin no Y3(y) - Bg sin no Y;(y)} Jake, £) ¢ 


dY,(y) dY,(y) 
+ Bz cos nd dy 





+ | Az Cos nd 


i dY,(y) 
+ Ds sin nd iy ——} Jn(ke, 0) y | el2ate 





+ C3 sin nd 
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fp (ro) af, (ro) ) 
a ey 











1 
LS 


l 
== 
ie 


Unalt.r.g.y) = ie kr, (A, cos n¢ Y;(y) + Bz cos n¢ Y;(y) 


dake, F) alg (r.4) 
+ Cz sin ng Y;(y) + D3 sin né Y;(y) eee ca 
d(kr, r) : 


- 2 te cos no Y3(y) + Dz cos no Y3(y) 


ofp (ro) 
ad 





- Az sin ng Y,(y) - Bg sin noYz(y)} Jake, r) 


— dY;(y) 
+ Bz cos n¢ dy 





+{A, cos nd 


dY; 3(y) dY,(y) 


el2alt 
+ D3 sin no dy ——} },( In{kr, t) r)| 








+ Cy sin nd (115) 


Np, 


Setting Unarttr.oye ) (114) equal to Unaltr.eye ) (115) and eliminating 


j2aft 
the common term —— yields 


Np, 
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- kp, {Ans COs nd ONCE ? + Bap cos no Yo, (¥p : 


OJalkr, r) af, (ro) 





+ Cop sin no Yoo(¥p ) * Dap sin no Yon(¥p - se 
kr, r) 


n ' - 
2 (ce cos Nd Yop(¥p,) + Day cos nd Yoplyg ) 


af, (ro) 
op 





- Aap sin nd Yon(¥p,) - Bay sin nob Yon(¥p J Jatkr, r) 


dY 2p(¥p,) dY 24(¥p,) dY 2p(¥p,) 


7 ace cos No 


+ Co, sin n 
y 2b > 


+ Bop Cos no yy 


dYop(¥p,) 
+ Dop sin nb = 2 Jntkr, #) 


dy 


=- ke, (A, cos Nd Y3(¥p,) +B; cos no Y3(¥g,) 


Aa(kr, r) fg (r.9) 





+ C3 sin no Y3(¥p,) + D3 sin no Ya(¥p,) dk) 8 
rf 


n : . 
- 2 1G; cos nd Y (Vp) + D3 cos nd Y3(¥p,) 


af, (r¢) 


- Az sin no Y3(¥p,) - B; sin no Yx(¥p,)) Jntkr, r) aera 
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dY,(¥g,) d¥3(¥p,) 


+ (A, cos nd dy + Bz cos nb ay 


d¥3(yp,) d¥(yp,) 
+ C3 sin nd a + D3 sin nm——) Jntkr, 0). (116) 








Factoring (116) yields 


dJalkr, r) af, (ro) dJalkr, ©) of, (9) 


: d(kr, r) ae: : d(kr, r) a 


af (r.4) af, (r.6) 


Nn + Ni = 
: 7-2 In{kr, r) Yoo'¥s, ) = Cob - a Jn(kr, ©) You(¥p, ) eo Dop 


dY 25(¥p,) dV 25 (¥p,) 


+ Jnlkr, 1) Aap + In{kr, 1) Gy Bap! cos no 


dy y 


d]q(kr, 0) af, (r.9) dJy(kr, 1) af, (1.9) 


at YolVp,) a. Cob Ke Sang) 
: d(kr, r) PI ( fo ) Z 


cance iat: 0) 
ar 2b 


fp (ro) ofp (r.9) 


+ 7 Jalkey ) Yoo(¥p,) sp Aaw * 9 Inlkrg 1) Yanl¥p,) gg Bao 
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+ 


dY 26 (¥p, ) dY 25 (Yq) 
* InlKe, 1) “Gy Cap + Jalke, 1) “Gy Dap] sin ng 





| An(kr, r) af (r¢) dJalkr, r) af, (r.4) 
= [= Ke, eee) eee 
2d ks, r) i li : d(k,, r) > 
; af, (1.4) 4 af, (r.6) 
- 72 dnlke, r) Y3(yp,) = C3 - 72 In(Kr, r) Y3(yp,) — D3 
dY3(yp,) dY3(yp,) 
i Jake, r) dy A3+ Jntkr, r) dy Bs | cos nd 
dJn(kr, r) of, (rd) dJalkr, r) ; of (r.4) 
Ky d{k,, r) "al¥p,) is i % et) 3a) “a . 
3 
af, (r.$) af, (1.9) 
Nn + I pal = 1 
+ 72 Jalkr, r) Y3(yp) | Az + 72 In{kr r) Y3(yp,) — B3 
d¥3(yp,) d¥3(¥p,) 
+ Jalkegt) —Go— Cp * Intkegt) —gy— Ds] sin np. (117) 
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Setting the respective coefficients of cos nd and sin nd equal and 


rearranging yields the following: 


dY 25(Vp, ) dY 24(¥p,) 


LJa(ksy£) gy Aan * Jak 1) Bae] 


dy 


dalkr, r) dJalkr, r) fg (rg) 


Sle aero ei fap ke ey 
* lke, 1) * dk, 1) 


afp (rd) 


n se : 
- 7 Late) Yael p,) Cop * date, Youl¥p,) Dao] —— 


dY3(¥p,) d¥3(Yp,) 


a2 Linke, r) dy A; + In{Kr, r) 


y Bs| 





dy 


dIn{kr, r) dJqlkr, r) 
- Ike, ner re ey) A3+kr, euee i B; 


d(kr, r) d kr, r 


af (r.9) 


or 


af, (r.4) 


nN + re 
- 7 [alkegt) Valv¥p,) Ca * Ja(key ©) Y3l¥p,) Ds “age, le 


and 
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dY 25 (¥p,) dYop(Yp,) 
Sang) a, — Coe rer 


dJn(kr, r) dJnlkr, r) af, (r.4) 
a Mamma SAVIOR ee ag — 


kr, r) d{kr, r) oh 


af, (r.4) 


nN + ES 
uF a OCs el You'¥p,) Aap + Jn{kr, r) Yool¥p,) Bop 2, 8 


dY3(yp > d¥,(Yp,) 
= [y(krg Og Oe a(n) eee Da 


dIn(ke, 1) dJalke, 1) af, (r.9) 
- | i ——<--$$- == 1 (yee Ke ee) Dee 
3 d(kr, r) ae 3 d( 7 r) re a 


af, (1.4) 


n + = 
+ [alkrg #) ¥9(¥p,) Ag * Jntkry) Yol¥p,) Bs] yp. 119) 


Again, (118) and (119) are valid only if the associated trigonometric 
function is not identically zero for all values of azimuthal angle $. Conducting 
an analysis of (118) and (119) similar to that of the previous boundary 


condition yields the set of six equations representing the sixth boundary - 
condition (BC # 6) 
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dY 25(¥p,) dY 25(¥p,) 


Ja(Ke, 1) —Go— Aap + Jake, 1) 


dy dy 


BO, (120) 


Yil¥p)) Ane +key Val) Ba 
9 d(kr, r) 


- Jake, 1) Y3(¥g_) C3 - Ja(kr, 0) Ya(vg,) D3 = 0, (122) 


dY oy (Vp > dY 25 (¥p,) 


Jn{Kr, r) dy Cop + In{kr, r) Dop 


dy 


dY3(¥g,) dY3(¥,) 
7 Intkr, r) dy C3 - Intkr, r) dy Dz=0, 


dake, r) dJalkr, r) 
: Yo(¥p ) Cop + Kr You (Vp ) Dap 
: ‘Po c : d(kr, r) 


dake, 1) dJa(ke, 1) 


- Ke, ——— Yal¥g,) Cy - Ke, ——— Yl¥g,) D3 - 0. 


d(k,, r) d ke, r 


and 
In{Kr, r) You(¥p,) Ao + Intkr, r) Yoo(Vp,) Bop 


- Jn(Kr 1) Y3(¥p,) Ag - In(kr, 1) Ya(vg,) Bg = 0. 


af, (re) 





Here, (121) and (124) are valid only if = 


(123) 


(124) 


(125) 


is not identically zero 


for all values of range, r, and azimuthal angle, Similarly, (122) and (125) 
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af, (ro) 
a 





are valid only if is not identically zero for all values of range, r, and 


azimuthal angle, 9. 


The seventh boundary condition is that of continuity of the normal 


component of the acoustic particle velocity at the boundary y = YB, This 


implies 
Unalt.r.oyp,) = Unalt.r.o.¥p,) 
where 
Una(tr.y) = Us(t.r.gy) * De roy). (126) 
Unalt.r.oy) = U(t.r.gy) « fi (Fy) (e277) 


A 
and n, (roy) represents the unit vector normal to the boundary at y = y, . 


Substituting (85) and (98) into (126), substituting (86) and (98) into 
(127), and performing the indicated dot products yields 
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Un3(t.r.@y) = Ike, ON cos nd Y3(y) + Bz cos nd Y,(y) 


dJqlkr, r) n 

+ Cz sin no Y;(y) + D3 sin no Y;(y)} ean 
at) 

3 


~ z {c, cos nd Y3(y) + Dz cos ng Yay) 


- Ag sin np Y3(y) ~ Bg sin ng Y4(y)} Jnl, 0) @ 


dY,(y) dY;,(y) 
+ Bz cos no 








+ LA3 Cos nb dy dy 
dY3(y) dY,(y) 


+ Dz sin nd 








+ Cz sin nd } Jake, ry | eizet 


dy dy 


af, (ro) 7 
¢ te ea [i = r 


or 
Np ; 


afp (to) .) 
a ¢+y 








Unalt.r.oy) - 2 kr, {A, cos nd Y,(y) + Bz cos nd Y,(y) 


d]n(ke, ¢) of, (ro) 





+ C, sin nd Y,(y) + Dy sin no Y;(y) on 
d(kr, r) 


- = ic cos nd Y3(y) + D; cos no Y,(y) 


af, (r¢) 


- Az sin ng Y,(y) - Bg sin ndY;3(y)} Intkr, r) =o 
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+ 


st i 
+ Bz COs no 








¢ he cos No 


+ 


dY3(y) dY,(y) 


el2alt 
2 D3 sin nd dy ——} },( al ky, f) ) 





+ C3 sin nod 





(128) 
B> 


& 


and 


3 + A 
Unglt.r.oy) - Ike, { A, cos nd + B, sin no} VY armas rt 
d(kr, r) 


+ : {B, cos Nd : Ag sin no} Y,(y) Intkr, r) $ 


dYz(y) 


A} j2ait 
dy yle 





+ { A, cos nd + B, sin no} Jake, 7) 





dJn(kr, 7) of, (0.6) 
Unaltirdy) = | - ke, [Ay cos no + By sin no} ¥4(y) ———— 
d(kr, r) 


af ptt) 
ap 





E 2 {B, cos nd - Ag sin no} Y,y) Infkr, r) 


dY,(y) l2att 


Dk cos nd + Bs sin no} Jnlke 0) Gy a (129) 
Bo 


a7 


Setting Unalt.r.o.yp,) (128) equal to Unaltro.Yp,) (129) and eliminating 
ei2alt 


the common term —— reveals 
NB, 


- ks, (A, cos nd Y3(¥g,) + Bz cos n¢ ¥3(Yp,) 


dJnlkr, r) of, (r¢) 
+ C3 sin nd Y3(¥p,) + D3 sin nod Y3(¥p,) emer —— 
( ry r) 


n ' : 
“3 AG: cos nd Y3(yp,) + D3 cos nd Y3(Yp,) 


af, al 9) 


- Az sin ng Y3(¥p, ) - B, sin no Ys(yp, )} Jntkr, .) a6 


dY; (Vp, dY3(yp,) 


+ Bz cos n¢ 








+ (A; cos Nd 


Pe dY3(yp,) 


err Jalkr, F) 


+ C3 sin no—~“qy + Ds sin nd y 


dJn(kr, 1) fp, (r) 


or 





=- kr, ivi: cos ng + By sin no} Ya(¥p,) 
a, f 


f,(r.¢) 


- 5 (By cos nd - A, sin no} we 2) Inlkr, gta x 
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dY,(¥p,) 


+ {A, cos nd + By sin no} Jnlkr, r) (130) 


dy 
Factoring (130) yields 


dJn(kr, ©) af, (rd) dJn(Kr, 0) fp (rd) 
ew on) eaten oie, cena iy) eee 


[- k, 
3 or d(kr, r) 


d(kr, r) 


ofp, (ro) ofp, (r.¢) 


nN + - 
- 27 Jal ke £) V3(¥p9) gg — C3 ~ 22 Jalkeg 2) Yelp.) —y— 


dY3(yp,) d¥3(yp,) 


eee kn 


y Bs| cos nd 


i Jn(Kr, dy 


dJa(kr, 1) af, (ro) dJn(ke, 0) af, (r.9) 
ValVpp) ee C3 - Kr, a. YY) eee ae 


+ [- k 
or 
d(kr, r) dike, r) 


"% 


af, (r.9) af, (r.9) 


+ Jatkeg) Y3(¥p,) —gp— Aa + 22 lrg 1) Y5l¥p,) —o— 


Dy. 


d¥3(yp,) dY3(yp,) 
¥ In{kr, r) dy 33 o In(kr, r) dy D3! sin nd 





a(kr, 0) af, (r.9) 
eT Che 
d(ky, 1) i 


ofp (r.¢) 


Nn + 
" [- kr, Ag - 72 In'kr, r) Y4l¥p,) ve By 


d¥4(¥p,) dJalKr, r) ofp, (ro) 
A,| COS Nd + E kr, “cay — 


a JalKr, r) r 
"4 


dy Ba 


af. {9 dY4(¥p,) 


r = Jalks, r) Y4(¥g,) a Ag+ Ialke, 0) ram sinnd. (131) 


Setting the respective coefficients of cos nd and sin nod equal and 


rearranging yields the following: 


dY3(¥p,) dY(y,,) 
Linke, r) "=a A3 + Jn(Kr, 0) ioe Bs| 


dnlke, 0) da(kr, 0) ofp, (r.¢) 
; Ike, cre 15) A3+ kr, lene B3 


or 
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af p(t) 


n : , 
: 2 Wiles r) Y3(Yp,) C3 + Jntkr, r) Y3(yp,) D; a 


dY4(yp,) dake, 0) af (r.6) 


Ag - ky, ——— Yal¥p,) Ag 


: Jake, 0) 
r) 


dy or 


"4 


af ptt.) 


- dale, 1) Yal¥g,) By ap ee 


and 


atke, r) Sve r 3 Jatkr, . dy 


y Ds 


| dJq(ke, 1) dake, 1) afy,(r.¢) 
- |k,, ———YVil¥p,) Cp + keg ——— Yl Vp.) Ds 
: dike, 1) P2 : dike, 1) P2 ) 


af, {r) 





n ‘ 4 
2 2 atkr, r) Y3(yp,) A3 + Intkr, r) Y3(yp,) B3 op 
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d¥4(¥p,) dJa(kr, 1) af, (r.6) 


cr, 0) 
' ofp, (ro) 
+ 72 Inlkr, 0) Y4l¥p.) Ag ae (133) 


Again, (132) and (133) are valid only if the associated trigonometric 
function is not identically zero for all values of azimuthal angle >. Conducting 
an analysis of (132) and (133) similar to that of the previous two boundary 
conditions yields the set of six equations representing the seventh boundary 
condition (BC # 7) 


dY3(¥p,) d¥,(Yp,) dY,(¥p,) 

Jn(Kr, r) ~ dy A3+ Jntkr, ¢) dy Bs z In{kr, ©) dy Ag=0 , (134) 
dJalkr, r) dJa(kr, 1) 
k Ree itane K Y,(yp,) Bg 
d(ke, 1) ae ; r,t) aii 
alkr, r) 
= Koa ae an, (135) 
4 d(kr, r) a? 
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Jn(ke, £) Y3(¥p,) Cy + Jn(ke, t) Y3(¥p,) Ds - In(keg 1) Ya(¥p,) Bg =-0, (136) 


dY3(yp,) dY3(¥p,) dY4(yp,) 
Jalkr, r) "Ww 8 + Jake, r) uae - In(kr, r) dy Ba = 0, (137) 
qd) (ket) OAK) 
kr, a Y3(¥p,) C3 * ke, 7 Y3(¥p,) D3 
dJalkr, r) 
- ke, too. G B, - 0, (138) 


and 


af, (r.9) 





Here, (135) and (138) are valid only if is not identically zero 


or 
for all values of range, r, and azimuthal angle, ¢. Similarly, (136) and (139) 


ofp, (r.o) 
a 





are valid only if is not identically zero for all values of range, r, and 


azimuthal angle, 9. 
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The final boundary condition is that of discontinuity of the normal 
component of the acoustic particle velocity across the boundary at y = yp. 


This implies 


U ngltirdY¥o) = UL, (t.t6.Vo) + ic cos nd + Gp sin ng} Jn(Ke, €) el2sft (1 40) 


where: 
U_, (try) = Uraltirgy) © Be(rby) (141) 
U a (trey) = Ung(tirey) * Molr.py) (142) 
L(r.oy) - y (143) 


(143) represents the unit normal vector to the boundary at y - y,, and G, 


and G2 represent amounts of discontinuity. The “prime” superscripts (in OP 
and Woe are used to indicate that these velocities are to be evaluated at y = 


y, (This was necessary since the notation U,2, was used in the evaluation of 
the velocity boundary condition at the surface, and the notation U,o, was 
used in the evaluation of the velocity boundary condition at the boundary y 


Substituting (83) and (143) into (141), substituting (84) and (143) into 
(142), and performing the indicated dot products yields 
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Ue Lt oy) Ike, {Ars cos no Y5,(V) + By, cos no Y5,(y) + Coq sin ng Y5,(y) 


: dJq(kr, 0) a 
+ Do, sin nd Y3,(y)) ee We ie. cos n¢ Y,,(y) 


d(kr, r) 


+ Dog cos no Yo,(y) - Aga sin no Y3,(y) - Byg sin noY>,(v)} Jnlkr, re 


dY dY,,(y) 
+ Bog COS No dy 








at Ao, COS Nd 


dY2,(y } dY9,(y) 
¥ j2aft, A 
+ Co, SIN Nd dy + Do, Sin nd dy i lat ket ry elemey 


si dY,,(y) 
+ Bo, COS Nd 








U og tt oy) = Eee cos nd a 


dY},(y) ae 
* Dag sin. nd “Ge See ais rel (144) 





+ Co, Sin no 


and 
U apt oy) - [k, ; Les cos Nd Ve Gy) + Bop COS Nb Y>,(y) 


dJalkr, r) 
+ Cop, sin no Yonty) + Dap sin nd Y>y(y)} a ecenenrs | § 
d(kr, r) 


n . : 
os oe cos Nd Yop(y) + Dap cos nd Yo, ly) 
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- Agy sin nd Y>,(y) - Bap sin no Yop(y)) Jalke, re 


dY>,(y) dY,,(y) 
+ Bop Cos nd dy 





; Kans cos nd dy 


dY>,(y) dY5,(y) 
+ Co, sin nd dy * Dop sin nd 








A : A 
a J} Inte, 0) ¥ | eit. y 


dY>,(y) dY>,(y) 
a Bop cos nd dy 





U op (tit ey) - en cos nd ay 


dY>,(y) dY54(y) | 
+ Dap sin nd ae In(ke, te", (145) 





7 Cop sin nd dy 


Substituting (144) and (145) into (140), dividing out the common 


terms, and evaluating the resulting expressions at y = y, yields 


dY3,(¥o) dY,,(¥0) dY 54 (Vo) dY a4 (Vo) 


{ dy Ar, + dy Boa - dy Adm - dy Bop) cos Nd 
dY54(¥o) dY5,(¥o) dan (va) dY 94 (Vo) 
ik dy Coa + dy Do, - dy Cop - dy DopJ sin nd 
- G; cos nd+ Go sin nd. (146) 


Setting the respective coefficients of cos n@ and sin nd equal and 
rearranging yields the following pair of equations representing the eighth 


and final boundary condition (BC #8): 
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dY5,(Vo) dY,,(y,) dY 54 Vo) Avery a) 
Arat Boa - 











Aap - 


dy dy dy 


and 


dyy,)  ‘dYaely,) d¥5, (ys) dY a (Vo) 


{ dy Co, + dy 


ay Pande Sy. 


Daa gy Con - gy — Da} = Go. 


(147) 


(148) 


Again, (147) and (148) are valid only if the associated trigonometric 


function is not identically zero for all values of azimuthal angle 9. 


F. SUMMARY OF BOUNDARY CONDITION EQUATIONS AND THEIR 


VALIDITY 


To summarize, the boundary condition equations which must be 


satisfied for our general waveguide model are as follows: 
p2l¥s) In(Kr, 1) Yoa(¥s) Ay, + P2l¥s) In(Kr, #) Yoal¥s) Bog 


i py (ys) In{ke, r) Y , (ys) B, - 0, 


which is valid only if cos ng is not identically zero for all values of 9. 


poly¥s) Intkr, r) Y5,(Ys) or + polys) Intkr, r) Y»,(¥5) D,. 


= p1 (Ys) Jalke, r) Y , (ys) A, =0, 


which is valid only if sin no 1s not identically zero for all values of 9. 
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(149) 


(150) 


Yoal¥o) An, + Yoa(¥o) By, - Yop(¥o) Ang - Yon(¥o)By,-0,. (151) 


which is valid only if cos nd is not identically zero for all values of 9. 


Yaal¥o) Coa, Voal¥o) Doa- Von¥o) Co, - Yop(¥o) D,, = 0 1 Mu tls2) 


which is valid only if sin nd is not identically zero for all values of 9. 
e2(Yp,) Jnlkr, r) Yon (¥p,) Ay, + p2lyB,) Jntkr, r) Yoo(yp,) B.., 


- 03(¥B,) Jn(ke, £) Y3(¥p,) Ag - 03(¥B,) Inlkr, 1) Y3(¥p,) Bz =O, (153) 


which is valid only if cos nd is not identically zero for all values of 9. 


é o3(YB,) Jntkr, r) Y3(yp,) Ge e3(¥B,) Intkr, r) Y3(¥p,) D3-0, (154) 


which is valid only if sin ng is not identically zero for all values of 9. 
p3(¥B,) In{Ke, r) Y3(yp,) A, + p3(¥B,) Jntke, r) Y3(yp,) B, 


- 04 ¥Bq) Jake, £) Ya(¥p,) A, 70. (155) 


which is valid only if cos n@ is not identically zero for all values of 9. 
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p3(¥Bp) Inlkr, 0) ¥3(¥B,) C, + 3(¥By) In(Kr, ©) Y3(¥B,) D, 


- p4l¥Bo) Jn(kr, 1) Ya(vp,) B, - 0, (156) 


which is valid only if sin ng is not identically zero for all values of 9. 


dY5,(Ys) dY,,(Ys) dY (ys) 


Jalkr, ¢) oe Aza * Jalkr, 1) ay) bas - Jalkr, r) Page -0, (157) 


which is valid only if cos nd is not identically zero for all values of 9. 


dY5,(ys) dY,(Ys) dY j(ys) 


Ja{Kr, ©) Gy i JnlKr, r) Dog = In{Kr, r) Ay = 0 ' (158) 


dy dy dy 


which is valid only if sin ng is not identically zero for all values of 9. 


dq(kr, r) dq(kr, r) 
Yoal¥s) Aza +t a Yoql¥s) Baa 
d(ke, r) dike, r) 


kr, 


alkr, r) ; 
“Ke CN tg By - 0, (159) 


fg(r.o) 


afe(r, 
which is valid only if cos n¢é and ro are not identically zero for all values 


of r and $. 
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In(Ke, r) Y5,(¥s) Coa a In({Kr, r) Yoa(¥s) Doa = Jn(Ke, r) tin (ys) Ay =Q, (160) 


wae afc(r.o) 
which is valid only if cos nd and ustise are not identically zero for all values 


op 
of r and 9. 


dJqlKr, r) da(Kr, r) 


kp, ———— Yogl¥5) Coa + Kr. ——— Yoogl¥5) Dog 
d(kr, r) d c r) 


dJnlke, r) ; 
- ke, ——— Flys) Ar - 0, (161) 
d(ke, r) 


(ro) 


afer, 
which is valid only if sin nd and oa are not identically zero for all values 


of r and ¢. 


In(Kr r) Y5,(¥s) Ada ate Jn(Kr, r) Y5,(¥s) Bo, = JnlKe, r) Y (ys) By - 0 A (162) 


afo(ro) 
which is valid only if sin ng and ae are not identically zero for all values 


of r and 9. 
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dY54(¥p,) dYo4(¥B,) 


Jn{Kry 1) ge Aap * Jake. 1) “Go Ba 
dY3(yp,) dY,(yp,) 
~ Jatke, 1) Gy Ag - Jake, 1) “Ge Bg - 0, (163) 


which ts valid only if cos ng is not identically zero for all values of 9. 


dYon(¥e, ) dYo4(¥p, ) 


In(Kr, €) Cop + Jn(Kr, £) mar Dop 


dy 


dY3(yp,) dY;(yp,) 


“ Jnlkr, 0) dy 8 = Jafkr, ©) mage 2 =0, (164) 


which is valid only if sin ng is not identically zero for all values of 9. 


Wnlkr, r) alk, r) 
dJalkr, r) alk, r) 
: kp. ———— Y3lyp,) A3 - ke, ——— Y3yp,) B3- 0, (165) 
d{kr, r) d(kr, r) 


7 


ofp (r.) 





which is valid only if cos n@ and are not identically zero for all 


values of r and 9. 
Intkr., r) You(¥p,) Cop + Jntkr, r) You(¥B,) Dap 


- In(Kr, 0) Y3(¥p,) Cy - In(ke, t) Ya(yp,) D3 - 0. (166) 


fg, (ro) 





which is valid only if cos n@é and are not identically zero for all 


values of r and 9. 


Jalkr, F) Jake, 1) 


kr, —. Yon(¥p,) Co, 3 kr, ———* Yoo(¥B,) Dap 
d{kr, r) d(ke, r) 


d)a(kr, 1) dJa(kr, 0) 
- ke, ———_— Yl ¥p,) C3 - Ke ———— Ylyp,)D3-0, (167) 
) d(ke, r) 


ial rf 


af, (r.9) 





which is valid only if sin nd and are not identically zero for all 


values of r and 9. 


Ue 


Jnlke, r) You(¥p, ) Ax» a Jntkr, c) You(¥p,) Bop 


- In(ke, 0) Y3(¥p,) Ag ~ Jn(ke,) Ya(vp,) Bg = 0, (168) 


afg (r.9) 





which is valid only if sin nd and are not identically zero for all 


values of r and 9. 


dY3(Yp,) dY3(yp,) dY;(yp,) 


Jake, r) aye A3+ Jn{Kr, r) dy Bs - Jnlke, 1) ee p A,=0, (169) 


which is valid only if cos nd is not identically zero for all values of 9. 


dY3(yp,) dY3(yp,) dY;(Yp,) 


Jntkr, r) Time: + Jntkr, r) ar a - Intkr, r) ay Bs =0, (170) 


which is valid only if sin nd is not identically zero for all values of 9. 


15 


ky Y3(yp,) A3 + Kr Y(yp,) B3 
* lke, £) * dlke, 1) 
da(Ke, r) 
= ky, ee apes. (171) 
d(kr, r) 
afp.(r.o) 





which is valid only if cos ndé and are not identically zero for all 


or 


values of r and 9. 


Jatkr, r) Y3(yp,) C3 + Jatkr, r) Y3(yp,) D3 ae Jalkr, fy) Y4(Yp,) Bg a 0 ’ (172) 


afp.(r.9) 





which is valid only if cos no and are not identically zero for all 


values of r and 9. 


dJalke, r) dJalkr, r) ; 
Y3(yp,) C3 1 kr, i Y3(yp,) D3 
d r, r 


Kr, 
d(ke, r) 


dJa(Kr ; r) 


SS) imme bi) i) Pele (173) 
d(ke, r) 
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afp(r.>) 





which is valid only if sin nd and are not identically zero for all 


values of r and 9. 


In(kr, 0) Y3(¥B,) Ag * In(ke, 0) Y3(¥By) Bg - In(Ke, 1) Yq(ypy) Ag =0, (174) 


afg.(r.) 
ap 





which is valid only if sin nd and are not identically zero for al! 


values of r and 9. 


BY alyy) dY5,(Vo) dY 54 (Vo) dY 54 (Vo) ; 
dy Atat gy Baan gy Aab- gy Band = Gi. 175) 


which ts valid only if cos ng is not identically zero for al! values of 9. 


dye) dY5,(¥o) dY 54 (Vo) dY5,(¥o) 


era Ca dy Dae-~gy  2-~ gy Dard = G2. (176) 


which is valid only if sin nd is not identically zero for al! values of 9. 
G. DIFFERENCES NOTED DUE TO ARBITRARY BOUNDARY SHAPE 


Before going on to verify that the set of derived general boundary 
condition equations reduces to a well-known and well-documented set of 
boundary condition equations for a very specific set of waveguide conditions, 


the interesting and somewhat unexpected appearance of Jak), i= 1,2,3.4 


75 


terms in some of the general boundary condition equations needs to be 
discussed. These terms cannot be eliminated in the general case of 
arbitrarily shaped boundaries because the radial component of the wave 
number is not constant, and in fact, it depends upon the orientation of the 


local normal vector to the surface. 


To show that this is true, we will begin with the specific case of a planar 
boundary. As shown in Figure 3, the vector wave number k may be resolved 
into its component vectors, k, and ky with respect to the coordinate axes r 


and y, that is, 


ben k, + ky, (177) 





Planar boundary 


y 
Figure 3. Planar Boundary Wave Number Vector Decomposition 
Evaluating the geometry reveals the following pair of equations to describe 


these component vectors: 


k -koos@f= kf, (178) 


and 


ky-ksiney- kyy, (179) 


76 


where k is the magnitude of the vector k, 6 is the angle between the vector 

Te ; ; —. AG 
k and the r-axis, r is the unit vector in the radial direction, y is the unit 
vector in the y direction, k, is the magnitude of the radial component, and ky 


is the magnitude of the depth component. 


Now we will explore the more general case shown in Figure 4. Again the 


k vector may be decomposed into its component vectors with respect to two 





arbitrary boundary 


¥ 


Figure 4. Generalized Boundary Wave Number Vector Decomposition 


very different sets of coordinate axes, the standard r and y axes, and the r’ 
and the y’ axes, which are oriented based on the local norma! vector. When 
decomposed with respect to the r and y axes, the components may be 


expressed as indicated in (178) and (179). However, when decomposed with 
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respect tothe r' and the y’ axes, the resulting components may be expressed 


as follows: 


k.-kcos(o+p)f" -koe , (180) 


and 


ky -ksin(0+B)y' - ky’, (181) 


where k is the magnitude of the vector k, @ is the angle between the vector 


A 
k and the r-axis, B is the angle between the r and r’ axes, r’ is the unit 
A 
vector in the r’ direction, y’ is the unit vector in the y’ direction, k,. is the 
magnitude of the component in the r’ direction, and ky. is the magnitude of 
the component in the y’ direction. 


Using the appropriate trigonometric identities, (180) and (181) may be 


rewritten as 


k..= k (cos 8 cos B - sin @ sin 8) r' (182) 


and 


k,. = k (sin @ cos B + cos @ sin B) y’ (183) 


Equations (182) and (183) may be simplified further by carrying out - 
the indicated multiplications and using (178) and (179), revealing 


78 


k,. = (k, cos B - ky sin B) r’ (184) 


and 


k,. = (k, sin B + ky cos 8) y’ ; (185) 


Equation (184) shows that the component of the wave number along 
the tangent plane (ie., in the r’ direction) at any point depends on the angle 
6, and therefore, on the specific point along the arbitrarily shaped boundary 
at which the vector is to be evaluated. Therefore, this component is not 


constant, and must be maintained in the boundary condition equations. 


H. VERIFICATION OF INITIAL RESULTS 


We will now show that the boundary condition equations derived in 
this section reduce to the well-known and well-documented set of toundary 
condition equations for the following classical waveguide problem: Assume 
that there are only three fluid media, not four, and that 

All sound speeds are constant, ie., 

* c;ly) = Ci. 
* Coly) = Co, and 


id c3(y) Baty. 


All ambient densities are constant, Le., 


*pyly) = py, 
* poly) = p2, and 
° p3ly) = p3. 
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All boundaries are planar and parallel, ie., 
: ys(r.o) i 0, 
*yo(r.o) = Yo, and 
Yp, (0.9) ne D. 


These conditions represent a waveguide made up of three layers. The 
flat boundary at y = 0 separates a semi-infinite medium (medium I: 
- os ys 0) and a finite medium (medium I]: 0 s y s D) of (perhaps) 
different specific acoustic impedances. The flat boundary at y = D separates a 
finite medium (medium I]: 0 s y s D) and a semi-infinite medium (medium 
Ill: Ds ys +) of (perhaps) different specific acoustic impedances. These 
conditions imply that the following arbitrary constants may be set equal to 
zero for the reasons indicated: 
«Bz (no wave reflected in negative y direction) 
¢ D3 (no wave reflected in negative y direction) 
¢ A, (medium not modeled) 
¢B, (medium not modeled) 
Also, n = 0 since plane, parallel boundaries remove angular dependence. 
These conditions also imply that the wave number k and the 
propagation vector component in the y direction ky are constant in a given 
fluid medium. In this case, the solution to (28) is known, and can be written 


as follows: 


Vy), eey SB a (186) 
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Thus, we may set the arbitrary functions ‘a (y) and Y (y) in our 


previous work as follows: 


P -jk 
Y' (y)= Aye hy”, (187) 
and 
Y Gieboe ay. (188) 
af ap, of 
Also, since all of the boundaries are plane, parallel! surfaces, —> sea —S 
) ns) Os) 
of 


By 
and ry are identically zero for all values of r and ». Thus, the following 


boundary condition equations have been invalidated for the reasons 


indicated: 


¢ (150) is invalid because sin nd (for n = 0) is identically zero for all values 
of ¢, 

¢ (152) is invalid because sin nd (for n = 0) is identically zero for all values 
of 4, 

¢ (154) is invalid because sin nd (for n = 0) is identically zero for all values 
of 4. 

¢ (155) is invalid because medium IV is not being modeled, 

¢ (156) is invalid because medium IV is not being modeled, 

¢ (158) is invalid because sin nd (for n = 0) is identically zero for all values 
of $, 


f(r), | 
¢ (159) is invalid because = is identically zero for all values of r and $, 
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is identically zero for all values of r and 4, 


afe(r, 
¢ (161) is invalid because both sin ng and asce are identically zero for 


afc(r) 
¢ (160) is invalid because ae! 


all values of r and 4, 


afs(r.o) 


a are identically zero for 


¢ (162) is invalid because both sin ng and 


all values of r and 4, 
¢ (164) is invalid because sin nd (for n = 0) is identically zero for all values 
of ¢, 


ofp, (r.4) 





¢ (165) is invalid because is identically zero for all values of r and 


ofp, (r.) 





¢ (166) is invalid because is identically zero for all values of r and 


ofp, (1.4) 





¢ (167) is invalid because both sin n@ and are identically zero for 
all values of r and 9, 


ofp, (r.) 





¢ 168) is invalid because both sin no and are identically zero for 


all values of r and 94, 

¢ (169) is invalid because medium IV is not being modeled, 

¢ (170) is invalid because medium IV is not being modeled, 

¢ (171) is invalid because medium IV is not being modeled, 

* (172) is invalid because medium IV is not being modeled, 

¢ (173) is invalid because medium IV is not being modeled, 

¢ (174) is invalid because medium IV is not being modeled, and 

¢ (176) is invalid because sin nd (for n = 0) is identically zero for all values 
of 9. 
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Thus, the original set of 28 equations in 17 unknowns has been reduced 
to a set of six equations in six unknowns. The next step will be to evaluate 
each of the remaining equations in turn so that they may be compared with 
the equations developed by Ziomek (1991) for this particular waveguide 


problem. 


The first of the remaining equations is 
p2lys) Jn(Kr, r) Y5_(¥s) ee e p2lys) In{kr, i) Yo,(¥s) B,. 


- pi(¥s) Jn(ke, ) Yy(ys) By = 0. (149) 


The Bessel function dependence of (149) may be eliminated by virtue 
of the fact that, in this simple waveguide problem, the radial component of 
the propagation vector is the same in all three media (implying that the 
Bessel functions may just be divided out). Recalling that the densities are 


constants and that the value of y at ys is identically zero, (149) becomes 


p2 Y>,(0) Aj, + P2 Y2,(0) B,, - py ¥,(0) By = 0. (189) 


Substituting the appropriate functions of y into the simplified 


expression (189), and noting that el? is equal to unity yields the following: 


p2 A, + 92 B,, - 9, By = 0. (190) 
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This equation is the same as that derived by Ziomek (1991, equation 


(3.9-25)) for this boundary condition. 


Conducting a similar analysis on (151) yields 


This equation is the same as that derived by Ziomek (1991, equation 
(3.9-34)) for this boundary condition. 
Equation (153) reduces to the following (after additionally noting that 


B; has been set equal to zero) 


-jk. D tjk_ D -jk_ D 
p2e 19 Ay, + p2e LS B,,  P3 & Y3 A3=0. (192) 


This equation is the same as that derived by Ziomek (1991, equation 


(3.9-44)) for this boundary condition. 


Equation (157) becomes 


KY Ara ~ ky, Baa+ ky By - 0. (193) 


This equation is the same as that derived by Ziomek (1991, equation ~ 


(3.9-30)) for this boundary condition. 
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Equation (163) reduces to the following (after additionally noting that 


Bz has been set equal to zero): 


-jk, D +jk D -jk D 
ke %2 Any - ky, %2 Boy-k,e %3 Ag-0. (194) 


This equation is the same as that derived by Ziomek (1991, equation 


(3.9-49)) for this boundary condition. 


ane 


5 (as suggested by Ziomek (1991, 


In order to evaluate (175), let G; - 


equation (3.9-35) and following)). Substituting yields 


“ik y tik, oy 
: y, 0 ; y,°0 
= Ko foe Ag | KY e “2” Boa 


(195) 


This reduces to 


ike. ¥ +ik oy -{k oy +jk Oy 
=—1e %2 ° Aoate 73 ? 


85 


This equation is the same as that derived by Ziomek (1991, equation 


(3.9-38)) for this boundary condition. 


Thus, we have shown that the theoretically derived set of equations for 
a general waveguide problem reduces to the set of equations expected for 
the classical waveguide problem. This provides us with the confidence to go 


on with the solution for the unknown arbitrary constants. 
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IV. SOLUTION FOR THE UNKNOWN ARBITRARY CONSTANTS USING 
SYMBOLIC ALGEBRA CAPABILITIES OF Mathematica 


In Section III of this thesis, we developed a set of 28 boundary 
condition equations in the !7 unknown constants. The purpose of this section 
is to generate a solution to this system of equations for the general 


Waveguide case. 


A review of (149) through (176) reveals that the coefficients of these 
unknown constants are, in general, complicated expressions involving depth- 


dependent densities, range-dependent n'" order Bessel functions, and as yet 


unspecified depth-dependent velocity potential functions (i.e., the ee and Y, 


terms). In order to maintain the generality of the generated solution, we will 
require either many long hours of tedious algebra involving manipulations of 
these complicated expressions (with the high probability of algebraic errors) 
or a computer program capable of conducting such manipulations directly on 
these symbolic expressions. Fortunately, Mathematica for the Macintosh 
computer (version 1.2.1 [33 (enhanced)) is the one such program available to 
us at the Naval Postgraduate School, and therefore will be used to generate 


the general solution desired. 


The first step in this process will be to program //athematica to solve 
for the unknown constants for a very specific set of waveguide conditions. 
By doing this, we will gain experience in using the program and confidence 


that the program output is reliable. For this work, we will use a three media 


87 


Waveguide with plane, parallel! boundaries. Using vector-matrix notation, a 


compact system equation may be written for this (or for any other) case as 


follows: 


Ax-b, 


(198) 


where A is the matrix of coefficients, x is the column vector of unknown 


constants, and b is the column vector of known constants. 


For the three media waveguide with plane, parallel boundaries, these 


vector-matrix quantities are defined as follows: 
A= Tai at2 a3 00 
0 422 423 424 a25 0 
0 0 QO 434 435 a36 
a41 a42 443 0 O 9 


0 0 0 a54 855 a56 
10 862 463 464 a65 0 


where 


ay =~ pyl¥s) YyCys) 

12 - Palys) Yoav) 

1,3 = pal¥s) Yoalys) 
apo = YoalVo) 
a23- Yoal¥o) 
ao4=- eelye) 


a2 5 = Yon(Vo) 


88 


(199) 
(200) 
(201) 
(202) 
(203) 
(204) 


(205) 


a3.4 = p2l¥p,) Yov(¥p,) 
a35 = pal¥p,) Yoo(¥p,) 


43,6 ~~ p3(¥p,) Y3(¥p,) 

















dY (ys) 
44.1 dy 
dY,,(ys) 
a4i2 = dy 
dY>,(Ys) 
avons > 
a54 7 dy 
dV op(¥p,) 
cima ae 
dY3(¥p) 
a567 7 dy 
dY5,(¥o) 
46,2 dy 
dY,, (Vo) 
a6 3 * dy 
dY 54 (Vo) 
a647 7 dy 
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(206) 


(207) 


(208) 


(209) 


(210) 


(211) 


(212) 


(213) 


(214) 


(215) 


(216) 


(217) 


dY 54 (Vo) 


a65-- Gy (218) 


x= [B, Ava Boa Aap Bop Az)" (219) 


where the superscript T indicates the transpose matrix operator (indicating 


that x is a column vector), 


b-(00000G,}" (220) 
where the superscript T indicates that b is a column vector, and 


G, =- se, (221) 

It should be noted here that we have defined the matrix A in a very 
specific manner. Each row of A represents one of the valid boundary 
condition equations for the specific waveguide being studied. These appear 
in the order presented in Section III. For the three media waveguide with 
plane, parallel boundaries, row 1 of A contains the coefficients found in 
(149). Row 2 contains the coefficients found in (151), and so on. Row 6 
contains the coefficients found in (175). For simplicity'’s sake, we have used 
generic elements, such as a; 9, to replace the more complicated expressions. 
In addition, zeros have been used to indicate that the appropriate unknown 
constants do not appear in a specific boundary condition equation. We will 


use this convention in the work which follows. 
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In this three media waveguide case, A is a six by six square matrix. 


Therefore, the solution to (198) may be written directly as 
x-A'‘b. (222) 


where the superscript -1 indicates the inverse matrix operator. 


Using Wolfram (1988) as a programming reference guide, a 
Mathematica “notebook” was created to solve this three media waveguide 
case using the solution technique expressed in (222). The A/athematica code 
required to perform this task is as follows: 

a~=i{(alcl, alc2, alc3, 0, 0, 0), 

(0, a2c2, a2c3, a2c4, a2c5, 0}, 
(0, 0, 0, a3c4, a3c5, a3c6}, 
{a4cl, a4c2, a4c3, 0, 0, 0), 


(0, 0, 0, a5c4, a5c5, a5c6), 
(0, a6c2, a6c3, a6c4, a6cS, 03}: 


b = (0, 0, 0, 0, 0, G1); 


x = (Inversefa]).b 


In developing this code, we continued to utilize the generic matrix 
elements described earlier. Two subtle differences in the notation used in 
the code from the notation discussed earlier need to be pointed out. First, we 
have represented the elements of matrix A (for example, a, 2) as individual 
variables (the corresponding variable name would be alc2). This slight 
deviation in notation was used because subscripting as defined by 
Mathematica would not have been useful for our purposes. In this revised 


notation, the small case letter ‘c’ represents the comma in the element name. 
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This deviation was required because the program uses the comma to 
separate individual array elements. The second notational comment refers to 
the fact that a lower case letter "a" was used to represent the matrix A. This 
was required to conform with Mathematica 's notational convention, which 


reserves names beginning with capital letters for built-in functions. 


Unfortunately, running this code resulted in halted execution due to a 
singularity error. We surmise that the problem occurred when the program 
was attempting to take the inverse of the matrix A. Luckily, Mathematica 
has a built-in function, LinearSolve, which evaluates (198) directly if the 


matrix A is a square matrix. Thus, the following revised code was written: 


a= ({alcl, alc2, alc3, 0, 0, 0}, 
(0, a2c2, a2c3, a2c4, a2c5, 0}, 
{0, 0, 0, a3c4, a3c5, a3c6}, 
(a4cl1, a4c2, a4c3, 0, 0, 0}, 
{0, 0, 0, a5c4, a5c5, a5c6}, 
{0, a6c2, a6c3, a6c4, a6c}, 03): 


b = (0, 0, 0, 0, 0, G1}; 


LinearSolvela,b] 


The revised code ran successfully. The output of this code is the desired 
vector x. When Mathematica functions such as Factor, Cancel, and Simplify 
were applied to the output, the same result was returned, indicating that the 
program was satisfied that the output was as simple as it could make it. 
Closer inspection of the output revealed that each of the six elements was of 


the form 


SZ 


num, 
ie ee 
1” denom 


(223) 
where x, represents the first element of the vector x (in the three media 
waveguide with plane, paralle! boundaries, this element is the unknown 
constant B,), num, represents the numerator expression for the first 
element, and denom represents the denominator. Fortunately, all of the 
elements of the output vector have a common denominator. This inspection 
also revealed that some algebraic manipulations could be manually 
performed to simplify the expressions somewhat. Thus, the robustness of the 


symbolic algebra functions of Mathematica is at best questionable. 


We will now present the results of the program for the three media 
waveguide with plane, parallel boundaries. The first step will be to simplify 
the results manually in order to generate generic expressions for the 
unknown constants in terms of the generic elements. Second, we will 
substitute (199) through (218) and (221) into the generic expressions to 
reveal general expressions for these unknown constants. Finally, we will 
assume constant speed of sound and constant density and show that the 
general expressions formed from the Mathematica output are the same as 
those derived by Ziomek (1991) for the classical waveguide case. This 
verification will be conducted in the following order: Ao, Bog, Aap, Bap, By, 
and A3. We will demonstrate this entire process for the unknown constant 
Ao, only, and simply present the results for the other five unknown 


constants. 


oe 


The first output expression to be explored will be the common 


denominator, denom 


denom = - a1 3 425 43.6 441 35.4 462+ 811 42,5 43.6 a4 3 a5 4 46 2 
+ Aj 3 82 4 83.6 Aq) a5 5 46.2 - Ay 22.4 43.6 24 3 a5 5 AG 2 
+ 413 425 23.4 841 25,6 AG 2 ~— 21,3 A2,4 Az 5 Aq a5 6 AG 2 
~ Ay 1 42,5 43.4 443 a5 6 262 + Ay) 22,4 AZ,5 Aq 3 AS 6 a6 2 
+ 412 42,5 43.6 Aq.1 45,4 46.3 ~ Ay 1 225 43.6 a4 2 a5 4 463 
~ 84,2 42,4 43,6 44,1 45,5 46,3 + Ay 1 82.4 43.6 442 a5.5 363 
~ 41,2 42,5 43,4 441 5,6 26,3 + 212 82.4 az 5 ag a5 6 a6 3 
+ Ay 1 42,5 3,4 Aq 2 a5 6 46.3 - Ay) 22,4 AZ 5 Aq 2 a5 6 AG 3 
~ 813 82 2 43.6 41 5,5 46.4 + Ay 2 B23 43.6 ag) 5,5 464 
~ 81,1 82,3 36 A42 45,5 46.4 + Ay 1 B22 43.6 a4 3 55 AG 4 
+ 413 42.2 43.5 a4 1 25 6 26,4 ~ 81,2 42,3 AZ 5 Ag 35 6 464 
+ 1,1 42,3 43.5 84.2 a5.6 6.4 - Ay 42,2 AZ 5 a4 3 a5 6 464 
+ Ay 3 42,2 43,6 Aq.) 45.4 26.5 ~ 41,2 42,3 az 6 ag a5 4 Ags 
* 1.1 82,3 43,6 44.2 45,4 465 - Ay 82,2 43,6 Ag 3 a5 4 G5 
~ 81,3 42,2 43.4 41 45,6 46,5 * 81,2 42,3 A344 A5.6 365 
- Ay 1 82.3 434 49 a5 6 G5 + Ay, B22 a3 4 a4 3 35.6 G5. (224) 


Factoring (224) reveals 


denom - a 1 425 443 46,2 (43.6 a5 4 - 43,4 a5, 6) 
- 813 425 a4 1 46,2 (23.6 5 4 - 3.4 5.6) 
+ 413 82.4 441 26,2 (23.6 45,5 - 23,5 a5 6) 
- Aj 1 2,4 44,3 262 (43.6 a55 - a3,5 a5 6) 
* 1,2 42.5 a4 1 46,3 (43.6 5,4 - 23.4 45,6) 
- 11 42,5 44,2 26,3 (43.6 45,4 - 23,4 a5 6) 
+ 1,1 49,4 4.2 463 (43.6 as.5 - az 5 a5 6) 
- 81,2 42.4 Aq 1 26,3 (23,6 a5,5 - Az5 35,6) 
+ M12 42,3 44,1 26,4 (43.6 a5.5 - a35 a5 6) 
- 13 42,2 a4 1 6,4 (23.6 25,5 - 23,5 5,6) 
+ 1,1 42,2 44,3 46,4 (43.6 a55 - 435 as 6) 
- Ay 1 82,3 24,2 26,4 (23.6 a5.5 - 33.5 a5 6) 
+ 13 2,2 841 265 [43.6 5.4 - 3.4 25,6) 
- Ay ,2 42,3 841 Ag5 (43.6 25,4 - 3,4 35,6) 
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* 14 42,3 4,2 46,5 (23,6 5.4 - 43,4 5,6) 
- 41,1 42,2 44,3 26,5 (43,6 45,4 - 43,4 5,6) . (225) 


Collecting common terms yields the following generic expression for the 


common denominator : 


denom = {a3¢ as 4 - 434 aso) lay a2 § 443 462 
~ 4,342 5 441 AG 2+ 21,2 425 441 463 - a1 42,5 44,2 463 
+ Ay 3422 441 465 ~ a1, 23 44 465 + Ay 403 ag 2 AGS 
~ Ay 1 42,2 a4 3 465) + 143.6 45.5 - agg as ¢} a, 3424 a4 1 462 
~— 81,1 42,4 a4 3 ago + Ay Ad 4 Aq 2 AG 3 ~ Ay 9 04 Ag AG 3 
+ 12 423 441 464 - 213 422 444 AG 4 + Aly A22 a4 3 864 
= ai ar 3 a4 2 a6 4 : (226) 


Substituting (199) through (218) into (226), and using the facts that 
ys = 0 and vee D reveals 


dY3(D) dY,,(D) 


= p3\D) eee Y3(D)) 


denom = (p2(D) 5,00) i 


dy 


dY3,(¥o) dY,,(0) dY {(0) dY3,(¥o) 


ona Yop(Yo) - P2(0) — narra W101) Yon t ye! 


x [p4(0) ¥,(0) ; 


dy 


dY (0) dY3,(¥0) dY5,(0) dY5,(¥o) 


+ p2(0) Go Yoa(0) —Go— Yan vo) - p,(0) Y,(0) Vulyo! 


dy dy 
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dY,(0) dY a4 (¥o) dY ,(0) dY54(Vo) 


+ p3(0) == Vanly,) Ying (0) dy p2(0) a T4(0) Youve) 


dy 
aah dY54(¥o) dY,,(0) dY>, (yo) 


+194 (0) (0) aeereemeniie (vie) pare i) NAO) “ty a 


dY,(D) dY5,(D) 


+ {p9(D) ¥>,(D) - p3(D) —— ¥3(D) 


dy 


dY (0) dY5,(y,) dY5,(¥o) dY,(0) 


x p2(0) Y,,(0) None) Je p; (0) Y,(0) dy Yon(Vo) 


dy dy | dy 
dY5,(0) dY5, (yo) dY (0) dY>,(¥o) 
+ p(0) Y,(0) Te er, Dan 650) aoreem Yee(O) ay YaplVo) 


dY ,(0) dY a4 (Vo) dY ,(0) AY 54,( Yo) 


Y5,(0) Yoq(¥o) - 92(0) —F— Vaalv¥o) Yoa(0) —— 


- p2(0) ——— dy 


dy dy 
dY5,(0) dY a4 (Vo) = AY 54 (Vo) 


p,(0) Y, (Qe Yo4\¥o) =a 


+ py(0) Y,(0) Y5,(yp) ae 


Tar 
(227) 


Simplifying (227) reveals 
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dY5-(y,) d¥,,(0) dY¥,.(0) d¥,(y,) 
denom = [p,(0) Y ,(0) 


dy dy dy) dy 
dY}(0) dY24(¥o) dY>,(Y0) 
- p2(0) =. dy Y2al0) - Y3,(0) ar) 
dY,(D) dY5,(D) 
x LYon(Vo) {o(D) VAD) ae p3(D) Scie Y3(D) 
dY3(D) dY>,(D) 
- ¥3u(vo) (o2(D) ¥oq(D) —— - px(D) —— YD) 
dY¥2,(0) dY2,(0) _ 
- [p,(0) Y,(0) Yoal¥o) ar an a OTOL 


dY,(0) 


00) ae 


a (pale) You(0) ~ ¥5q(0) Yoal¥o) | 


is Yo) dY,(D) dY>,(D) 


x |—— {9(D) D) Y3,(D) dy 7 


—— yi(p)}]. (228) 
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Equation (228) is the final genera/form of the Mathematica output 
common denominator. \f we now assume constant speed of sound and 
constant density in a specific medium, the depth-dependent functions in 
(228) become complex exponentials and the denominator becomes ( using the 


subscript c to indicate the constant speed of sound assumption) 


-jk OD 
denom,=--j2k. e % 
M5, 
{( V( ee 
x Up, ky, - p2 ky P2 Sp ~ p3 Ay, e 2 
r None | 
+\p, es, + p2 Ky P2 oy + p3 Ky, e 2 , (229) 


Now that the denominator has been simplified, well concentrate on 
obtaining expressions for each of the unknown constants in the order stated 


above. The first constant is 


num 
2a denom ' (230) 
where 
NUM = G; |- a1 .3 A25 A364 aq) a5 4 + Ay 1 AQ5 AZ6 AG 3 AS 4 
+ Ay 3 82.4 23.6 Aq 45.5 ~— 211 82.4 43.6 a4 3 a55 
+ 13.425 a3 4 441 45.6 - 41.3 42.4 35 Ag O56 
~ 1,1 825 43.4 4.3 45.6 + Ay 1 24 A3,5 44,3 35,6) - (231) 
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Factoring (231) and collecting common terms yields the following 


generic expression for the numerator of Ad! 


NUM» = G; lay. 25 443 - 44,3 42,5 4,1) (23,6 45,4 - 83,4 5,6) 
+ (ay 1 82,4 4,3 - 41,3 82,4 84,1) (23,5 25,6 - a3,6.25,5)] . (232) 


Substituting appropriate expressions into (232) yields the following 


general erypressionfor the numerator of Aga: 


ne dY,,(0) dY ,(0) 
num, = 5 (p4(0) ¥;(0) —— - 210) 5 Y;,(0)) 
dY3(D) dY5,(D) 
x [¥ag(v0) Lo2(D) ¥34(D) —— - p(D) YD) 
dY,(D) dY5,(D) 
- Yu(¥o) {oo(D) Ypg(D) ——- ps(D) ——¥3(D)}]. (233) 


Again, making the constant speed of sound and constant density 
assumptions and substituting the appropriate depth-dependent expressions 
allows us to write the numerator of A», as (using the subscript ¢ to indicate 


the constant speed of sound assumption) 
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Thus, the gerera/ result for Aza is 


dY,,(0) dY ,(0) 


-k s 2 
Aras 35 (0110) ¥4(0) —G— (0) Ge 


dy 


dY3(D) dY>,(D) 


x [Y5p(¥0) (o2(D) ¥3,(D) aes 


dY,(D) dY,,(D) 


- ¥34(Vo) {09(D) Y>q(D) apes) rae yx(D)}] / 


AY 53(¥o) dY9,(0) dY5,(0) dYo,(¥o) 
[o,(0) Y ,(0) ‘a 











dy dy dy 


dY,(0) dY>,(yo) dY5,(Yo) 


Ve 10) =e (Ol 


dy dy 
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dY,(D) dY,,(D) 


x [¥zpl¥o) LentD) ¥o4(D) a - ex(D) —— ¥3(0)} 
dY,(D) dY>,(D) 
dY2,(0) dY3,(0) 
- [p,(0) Y,(0) (Y3,Cy9) a a 
dY (0) 
- p20) —— (Valve) Yoal0) - ¥34(0) Yoal¥o)) | 
= Yo) dY,(D) dY5,(D) 
»{(——— { p(D) D) Y5,(D) ae p3(D) dy Y3(D) 
AY 54(Vo) —— d¥,(D) dY5,(D) 
ar {p2(D) Yop(D) dy p3(D) la y3(D)}] 


For constant speed of sound and ambient density, A, reduces to 
a 
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jk y  -jk. D jk y +jk. OD 
x [los k, -p2 ky emu, (een a + (p2 k, + p3 ky Je 72" 6 ae | / 
2 3 3 2 
My” {( )( ee 
bikers 3 py Ky poy, pak, pak, e 
) en (236) 
+lo, K +pok po k +p3;k Je “2 , 23 
{ ¥ 2 3s ¥3 3 Y5 


We must now verify that (236) simplifies to the well-known and well- 


documented expression for this unknown constant. Eliminating the common 


-jk D 
terme 3 and dividing numerator and denominator by (p \ ky + P2 ky ) 
2 i 


yields 
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(py k 


Yo ic ie 


Ada ( ky 2 k, Je 
4n kk (p, kK +p.k ) 
Yt OY) rf 


If we define a reflection coefficient at the boundary between medium 
two and medium one as suggested by Ziomek (1991, equation (3.9-54}) as 


follows: 


Ro} < (238) 


then, (237) becomes 


+)k be “ik, D 


jk 
Ag Ra Mosk, - wrk, De 90 ee OD ms 
4x k 
% 
-jk yy +jk D 
(po k +p3k Je "2°e %2 ted 


“ity D rik, D 
[Re (pak, -p3k, Je "2 +{ogk, *o3k, Je "2 J. (239) 
a 2 3 2 


Rearranging the denominator and dividing both numerator and 


denominator of (239) by (p. k, + 93 ky ) reveals 
3 2 





+jk, y -jk, D -jk, y, +k, D 
rs OO) hes ee ae 6D | 
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+jk D D 
le %2  -R,,——————e "2 ]. (240) 


If we define a reflection coefficient at the boundary between medium 
two and medium three as suggested by Ziomek (1991, equation (3.9-55)) as 


follows: 


Ry3 = —————_ (241) 


then (240) becomes 


| k ky oy, -ik, D -iky y, tik, =D 
Ay, =~ Ray [Rg € te 2 ee 2 ee 7 


+jk OD 
(eae) Ree a oh (242) 
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+jk D 
Dividing the numerator and denominator of (242) bye 2 


reveals 
sik oy. -j2k oD -jk oy 
ip ee, Ro3 e OL GE TET yp, ol / 
4n k 
Y2 
-j2k D 
1 - Roy R23 e 2 (243) 





-jk, y 

Factoring e 2 ° out of the numerator and rearranging yields 

ik +j2k. y  -j2k. OD -jk sy 

Aga = i Ro 1+R3e are uy le Mao, 

4n k 
¥9 
-j2k D 
Beanies 2 I. (244) 


Rearranging (244) provides us with the following desired expression: 
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-j2k, (D a) EA, 


[1+ Roe %2 e %2 
_ iat aes 
te ee (245) 
-j2k. D k 
[1 - Ro; Rog e se 2 


This is exactly the equation derived by Ziomek (1991, equation (3.9-50)) for 


this unknown constant for the classical waveguide case. 


For the unknown constant B»,, Afathematica provided the following 


numerator: 


mums = Gy |ay.2 42.5 43,6 441 45,4 - 44,1 42,5 43,6 a4,2 5,4 
~ 41,2 42,4 43.6 44.1 455+ 411 22,4 A3.6 a4 2 a5 5 
~ 84 2 42,5 43.4 ag} 45.6 + 41,2 42,4 a3 5 a4) A564 


* 411 425 83.4 4.2 45,6 - ay 1 2.4 435 42 35,6) - (246) 


Factoring (246) and collecting common terms yields the following 


generic expresstonfor the numerator of Bos: 


num, = Gy lays 25 84.2 - 21,2 82,5 4,1) (23,4 25,6 - 83,6 25, 4) 
: (ay 42.4 44.2 ~ 41,2 42,4 aq) (4364535 - 435 4a56)]}. (247) 


Substituting appropriate expressions into (247) and using (228) yields 


the following genera/erpression for By, 
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Bra 5 (p px(0) ¥;(0) 5 — pal0) —— ¥},(0) 
dY,(D) 
x [Yps(¥0) {e3(D) —— Y}(D) - p(D) Yo4(D) 
dY3(D) aa 
+ You(¥o) {p2(D) Yop(D) Gy - e3(D De -— a(D) I} / 


dY5.(¥o) dY2,(0) dY2,(0) dY.,(y) 
dy dy) ed dy 
qY (0), dite) dY,,(¥ 0) 


- p2(0) dy ay Y2al0) - Yy,(0) ae) 


[p,(0) Y,(0) 


= sel 
x [¥jp(¥o) Lp2(D) ¥3,(D) —— - p3(D) —— Y3(0) 


a dY5,(D) 
- Yu(¥o) {o2(D) ¥pg(D) Ge - pa(D) —— Y3t0)}] 


dY2,(0) d¥2,(0) 
- [p,(0) Y,(0) (¥3.(y0) = “ay Y2al¥o! 
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= GAD) aren 


Ene, dY,(D) dY5,(D) 


Fy boot) Yon) ——- est) —— YH]. (248) 


For constant speed of sound and ambient density, (248) becomes 


ky eye 
Boa On (py Ky, + po ye 3 


“ky y, -jk, D -jk, y, +jk, D 
xilteak -pok ewe? 6 wo +(p.k + 93k ee 2 ee | / 
Yo Y, ome Yo 


-jk. D 


-jk_ OD 
y : meeane y 
Aan ky, p2ky Mork, esky Je %2 


-j2k_ e 
Y9 


+jk OD 
Corky +e2k, York, sesk, De Y. } (249) 
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Using the definitions of Rp; and Ro3 presented earlier, (249) may be 


reduced to 


-j2k (D -y ), -jk, y 
[1 + Rog e %2 le eye: 





vas 
Boa 4n 


F - Roy Ro3 5 tne 2 


ee Roy ; 
Multiplying (250) by Ro B,, may be written as 
Zio aec 


l 


Boa eRe aca 


k,. 


(250) 


(251) 


This is exactly the equation derived by Ziomek (1991, equation (3.9-51)) for 


this unknown constant for the classical waveguide case. 


For the unknown constant A», Mathematica provided the following 


numerator: 


numg = Gy [- ay,3 22,2 03,6 24, 45,5 * 41,2 42,3 43,6 24.1 45,5 


~ Ay 1 42,3 43 6 Ag. 2 A565 + Ay 1 AD 2 436 Ag 3 A555 
+ A13 422 A355 aq 45.6 ~ a2 A2,3 435 ag 45.6 


+ Ay 1 82 3 835 Ag 2 agg - Ay | AD 2 AZ 5 Ag 3 ASG) . 


(252) 


Factoring (252) and collecting common terms yields the following 


generic expressionfor the numerator of Aop: 
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numg = G; lay, 42.2 444 - Ay,2 423 4,1 + Ay, 22,3 42 - 
ait a2 2 a4 31 (az 5 as 6 = az 6 as 5) Fe (253) 


Substituting appropriate expressions into (253) and using (228) yields 


the following genera/ expression for Aop: 


dY ,(0) dY,,(0) 


-k 2 ; 
Am" 35 LYz4lve) (e910) dy Y2al0) - 94(0) Y,(0) —— 


dY ,(0) dY,,(0) 


“Yoal¥o) (ox(0) dy Y2al0) - 9410) ¥4(0) ~— 


dY,,(D) dY,(D) 


Se aD ae (Ol 


dY5,(¥o) dY2,(0) dY,,(0) dY,,(¥o) 
{(0) dy dy dy dy 











[o,(0) Y 


dv (0) dv, (y,) GNA) 
- p2(0) Syn agg 1289) - Y,,(0) dy 





Lit 


dY,(D) dY5,(D) 
x [Yp(ve) Co2(D) Yp4(0) —— - 05D) — — ¥3(0)) 


dY;(D) dY>,(D) 


- You(ve) Lop(D) ¥py(D) ——- e(0) —— Y3t0)}] 


dY5,(0) d¥,,(0) 


- [p4(0) ¥,(0) (y3,(y,) Y>_(Vo)) 


mas oa 


dY,(0) 
- p2(0) (Y3,( 79) Yoq(0) - ¥3_(0) Yoal¥0)) | 


= Yo) si ik 
x | — {p,(D) ¥},(D) —— - p3(D) —— ¥3(D) 


dY 54 (¥o) dY,(D) dY5,(D) 


— (o9(D) Yoq(D) -p,(D) ——YiD)}]. (254) 


dy dy 


For constant speed of sound and ambient density, (254) becomes 


“ik, se tik, D 
to; k +pok Je “2° e 72 
as, us 
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+ (py ky Pa ky Je i ie Iloaky + ak, We 
eee Eek Conk amen eran 
= Cys = : G 
J Y Pr ¥» P2 yy P2 , P3 Y5 
tik, D 
q (p, sp + (5 ay ) (p> ky, uae or k, Je 2 } . (255) 


Using the definitions of R2; and R23 presented earlier, (255) may be 


reduced to 


Piivceensee se See | 
Ary.” 4x | (256) 


-j2k Dyk 
1- Ry Rye "2 | Y 


This is exactly the equation derived by Ziomek (1991, equation (3.9-52)) for 


this unknown constant for the classical waveguide case. 


Eis 


For the unknown constant Bz, Mathematica provided the following 


numerator: 


nums = G, la; 3 a22 43,6 44,1 45,4 - 31,2 22,3 43,6 24,1 95,4 
+ Ay 1 42,3 43.6 442 45.4 ~ Ay 4 42,2 43.6 443 a5 4 
~ 81,3 82,2 43.4 44.1 45.6 + 81 2 42,3 a3 4 a4) 25,6 
- Ay, 82,3 3,4 24,2 85,6 * 21,1 22,2 43,4 84,3 5,6) . (257) 


Factoring (257) and collecting common terms yields the following 


generic expressionfor the numerator of Boy: 


nums = G, la; 3 42.2 444 - 4,2 42,3 44,1 + 84,1 42,3 A4,2 - 
1,1 42,2 24,3] (43,6 45,4 - 23,4 45,6) . (258) 


Substituting appropriate expressions into (258) and using (228) yields 


the following genera/ expression for Boy: 


dY (0) dY,(0) 


ety ; 
Bap 5 [Yaa Ve) (p2(0) —G— ¥34(0) - 91(0) ¥4(0) 


dy dy 


dY ,(0) dY,,(0) 


- Y5,(¥o) \p2(0) ——— Y,,(0) - py(0) Y,(0) 


dy dy 


dY;(D) dY 94 (D) 


- p3(D) ———;(0)} / 


x (po(D) ¥3,(D) 7 


dy 
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dY34(Vo) dY>,(0) dY5,(0) d¥3,(vo) 


[p1(0) ¥,(0) i ae = a 
d¥,(0) d¥3,(¥o) dYs.(y5) 
- p2(0) ray ay Y2al0) =F pal) ie 
dY;(D) dY>,(D) 
x L¥zulv0) Lo2(D) ¥39(D) —— - 03D) —— YG) 
dY3(D) dY9,(D) 
- Y3o(¥q) (p2(D) Y¥>p(D) ay - st) a YG) | 
dY,,(0) dY,,(0) 
- [p,(0) Y,(0) (¥3,Cy9) Ady ee ae) 


dY ,(0) 
- 020) se (¥;a(v0) Yoa(0) ~ Ypa(0) Yoa(vo)) 


dY a4 (YQ) dY3(D) dY5,(D) 
« LF Co(D) ¥34(D) —— - 03(D) —— YG) 
dY (ys) dY,(D) dY,,(D) 
-—F— Coal) 540) — - 030) Ge — ¥300)) 


15 


(259) 


For constant speed of sound and ambient density, (259) becomes 


k. -ik, D ay > 
aces 2 Ito, k 


+k a 
->e 
2be 2x 


pez.) Yo 


-jk, y 
+ - ia - 
(py Ky Pa ky. Je %2 “(o, ky ~P2 ky) / 


ik, D -jk 
= y 7 7 y 
jez SY, e3 (lie Sy P2 ky ) (p> ky. P3 T ) Cur 2 


tjk OD 
+ (pk, +pok y,) be k ¥; tesk, Je Y) }]. (260) 


Using the definitions of Ry; and Ry3 presented earlier, (260) may be 


reduced to 
-jk. y +ik oy 
0 Y9 0 


Ro, e Y2° +e 


k -j2k, D 
2b a tite (261) 
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or, rewriting, 


“2k, D 
Bop. = Ar», Ry3 ( 2 ‘ (262) 


This is exactly the equation derived by Ziomek (1991, equation (3.9-53)) for 


this unknown constant for the classical waveguide case. 


For the unknown constant B,, Mathematica provided the following 


numerator: 


num, - G lays 42,5 43,6 44,2 45.4 ~ 1,2 42,5 43,6 4,3 a5,4 
~ 81 3 42,4 83.6 84,2 45.5 * 41,2 42,4 83.6 4,3 55 
~ 813.425 43.4 a4 2 5 6 + Ay 3 a2 4 A35 Ag 2 a5 6 
+ 412 825 03.4 943 85,6 - 21.2 22,4 235 843 a56). (263) 


Factoring (263) and collecting common terms yields the following 


generic expressionfor the numerator of B;: 


hum, = Gy lar 495 44,2 - 41,2 82.5 a43} {a3 6 a5 4 - 83.4 a5 6} 
(ays A2.4 44,2 ~ 81,2 82,4 ay 3) (a3 5 a5 6 - 435 a55)] . (264) 


Substituting appropriate expressions into (264) and using (228) yields 


the following genera/ expression for B: 


dY3,(0) dY,,(0) 


eb : é 
: ay Y2af0) 


By - 5" x(0) (¥3,(0) 





dy 
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dY3(D) dY>4(D) 


x [¥pp(v) (oo(D) ¥p4(D) —— - 05D) —— ¥300) 
las dY,(D) 
+ Youve) (p3(D) ——— YH(D) - 2D) Ypp(D) Sel / 


dY5,(¥9) dY2,(0) dY5,(0) dY,,(yo) 
dy dy dy dy 
dY,(0) dY5,(yo) dY>,(¥ 9) 


- p2(0) “ay ay Y2a0) - Yo, (0) rae) 


[p,(0) Y,(0) 


<< al 
x [Y59(v0) (o2(D) ¥3,(D) — - px(D) —— YD) 


al aca 


dY,,(0) dY3,(0) 


= lent) Y (0) (Y3,(y,) ie ay Y2al¥o) 


dY (0) 
- p9(0) —— —— (Waly Y5_(0) - ¥3_(0) Yoa(vo))] 
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dY a4 (Vo) dY3(D) — 





«lL — oat (D) Y34(D) ——- pg(D) Y3(D)} 
dY 54 (Vo) ——— d¥3(D) dY,(D) 
z (D) Yo4(D) —— - ps(D) —— ¥3(0)}]. (265) 


For constant speed of sound and ambient density, (265) becomes 


- “ik, D “ik, ee “Iky D 
B ->-(2p.k fe 3 (ae = po k le 2 2 
I¢ Yo Y, 


o 
t 


D 
ie k, 72 k, Vee ky 3 k, Je 2 


] 
+(p, ky, +19 ky, ) (p. ky +13 ky Je 2 (266) 


Using the definitions of R., and Ry3 presented earlier, and defining the 
transmission coefficient at the boundary between medium two and medium 


one T», (see Ziomek (1991, equation (3.9-58))) as follows: 


19 


(266) may be reduced to 


jak yy. -j2k. OD 
ae awe ee 3 Oe: 


By ” Ax 
-j2k D k 
[1 2 Roy Ro3 e Yo 


R 
Multiplying (268) by Ror allows us to write By in the form 


By ~ Boy Tar. 


J_ —_________________ ,. f,,, 


(267) 


(268) 


(269) 


This is exactly the equation derived by Ziomek (1991, equation (3.9-56)) for 


this unknown constant for the classical waveguide case. 


Finally, for the unknown constant A3, Mathematica provided the 


following numerator: 


num, = G, [- A, 3422 435 44 ag 4 + Ay 9 AD 3 ag 5 Ag a5 4 


~ Ay 4 823 AZ 5 a4 2 ag 4 + Ay] AD 2 AZ 5 AG 3 As 4 
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Rains 42 2 aera nines 5 34 2 42,38) 4 2451 455 
+ Ay 823 3,4 942 355 - 341 B22 A34 Ay3 5,5). (270) 


Factoring (270) and collecting common terms yields the following 


generic expressionfor the numerator of Az: 
numg ~ Gy [ay 3 a2.2 aq - Aj,2 22,3 241 + Ay 1 0,3 4,9 


- Ay 1 422 a43)! (a3 4 ag 5 - 435 as 4) : (271) 


Substituting appropriate expressions into (27!) and using (228) yields 


the following genera/erpression for A3;: 


Yen D) > a¥ee(D) 


-k ; 
A3 = ae p2(D) (y3,(D) i ay Ya0(D) 
d¥ ,(0) dY},(0) 
x L¥pal¥o) (p2(0) Se (0 iO ON re 
dY,(0) dY3,(0) 


- Yzalvo) Le2(0) —— ¥3,(0) - (0) ¥(0) — —}}] / 


dY>,(¥,) d¥>,(0) dY¥,,(0) dYp,(yQ) 


[p,(0) Y;(0) dy dy dy dy 
d¥}(0) dV 2,(¥o) dY,(¥o) 
p2(0) ores ¥;,(0) - ¥;.(0) ——) | 
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dY;(D) dY5,(D) 


x [¥pp(v0) {o2(D) ¥54(D) —— - ox(D) —— ¥300) 


dy 


dY,(D) ei 
- Ypp(¥o) Lo2(D) Yzq(D) —— - e(D) — — ¥3(0)}] 


dY5,(0) dY3,(0) 


- Lp1(0) ¥,(0) (Yaq(¥0) —— - a Vaal 
dY (0) 
- 9210) —G— (¥5a(v) ¥oq(0) - Y5q(0) Youle) 
dY 54,(Vo) dY3(D) dY5,(D) 
«LF — Coat) ¥3,(0) — - 030) 30) 
WEAe) dY3(D) <i 
ay Lent) YD) -- eD —E—YD]. (272) 


For constant speed of sound and ambient density, (272) becomes 


A an (2 P2 ky, ) tor k, *p2k, Je es 
2 | 


3e = 


eZ 


2] 
+ (p, Ky, + po ky V(p. Ky, + p23 ky Je : F (273) 


Using the definitions of Ry; and Roz presented earlier, and defining the 
transmission coefficient at the boundary between medium two and medium 


three T.; (see Ziomek (1991, equation (3.9-59))) as follows: 
Z po kK 
2 Y> 


T93 eas ae aT (274) 


Se ana aes 
Yo 3 


(273) may be reduced to 
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sik y +jk oy 
Roy e Y2°%+e 9 0 | ae 
K 
A. ° 9. 23 . (275) 
jr Dyk. <i D 
[1 7 Ro Ro3 e ‘2 Y9 ea 
(275) may be rewritten as 
-jk_ D 
ne, 
Ay. = Ar», T93 an To (276) 
-jk. D 
<3 


This is exactly the equation derived by Ziomek (1991, equation (3.9-57)) for 


this unknown constant for the classical waveguide case. 


To summarize the results for the three media waveguide with plane, 


parallel boundaries, generic expressions for the numerators of the unknown 
constants are given by (232), (247), (253), (258), (264), and (271). The 


generic expression for the common denominator is given by (226). The fact 


that these constants have a common denominator may be useful when 


constructing a modular program for numerical calculations involving these 


constants. For the genera/ case, in which speed of sound and density are 
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arbitrary functions of depth, the unknown constants are given by (235), 
(248), (254), (259), (265), and (272). Again, it should be noted that the 
denominators in these expressions are exactly the same, offering the same 
advantages in modular programming. For constant speed of sound and 
density, the unknown constants are given by (245), (251), (256), (262), 
(269), and (276). 


The fact that Afzthematica's output can be shown to be equivalent to 
results derived in a completely independent manner gives us some 
confidence in the answers arrived at by Mathematica 's LinearSolve function 
and symbolic processing. We will now program Afzthematica to solve the 
four media waveguide problem assuming all boundaries are planar and 
parallel. By obtaining constant sound speed, constant density expressions for 
the resulting eight unknown constants and making judicious selections of 
transmission and reflection coefficients, we can mathematically ‘remove’ the 
fourth medium and verify that the four media case correctly reduces to the 
three media case. This will offer us more confidence in the program's output, 


and provide some meaningful results. 


For the four media waveguide problem, the vector-matrix quantities of 


(198) are defined as follows: 


Ie 


where 


A-[ai1 a2 a3 0 0 0 0 0 
O a22 423 a24 a25 0 O 0 
0 0 0 434 435 436 437 O 
0 0 0 0 DO 446 34,7 44,8 
a5; as2 a53 0 0 0 O 0 
0 0 


464 465 466 467 0 
OnceO. Om One O cama, cman ans 
O ago ag3 ag4 a5 0 O O 


ayy ~~ pilys) YyCys) 
ay.2 - pal¥s) Yo_(¥s) 
a13 - Palys) Yoal¥s) 
2,2 = YalVo) 
42,3 = YoalVo) 
a2.4 = - Youle) 
425 - - Yop(¥o) 


43.4 - p2(¥p,) Yoo(Yp,) 
43.5 - p2l¥p ) Yon(¥p ? 
a3,6 ~~ p3l¥p,) Y5(¥p,) 
a3. ~~ p3l¥p,) Ya(¥p,) 


a4,¢ ~ 3(Yp,) Y3(¥p,) 
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(277) 
(278) 
(279) 
(280) 
(281) 
(282) 
(283) 


(284) 
(285) 
(286) 
(287) 


(288) 


4,7 = 03(Yp,) Y3(¥p,) (289) 











a4g =~ p4lYp,) Y4(¥p,) (290) 
dY (ys) 
a5, =~ dy (291) 
dY5,(ys) 
2 dy (292) 
dY,,(ys) 
a5 3 = dy (293) 
dY 25 (¥p,) 
46,4 = dy (294) 
dYonl¥p ? 
agar gy (295) 
d¥3(yp,) 
Th a areal (296) 
dY3(yp,) 
a dy (297) 
dY3(¥p,) 
dY3(yp,) 
a7 = dy (299) 
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dY4(¥g,) 





a7 nr (300) 
dY3,(Yo) 
232 ea (301) 
dY,,(¥o) 
a5 ea (302) 
dY 54 (Vo) 
ag4™- dy (303) 
dYo4(¥o) 
ag5=- dy ; (304) 
i= [B, Ara Boa Aon Bop A; B3 A," (305) 


where the superscript T indicates the transpose matrix operator (indicating 


that x is a column vector), 


b-|00000006,|' (306) 


where the superscript T indicates that b is a column vector, and 


k 
Claas 


(307) 


Since the matrix A is an eight by eight square matrix for this four 
media waveguide problem, the solution to the system can be written as 
described in (222). The Mathematica code required to perform this task is as 


follows: 
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a= (ialcl, ale2,atcs, 0, 0, 0,0, 0), 
(0, a2c2, a2c3, a2c4, a2c5, 0, 0, 0), 
(0, 0, 0, a3c4, a3c5, a3c6, a3c7, 0}, 
(0, 0, 0, 0, 0, a4c6, a4c7, a4c8), 
{aS5cl, aSc2, a5c3, 0, 0, 0, 0, 0), 
(0, 0, 0, a6c4, a6c5, a6c6, a6c7, 0), 
{0, 0, 0, 0, 0, a7c6, a7c7, a7c8}, 
(0, a8c2, a8c3, a8c4, a8c5, 0, 0, 0)): 


be3(0; 0,0, 0,0,,0, 0,G1); 


x = (Inversel[a]).b 


This code resulted in halted execution due to a singularity error in the 
evaluation of the matrix inverse. Luckily, the LinearSolve function can be 
used to solve (198) directly for this problem because A is an eight by eight 


square matrix. The following modified Mathematica code was developed: 
a= {{alcl, alc2, alc3, 0, 0,0, 0, 0), 
(0, a2c2, a2c3, a2c4, a2c5, 0, 0, 0}, 
(0,0, 0, a3c4, a3c5, a3c6, a3c7, 0), 
{0, 0, 0, 0, 0, a4c6, a4c7, a4c8}, 
{a5cl, aS5c2, a5c3, 0, 0, 0, 0, 0), 
(0, 0, 0, a6c4, a6c5, a6c6, a6c7, 0), 
(0050000. a/cG, a/c7, a7c3). 
(0, a8c2, a8c3, a8c4, a8cS, 0, 0, 0)): 


b = (0,0, 0, 0, 0, 0, 0, G1); 


LinearSolvela,b] 


This revised code ran successfully, yielding results in the same format 
as the three media waveguide problem discussed earlier [see (223). Again, 


when Mathematica functions such as Factor, Cancel, and Simplify were 
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applied to the output, the same result was returned, indicating that the 
program was satisfied that the output was as simple as it could make it. This 
fact confirms our suspicions that the symbolic algebra routines contained in 


Mathematica \ack sufficient sophistication for this application. 


We will now present the results of the program for the four media 
waveguide problem. The analysis of these results will be conducted in a 
manner similar to the analysis conducted for the three media waveguide 
problem. That is, the results will be simplified manually in order to generate 
generic expressions for the unknown constants in terms of the generic 
elements. Then, (277) through (304) and (307) will be substituted into the 
generic expressions to reveal general expressions for these unknown 
constants. Next, we will assume constant speed of sound and constant 
ambient density, and develop a set of expressions for the unknown 
constants. Finally, we will make some judicious assumptions regarding the 
reflection coefficient (at the boundary between medium three and medium 
four), the transmission coefficient (at the boundary between medium three 


and medium four), and the location of the boundary YB, (thereby 


mathematically removing the fourth medium) to show that the constant 
speed of sound, constant density expressions for the four media waveguide 
problem reduce to the results of the three media waveguide problem 
already verified. This verification will be conducted in the following order: 
Ava: Bag, Aap. Bap, Az, Bz, By, and As. We will demonstrate this entire process 
for the unknown constant A», only, and simply present the results for the 


other seven unknown constants. 
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The first output expression to be explored is the common denominator, 


denom 


denom = - aj 3 2,5 43.7 44,8 45,1 46,4 47,6 8,2 
* 11 92,5 83,9 846 95,3 86.4 47,6 46,2 * 413 42,4 43,7 848 45,1 465 47,6 48,2 
- 41,1 42,4 43,7 44.8 5,3 26,5 87,6 48.2 * 41,3 42,5 23,4 Aq a5) 46,7 47,6 ag 2 
~ 41,3 42,4 43,8 A4g 45,1 26,7 47,6 482 ~ 21,1 42,5 3,4 448 5,3 46,7 47,6 Ag 2 
+ Ay 1 42,4 43.5 44g 45 3 46.7 47.6 ag 2 + ay 3 425 a3.6 448 a5 1 a6 4 A77 Ag 2 
- Ay 1 42.5 43.6 A4.g 453 a6 4 77 Ag 2 - 41,3 42,4 43.6 agg a5) 465 277 Ag 2 
+ 411 42,4 43.6 44g a5 3 46,5 477 ag 2 - ay3 42,5 3.4 agg a5) 466 47,7 ag 2 
+ @y 3 424 435 44.8 451 46.6 47,7 ago + Ay.) 42,5 43,4 agg a5 3 G6 a77 ago 
- Ay 1 42.4 43.5 agg a5 3 46.6 47,7 Ago t+ Ay 3 42.5 43,7 a4 6 a5 | AG 4 a78 Ago 
- 41 3 42.5 43.6 447 451 864 a7 g Ag 2 - 441 42,5 43.7 a4 6 a5 3 46 4 a7 8 4g 2 
+ 4 1 425 43.6 44,7 45,3 46.4 a7 g Ag 2 - ay 3 a2 4 237 a4 6 a5 | 46,5 a7 8 ag 2 
+ 413 424 43.6 44,7 45) 465 a7.g ago + Ay 1 42.4 AZ7 44.6 a5 3 G5 a7 g ag 2 
~ Ay) 42.4 A36 a4 7 a5 3 4g5 a7 g ago + Ay 3 AD 5 Az 4 A47 a5 466 a7 g ag 2 
~ 41,3 42.4 435 847 a5 1 466 27.8 ago ~ 4,1 425 434 47 a5 3 d66 a7 9 ag 2 
+ Ay) 82.4 43,5 447 a5 3 46.6 47g Ag 2 - ay 3 825 43.4 a4.6 a5) 46,7 a7.8 ag 2 
+ Ay 3 404 435 44.6 45.) 46,7 a7 g ago + ay | 425 43.4 a4 6 a5 3 a6 7 AZ 8 Ag 2 
~ Ay) 42,4 43.5 84.6 a5 3 467 a7 g ago + Aj 2 Ar 5 43.7 448 45) 46,4 47,6 ag 3 
~ Ay 1 42,5 43.7 agg 45 2 a6 4 a7 6 ag 3 - Ay 2 a2 4 A37 agg a5 A65 AG ag 3 
+ Ay 1 82.4 43,7 a4.g a5 2 46,5 47.6 Ag 3 - A1.2 A255 43.4 a4. a5) 267 a7 6 ag 3 
+ Ay 2 42,4 43,5 A4g 45) 267 47.6 ag 3 + Ay) A255 A324 agg a5 467 a7 6 ag 3 
~ 8y,1 82,4 435 24g 45 2 AG 7 a7 6 Ag 3 - Aj 2 A205 43.6 agg a5) 46,4 A7,7 ag 3 
+ Ay 1 425 43.6 44.8 45 2 AG 4 a7 7 ag 3 + Ay 2 Ad 4 AZ 6 agg a5) AGS A7,7 ag 3 
~ Ay 1 424 43,6 A4.g a5 2 465 a7 7 ag 3 t Ay 2 a2 5 A324 agg a5 AG 6 a7 7 Ag 3 
41.2 42,4 43,5 44,8 45.1 46,6 47,7 Ag. - ay) B25 43.4 448 a5 2 466 a7,7 Ag 3 
+ 411 42,4 435 a4.g a5 2 46.6 A7,7 Ag .3 - Ay 2 25 43,7 a4 6 451 AG 4 ADs as 2 
* 1.2 42:5 43.6 44,7 45,1 46.4 a7, Ag3 + 411 425 43.7 44 6 45,2 264 a7 g ag 3 
~ 84,1 42,5 43.6 44,7 45,2 26,4 a7,.g Ag 3 + Ay 2 A24 a3.7 a4 6 a5 46.5 a7 8 ag 3 
~ 41 2 40,4 43.6 47 a5 1 46,5 47g agg - ayy Ar 4 a37 a4 6 a5 2 465 47g ag 3 
+ 11 42,4 43.6 44,7 a5 2 465 a7. Ag 3 - ay 2 25 34 47 aS 466 AZ g ag 2 
+ 41,2 49,4 43.5 A247 a5. 46.6 A7.g Ag. + 411 42,5 43.4 44.7 a5 2 46.6 47.8 ag 3 
~ 41 1 42.4 435 447 a5 2 46,6 a7 g ag 3 t+ Ay 2 A25 a3 4 a4 6 a5 | 467 a7 g Ag 3 
~ 81,2 42,4 43.5 a4 6 45,1 467 A798 483 - 411 425 43.4 4 6 a5 2 467 a7 8 ag 3 
+ Ay 1 42.4 43,5 44.6 45 2 467 A7.g ag 3 - ay 3 a2 9 a37 a4g a5) 465 47 6 Ag 4 
+ 12 42.3 A3.7 aq .g a5 465 a7 6 Ag .4 - Ay) 423 437 a4 g a5 2 465 a7 6 ag 4 


+ A, | 42,2 437 Aq.g a5 3 Ag5 a7 6 ag 4 + Ay 3 AD 2 35 agg a5 1 A67 A76 ag 4 
412 22.3 43.5 a4.g a5) 46,7 47.6 4g4 * Ay) 223 43,5 agg AS 2 AG 7 a7 6 Ag 4 
A141 42,2 23,5 agg as 3 A67 A7.6 ag 4 * Ay 3 A2 2 436 a4g 51 AG 5 a7 7 ag 4 
~ 41,2 42,3 43.6 24g a5 465 a7 7 agg + Ay) 42,3 43.6 44.8 45,2 AG 5 a7 7 Ag 4 
411 42,2 43.6 a4.g as 3 46,5 477 4g 4 - 21 3 422 A35 agg As 1 466 A777 Ag 4 
+ 12 42,3 435 44g a5 1 466 477 ag 4 - ay 1 42,3 43.5 a4 g as 2 Ag 6 a7 7 ag 4 
+ Aj 1 42,2 43.5 agg as 3 AG 6 a7 7 agg + Ay 3 a2 2 a37 a4 6 a5) AG 5 a7 g ag 4 
A) 2 42,3 43.7 A464 a5 1 46.5 A7.g ag .4 - 41 3 a2 2 43.6 A477 a5 | A655 arg ag 4 
* 12 42.3 43.6 a47 a5 1 265 A78 ag 4 + Ay 1 A 3 A37 a4 6 a5 2 Ags a7 g ag 4 
— Ay 82,3 3.6 44,7 45 2 46.5 47,8 agg - 21,1 42,2 a3.7 Ag 6 a5 3 AG 5 a7 8 ag 4 
+ 11 42,2 43.6 A4.7 a5 3 46.5 a7 9 ag 4 + A 3 G22 35 a47 a5 | A664 a7 ag 4 
~ 41,2 42,3 83,5 ag 7 as 46.6 47,8 ag 4 + Ay 1 42,3 43,5 a4 7 a5 2 AG 6 a7 8 ag 4 
411 42,2 43.5 24.7 a5 3 46.6 a7 g ag.4 - Ay 3 A222 Az 5 a4 6 AS | AG 7 ATs ag 4 
+ Ay 2 42,3 43.5 aq 6 a5 1 46,7 a7 8 agg - Ay 1 A23 a3 5 4G a5 2 a6 7 a7 8 ag 4 
+ A, 1 82,2 435 A464 45,3 AG 7 7g agg + Ay 3 A2 2 AZ 7 Ag g a5) 464 a7 6 ag 5 
~ 41,2 42,3 43,7 4g 451 46,4 47,6 Ags + 41,1 42,3 43.7 4g a5 2 a6 4 A764 ag 5 
~ 8,1 42,2 23,7 24.8 45 3 464 47.6 485 ~ a1 3 22 434 a4 g a5 a6 7 a7 6G ag 5 
+ 1.2 42,3 43.4 A4.g a5 1 46,7 a7.6 Ags - 411 42,3 3.4 A4g AS 2 d67 A765 ag 5 
+ Ay 1 82,2 A3,.4 A4g a53 467 27.6 485 - 1,3 Ad 2 A364 agg 5 1 864 a7 7 Ags 
+ Ay2 42,3 43.6 44.8 45.1 46,4 A777 ag5 - ay 1 A2 3 a3 6 agg a5 2 Ag 4 a7 7 Ags 
41,1 42,2 43,6 a4.g a5 3 46.4 A77 gs * Ay 3 G22 AZ .4 agg a5) 466 47,7 a5 
~ 1,2 42,3 43.4 44.8 45,1 466 47,7 Ags + Ay A23 43,4 agg a5 2 466 77 Ag 5 
Aj 1 42,2 23.4 Aq .g a5 3 466 77 Ags - A, 3 A 2 AZ7 a4 6 a5 1 464 a7 8 Ags 
+ Ay 2 42,3 43,7 a4.6 451 a64 a7 8g Ags + A132 Ar 2 a3 6 A475) 464 a7 8 Ags 
1,2 42,3 43.6 44,7 45,1 26,4 a7, Ags - Ay 1 42,3 43.7 a4 6 A502 46.4 a7.8 ag5 
+ 1.1 42,3 43,6 24.7 25,2 46.4 a7 g Ags + Ay 1 A222 a7 A446 a5 3 AG 4 a7 g Ags 
~ 411 42,2 43.6 a47 As 3 46.4 7g Ags - ay 3 22 AZ 4 Ag? aS 466 a7 8 4g 5 
+ 1.2 42.3 43,4 a4.7 a5 1 46.6 a7.8 ag5 ~ ay 1 42,3 a3 4 a47 aS 2 AG 6 a7 8g ass 
41,1 42,2 43,4 A477 45 3 46.6 a7 g 4g5 + 413 A22 AZ 4 A466 a5 46,7 a7 8 ag 5 
412 42,3 234 a4.6 45,1 267 A7.g ag5 + Ay 1 2,3 3 4 a4 6 45,2 26,7 a7 9 Ags 
- Aj,1 82,2 43.4 24.6 453 26,7 A798 Ags. (308) 


+ 


+ 


It is readily apparent that (308) is a rather complicated expression 
involving 128 terms. In order to simplify the algebra somewhat, we have 


divided the expression into four distinct subexpressions, such that 
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denom = denom, + denomp+ denomc + denomy. (309) 


where denom,g is the sum of terms involving the element ag 9, denomg is the 
sum of terms involving the element ag3, denomc is the sum of terms 
involving the element ag 4, and denomp is the sum of terms involving the 
element ag 5. This simplification allows us to work with 32 terms at a time, 


and yields the following generic expressions for the various parts: 


denoma = ag 9 ee 5.1 - 444 as.3} lap.5 lt a4.5 7.6 - 84,6 47,5) 
x (a5. 46,7 - 43,7 a6 4) + (a4. a7. - 44.8 ay) A34 866 ~ 43.6 46,4 
— 42.4 lass a76 ~ 44,6 a75) {3,5 46,7 — 43,7 465 
+ (a47 47,8 - 44, 27,7) 1435 46.6 - 43,6 365 | (310) 


denomp = ag 3 fa; 45,2 - ay.2 as.) laa,s [a4 a7.6 - 44,6 47,5) 
‘s (a3, 46,7 — 43,7 46,4) + (a4 a7. - 448 a77) 43.4 46.6 ~ 43.6 46,4 
- 424 (ass a7 6 ~ 44,6 ays) 435 467 - 43.7 465 
+ (ay7 arg - a48 a77) (ars 466-436 46s5}), (311) 


denomce = ag 4 lags 47,6 - 44.6 27,8) (a3.5 46,7 - 437 a6} 
- (a49 a7. - 44.8 a7 7) 235 46.6 ~ 43.6 465 
x (eae (a, a5 - ay as 3) + a2 3 (ayy a52- 412 as.) (312) 


and 
denomp = ag 5 lass 47,6 - 44.6 47,5) (a5.4 46,7 - 43.7 aca) 


* (a, a7 .g - 448 a7 7) A234 466 - 43,6 46,4 
x tan (ayy a5 3-7 a1 3 as 1) + a2 3 (ay 9 a5, 7 ayy a5 2) : (313) 
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Substituting (310) through (313) into (309) and combining like terms 


yields the following generic eypressionfor the common denominator : 


denom - (ass 47.6 - 44,6 47,8) { 23,4 6,7 - 43,7 a6.) 

; (a4 a7,.8 - 448 a7 7) 43.4 46,6 ~ 43.6 46,4 

x (ao 5 ag 2 - a2 2 ag 5) (a, 3 as 1 - ary as 3) 

: (a2; 48.5 ~ 42,5 ag.3) (a, 2 45.1 - ay as,2) 

; L(a4.s a7,.6 ~ 44,6 a7,g) 43,5 46,7 - 43,7 465 

+ (a4.7 87,9 - 44.8 47,7) 43.5 46,6 - 43.6 465 

x (ao 4 43,2 ~ 42.2 ag.4) (ay; 45.4 - 44,1 as,3) 

+ (ap. 45.4 - 42,4 43,3) (ay2 25,4 - 24,1 5,2)1). (314) 


Substituting (277) through (304) into (314), using the facts that ys = 0, 


YB, = D;, and YB. = Do, and performing some algebraic manipulations results 


in the following genera/erpression for the common denominator : 


dY, (D2) dY3(D2) 


denom = [ (04(0,) Y3(D,) a - p4(Dy) dee 


5 Yi(D,)) 


dY54(D,) dY,(D,) 


(AWS) ahh) a 


dy 


dY,(D2) dY;(Dz) 


+ p4(D2) dy Y,(D,) = 3(D2) Y,(D) eS 
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dY3(D,) 
dY>,(D,) 


V0) = oD Ya, Dy) aa 
3 


dy 
x p3(D,) 


dy 


de yell 1) 


ay Yanl¥o) 
x [v3.09 qe sce 


dY,,(0) dY ,(0) 
2a 


 92(0) “Go ¥oq(0) 
x {p,(0) ¥}(0) 


dY>,(¥o) 
dina.) 


Oo) ei Ne tye) 


dy 
dy 


dY,,(0) dY ,(0) 
2a 


: dY3(Dz) 


4(D ) 
[( ) Y3(D ) = p4(D2) 2 
( IEA onl 4 
= 3 D v 


dy 
dy 


dY4(D,) 
dY>,(Dy) 7 
x 03(D,) dy 3 


dY ,(D,) 
dY,(D2) 


4 


dy 
dy 
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dY5,(Dy) dY3(D,) 


x p3(D,) Y,(Dy) - p2(D,) Yop(D1) i. 


dy 


GY on (¥o) ley y.) 


x (Y3.(y,) dy Y3u(Vo)) 


ae 
dY,,(0) dY ,(0) 
x (o1(0) ¥4(0) ——- 20) —— ¥p,(03) 
dY5,(¥o) dY 54 (Yo) 
+ dy Yop(Vo) = LeAUA) iy * 
dY,,(0) dY ,(0) 
« (6,(0) 0) — - en —¥3,(0)) |. (315) 


If we now assume constant speed of sound and constant density in a 
specific medium, the depth-dependent functions in (315) become complex 
exponentials and the denominator becomes (using the subscript c to indicate 


the constant speed of sound assumption) 
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-jk D 
denom, - 2 ky. e | Yas [(, Ky, + p2 a ) (5, ky, =p ky.) 


3 (5, ky, - p2 ky ) (0; Ky, + p2 k, 


~jk D, ~iky D, *iky 


7 Pi 
x (po, k - p3 k Je (SE aie ak igh se 
"3 %4 


Mork, -erk, ) esky -erk,, | 


x (oak, resk, Je 2° ¢ 
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a De aD See AD 
Ye 9% 27a | . 


(316) 


Now that the denominator has been simplified, we'll concentrate on 


obtaining expressions for each of the unknown constants in the order stated 


above. The first constant is 


where 


Ada 


num» 
denom ’ 


NuMy = Gy [- 413 42,5 43.7 44.9 45,1 46,4 27,6 


411 42,5 43,7 44.3 45.3 a6 4 a7 6 + 


~ 411 42,4 23.7 44g 45 3 46 5 a7 6 * 


ay 


+ 


~ Ay 1 82,4 a3,5 a4 6 a5 3 467 a7 gs]. 


41,3 42,4 43.5 44g a5) 46,7 a7 6 - 
4141 42,4 43,5 agg a5 3 G7 a7.6 +t 
41,1 42,5 43.6 44.9 45,3 6 4 a77 - 
411 42,4 43.6 44g a5 3 465 877 - 
41,3 42,4 43,5 44g a5 1 46,6 77 + 
411 42,4 43,5 a4qg as 3 66 a7 7+ 
1,3 42,5 43.6 447 451 A6.4 a7 - 
41,1 42,5 43.6 44,7 45,3 46.4 47,8 - 
413.424 43.6 A4.7 451 AG 5 A798 + 
411 42,4 43,6 44,7 A523 G5 a7 + 
413 42,4 43.5 84.7 a5 1 46.6 47,8 - 
411 42,4 43,5 a47 a5 3 Ag 6 a7 8 - 
@1 3 42,4 43,5 a4.6 45 1 Ag 7 a7 9 + 


413 42.4 437 agg a5) 465 47.6 
413 42,5 43.4 a4 g a5 Ag. 7 a7 6 
41,1 42,5 43.4 44g as 3 ag 7 a7 6 
41,3 42,5 43.6 44.8 45 1 46.4 477 
413 42,4 23,6 44,8 a5 1 46,5 a77 
41,3 425 43.4 44g a5 1 46.6 A772 
411 42,5 43.4 a4.g a5 3 46.6 AD 7 
41,3 42,5 43.7 44.6 a5.) Ag 4 a7 8 
41,1 42,5 43,7 24,6 45,3 46,4 478 
41,3 42,4 43,7 44.6 a5) 46,5 a78 
411 42.4 43,7 84.6 a5 3 465 a7 
41,3 42,5 43,4 44,7 a5.) 466 47,8 
41,1 42,5 43,4 447 as 3 46.6 a7 8 
81,3 42,5 43,4 44.6 45,1 46.7 47.8 
A1,1 42,5 43.4 44.6 45 3 467 a7 8 


(317) 


(318) 


Factoring (318) and collecting common terms yields the following 


generic expression for the numerator of Ada: 


NUM, = G, (a3 asi - 411 as,3) 


lay 5 [lags 47.6 - 44.6 ays) 


x (as 4 46,7 - 43,7 ag 4J * (ay9 a7.3 - 44.8 a7 7) 434 466 - 43.6 464 
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- 424 lass 476 - 446 a7 s) leas 467 ~ 43,7 46 Fi 
: (a47 a78- 448 a7 7) 435 46.6 - 43.6 465 | (31) 


Substituting appropriate expressions into (319) and using (307) yields 


the following genera/ expression for the numerator of Ada: 


i —— dY»,(0) dy,(0) 
num, = > (o4(0) ¥,(0) Tae PAD ae Ya) 
dY,(Dq) dY3 (D9) 
x [. Yopl¥o) I (o3(D») Y3(D,) dy 4 p4(D>) dy Y,(D2) 
dY¥2y(D1) dY,(D,) 
x 3(D,) dy Y3(Dd) : p2(D,) Yop(Dy) ap 
dY3(D2) dY (D2) 
= (04(D>) ay Ya(De) = p3(D) Y,(D9) dy 
dep dY3(D,) 
x Lp3(D,) ay YatDy) - p2(D1) Yo, (Dy) ora 
dY (D3) dY3 (D3) 
+ Y¥zu(vo) [(o5(D2) ¥3(D2) —G— - 04(D2) —— ¥4(02)) 


{39 


dY 9,(Dj) dY(D,) 


¥ LpxD,) Y,(Dy pay anon 


ma a 


dY, (D2) dY,(D2) 
+ (p4(D2) ay Y4(D2) - p3(D2) Y3(Dz) ie 


dY >, (Dj) dY, (Dy) 


x Co5(D)) —— ¥3(D)) - 6211) Yoq(04) =i _ (320) 


Once again, making the constant speed of sound and constant density 
assumptions and substituting the appropriate depth-dependent expressions 
allows us to write the numerator of Ao, as (using the subscript c to indicate 


the constant speed of sound assumption) 
num, - k -pok Je 4 
200 2n (6, Yy P2 y, ) 


-jk D, +jk y -jk D, +jk, D 
lla +pok ) (3 k apa K Je Yo $e %2%e %3 26 Came 
Yo Y2 V4 Y2 


alt +k, siky Dd, Mei 


losk, -erky )(esky +oak, )e cp 2) eS te 
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+{k D, -jk 
oe ity Yo 


+t D, -jk D 
2 ] 
MM ho mk lem te aoe te oll: 


(321) 


Thus, the genera/ resu/t for Aza is formed by dividing (320) by (315) 


and is given by 


me —— d¥,,(0) dY ,(0) 
Ara oe (,(0) Y,(0) Lie: p2(0) neo) 
; dY, (D2) dY, (D9) 
x [. Ge) I (o4(D,) Y3(D) arte p4(D.) rae Y,(D2) 
dYoy(D1) dY (Dj) 
x 03(D,) ay ¥3(D1) e p2(D;) Oy) ae 
dY3(D2) dY;(Do) 
a) Gy Yala) - 93(Da) ¥3(D2) —o— 
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+ 


dY 4 (D,) dY,(D,) 
x (edn a 30D 3 e2(D,) Vay) ae 


ns = 
+ Y5p(vo) [(o5(D2) ¥5(Dz) —— - pa(Da) —y— Y4(D2)) 


dY5,(D,) dY,(D,) 
x (p,(D)) By VHD) ~ p2(Dy) Yog(D1) = 


= = 
+ ((04(D2) —— Yi{Da) - (Dp) ¥3(D2) —— 


dY54(D1) dY3(D,) 
« (63(01) —— Y4(D)) - p2lD) YoxlD) II 


dY, (Do) sat 
[|(,,(0,) (D5 eer dy i pa (D sae Yi(D,)) 


dY5,(D1) dY,(D,) 
x p3(D,) ay 3D) z p2(D,) aD) Tih 


— ain 
+ (p4(D2) —— ¥4lDa) - 0(D2) ¥5(D2) ——) 
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dY54(Dy) dY (Dy) 
x | p3(Dy) ay YD) - po(Dy) Yop(Dy) ar 


NEG CVO 


x [(Y3,¢v0) << ne) 


dY5,(0) 


x {o1(0) ¥(0) —F— - p20) Vp, (00) 


dY5,(¥o) dY54(¥o) 


t; dy Yop(Vo) a Yoa{¥o) dy 


dY,,(0) dY ,(0) 


x (p4(0) ¥;(0) = - p(0) 


: v7,(0)}| 


dy 


dY (D2) dY (Dy) 
: I (p5(,) Y3(Da) ae p4(D>) ae Yi(D,)) 


& 


dY (Dj) dY,(Dj) 


Fy VHD) - 92101) ¥34(D)) ——} 


dY3 (D2) dY (Da) 


: p4(D) dy Y,(Dq) is 03(D) Y,(D9) 


y dy 


143 


dY5,(D)) dY3(D,) 


x e3(D,) ay 30 - p2(D,) Yop(Dy) | 

dY54(Vo) dY5, (Yo) 

x [(Y3.(y9) co % dy Yin(¥o)) 
dY,,(0) dY ,(0) 

x {o1(0) ¥}(0) ——- 20) —— ¥,(00) 
dY59(Vo) WY ap (¥o) 
| a Yop(Vo) - Yoal¥o) dy 
dY5,(0) dY,(0) 
x {p4(0) (0) —— - p20) —— ¥5q(0) I. (322) 


For constant speed of sound and ambient density, A, reduces to 
a 


LS ity» 
Ara” On lok, -p2k, )e 4 


2 


-jk, D, +ik, y. -ik,, 0, <iieeD 
1 0 y, 2 see || 
x (e ky, +p ky.) (3 Ky - P4 ky.) ee e "2°e "3 “emma 
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-jk. D j 
yD1 By, Yo. cit 
Ye % he ee 


Sak +p2 kK )( 

: +jk D, -jk 

y, Pa ky le ie 1 ay vik, D, -ik, D 
ea ee Bey 
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-jk D, -jk. D, +jk. D 
Je Memes, Sears! 


-jk D,  +jk. D 
Je ee 
V4 


-jk OD 
«(oak +a 3% i]. (323) 


-jk_ D 
Eliminating the common term e "4 : dividing numerator and 


denominator by bs KY + po ky ), and using the definition of Ro, given by 
2 I 


(238) yields 


«Ulosk, soak, (oak, nek, De 2 € 2 € "3S Ae 
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-jk D +jk D, -jk D 
1 ¥, 76 v3] 


- Roy losk, -e2ky ) esky + eak, Meee 
4 


(324) 
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Dividing numerator and denominator of (324) by a ky + Po ky ), and 
2 4 
using the definition of Ro3 given by (241) yields 


-jk, “iky D, ve Le “ik, D, +iky D, 
= Se = 2 2 3 3 
2a, Ro \e Ky P4 ky) e e e € 
4nk 
y 


A 
2 


a, D, +k, % MS D 


aR peak +p4k Je 2° ee ode Cae a) IC™ 25 
%4 Y3 


+k D, -jk, y. -jk D, +jk. D 
"Ras losky -eak, De Y9 le 2) ee aT ee 


rik. D, -ik y +i. Dy -jk 1 
+p,k Je yg Yo Pe 3 eg 
Page 


*iky D, “ik, D, tik, D, 


[a3 ny ky, - P3 k, ) e 72 Co ee 3 
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-jk D, <+jk D, -jk. D 
- Ro; Ro3 (on +o3k Je ee ae! 
Y2 V4 


Dividing numerator and denominator of (325) by ie ky + p3 ky } 
3 4 
defining a reflection coefficient at the boundary between medium three and 


medium four, R34, as 


Co SS (326) 
ia kK + p3 kK ) 
4 Y, 3 Ky, 


and multiplying through by the -1 appearing in the numerator yields the 
following desired expression: 
+jk Pied i= ¥ dee-ike AD, 2D.) 

Ady. aR [Ree yr lane - ed 


Yo 
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=jk) (Deny eee (Dee) 
+Ro3e — — e ¥3 2 


sik. (D,- y_) -jk. (D,- D,) «jk. (D,- y ) «jk. (D,- D,) 
Ro Regie y, | “o ees 2 ae Sine aligyae Ee 2 PI / 
-je (D,- D,) +t D +k (D,- D,) +ik. D 
[RARE ue : e Lo) lve 23 ‘ c ye 
-jk, (D,- D,) -jk, D sik (D,- D,) -jk D 
~RopRyg€ 22 °C UMUReER Res cei ']. 327) 


Now that we have derived the equivalent classical expression for the 


four media waveguide problem for the unknown constant Aga we must 


verify that it reduces to the proper expression if the fourth medium is 
removed mathematically. In order to conduct this evaluation, we must 


assume the following: 


D, = D, = D (that is, Logline we) (328) 
and 


R34 <0. (329) 


Substituting conditions (328) and (329) into (327) reveals 
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Equation (330) is equal to equation (242). Since (242) has already been 


verified, we can conclude that our derivation of the solutions for oe and 





(ae 


+jk. (D-y_) 
+e Yo A 


| / 


-jk. D 
- Ray Roz e ee i 


hence, Ang. for the four media waveguide problem ts correct. 


Continuing with the analysis of the four media waveguide problem code 


results, for the unknown constant B,, Mathematica provided the following 


numerator: 


oe 


+ 


num; = G, [a),2 425 43.7 44.5 a5 1 46.4 476 


Aj 1 42.5 43,7 44g a5 2 a6 4 a7 6 - 
41,1 42.4 43,7 44.8 45.2 46.5 47.6 - 
412 42,4 435 44g 451 46.7 476 * 
Aj 1 42,4 435 44g a5 2 467 a7 6 - 
41,1 42,5 43.6 44.8 45,2 46.4 47.7 + 
41 | 42,4 43,6 44g 45 2 a65§ a7 7 +t 
412 42,4 43,5 44.8 45.1 46.6 47,7 - 
A, 1 42.4 43.5 44g as 2 46.6 a7 7 - 
41,2 42,5 43,6 44,7 451 464 a7.8 + 
411 42,5 43,6 44,7 45 2 agg a7 g + 
412 42.4 436 44.7 a5) 46.5 a7 g - 
41.1 42,4 43,6 44,7 45.2 465 a7, - 
412 42,4 43.5 44.7 451 46,6 47,8 * 
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@12 42.4 43.7 44g a5 1 465 47.6 
A} .2 42,5 43,4 44g 451 46,7 47.6 
A11 42.5 43.4 44.8 45 2 467 a7 6 
@1.2 425 43.6 44.8 45.1 46.4 477 
Ai 2 42,4 43,6 44.8 451 46,5 47.7 
@1 2 42.5 434 44g a5) 266 477 
411 425 43,4 44.8 45 2 46.6 47,7 
Aj 2 42.5 43.7 446 45 1 464 a7 8 
811 425 43.7 44.6 45 9 464 78 
A 2 42,4 43.7 44.6 451 465 47.8 
411 42,4 43,7 44.6 45 2 465 47.8 
412 42,5 43.4 44,7 45,1 46,6 47.8 
41.1 42,5 43.4 44,7 45 2 26.6 47.8 


~ Ay 1 42 4 435 ag 7 as 2 466 a7 gt ay 2 aD § AZ 4 AG 6 a5 AG 7 a7 8 
~ 81.2 42 4 a3 5 ag 6 ag AG 7 AD g - ayy AD 5 434 44.6 AS 2 467 a7 


+ 411 42,4 435 84.6 45 2 467 a7 9] . (331) 


Factoring (331) and collecting common terms yields the following 


generic expression for the numerator of Bo,: 


num; ~ G, (ay 4 a5 9 - ay. a5 1) lays ane a7 6 - 44.6 47.8) 
“ (5.4 46,7 ~ 43,7 46,4) * (a, a7,.g ~ 44.8 a7 7) 43.4 46,6 ~ 43.6 464 
~ 42.4 (a4 47.6 ~ 44.6 a75) 43,5 46,7 ~ 43,7 465. 
+ (a47 47,8 - 44,8 47.7) La35 46,6 - 43,6 465] (332) 


Substituting appropriate expressions into (332) and using (315) yields 


the following genera/ expression for Boa: 


f dy} (0) dY;,(0) 
Bra ~ 5 (p2(0) ——Ypa(0) - (0) ¥,(0) — —) 
dY, (D>) dY3(D,) 
x le Yop(¥o) [(o3(D,) Y3(D,) dy = p4(D2) dy | Yi(D,)) 
dY2p(D)) AY (Dy) 
x {p3(D,) dy Y,(D;) = p2(D,) Yop(Dy) ee 
dY,(D2) _ d¥, (D2) 
+ (04(D2) —S— ¥4l2) - es(D2) ¥5(D2) ——) 
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a ; Die 
) = p(D,) Yop 
Y3(D 
dy 
x p3(D,) 


) dY3 (D9) Hen 
dY ,(D2 mi 


dy 
dy 


dY3(D,) 
ap(D1) 
Yop 


rae (D,) Sea 
)~ p(D1) Yop 
Y3(D,) 


dy 
) 
x p3(Dy 


dy 


4 ( 


2 dy 
p4(D>) a = 


dy 
) 
x Lo3(Dy 


dy 


Y,(D) 
(D2) ie? hale 
- P4 


dy 
) 


dy 


i 


er } 
peta) aie 
) o p(D;) Yop 
Y,(D, 
dy 
x p3(D,) 


dY, (D2), 
(Do) 
dY( 


2 

p Y,(D2) - p3(D2) ¥ Fr 
4 o 

+ 4 


dy 
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aye Dn) dY,(D,) 
ier Sy V(r) - 02(D1) Y>9(D1) —— 


Yon (V) amas (70) 


x Y>a(Vo) dy Sa Teli) 
a dY ,(0) 
x [p4(0) ¥;(0) ——— - p,(0) —— Y;,(0) 
dYoa(¥o) ——— d¥an(¥o) 
+ dv Yon(¥o) = Yoql¥o) dy 
_— d¥,,(0) dY (0) 
x (,(0) WING p2(0) a 2.1005] 
dY ,(D) dY (Dy) 
- [(o5(D2) ¥3(D2) GE - pal) —— Y40D2) 


GY 95(D)) ; dY3(D,) 
x {p3(D,) ay Y3Dy) - p2(Dy) Yop(Dy) —— 


= ha 
zs p.(D5) so Y,(D9) - p3(D2) Y,(D2) ——— 
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dY4(D,) dY3(D,) 
x Lp3(D,) wa 0 - p2(D,) Yo,(D1) dy 


AYotly,) dYosy,) 


x [(y3,y¢) ee a 


Yoneye) 
dY,,(0) dY ,(0) 


x (p,(0) ¥;(0) sear (0) aarti) 


dN sty) dY 54 (Vo) 


+ ay Yano) - Y,,(¥o) 


dy 
dY5,(0) dY ,(0) 


eh ea ile (333) 


x (9 4(0) ¥;(0) re 


dy 


For constant speed of sound and ambient density, (333) becomes 


-jk D, +ik y  -jk. D, +jk D 
eee -p,k Je Yo le HD @ oe ¥, | 
Y4 Y, 
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2 


“ec 0, ciety th >. See 
-(osk +p2k, )(osk, soak, Je %'e %2%e %e % ‘| / 
2 3 4 3 


D, “ik, D, rik, D, 


rik, 
2 es e 3 


x (5, Ky, - p3 k e 


+ (p, Ky, + pe Ky ) (ee Ky, + P2 ky) 
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: (5, Ky, =e ay ) (5, Ky, - P92 k, 


-jk D, «jk. DD, -jk 41] 
Yo te 3 2 e 3 IY (334) 


Using the definitions of Ra,, R23, and R34 presented earlier, (334) may 


be reduced to 


4nk 

V9 
je. (D,- y.) «jk. (D,- D,) 
+Roge 72 eae ae 


ay 


“ik Dk (D> DY) Bite D. Se D,- Do) 


e158) 


Using conditions (328) and (329), (335) may be reduced to 


pee [ee Pr 5! 


ze (D,- y ) 
2a¢ +Roge 72 : ]/ 


+jk D -jk D 
[. ap, - Roy Roz € a2 |. (336) 


-jk, y 
Factoring e *2 ° out of the numerator of (336), and dividing the 


+jk. D 
numerator and denominator of (336) bye *2 reveals 


158 


Equation (337) is equal to equation (250). Since (250) has already been 


verified, we can conclude that our derivation of the solutions for B,. and 


-j2k. D 
| 


[1 - Ra Rage”? 


hence, for Boa. for the four media waveguide problem is correct. 


For the unknown constant A, Mathematica provided the following 


numerator: 


oe 


Zt 


numy - G; E @1 3422 837 aqg as) 265 476+ 


41,2 42,3 43.7 44g a5 1 ag5 476 - 
41.1 42,2 43.7 a4.g a5 3 465 a7 Gg + 
A} 2 42,3 43.5 aq.g 45 1 467 a7 6+ 
41.1 42,2 43.5 aq.g a5 3 a6 7 arg + 
@1 2 42,3 43.6 44g a5) 465 877 + 
411 42,2 43,6 44,8 45,3 46.5 a7,7 ~ 
412 42,3 435 aq.g a5 1 466 77 - 
41,1 42.2 43.5 aq.g ay 3 ag 6 a7 7 +t 
1.2 42,3 43.7 44.6 a5 1 465 a7 § - 
41,2 42,3 43.6 44,7 45,1 465 4718 + 
41.1 42,3 43.6 44,7 a5 2 46.5 a7 ¢ - 
41.1 42,2 43.6 44,7 453 a6 5 a7 g + 
Aj ,2 42,3 43,5 44.7 a5) 46.6 a7.9 + 
41.1 42,2 43.5 44,7 45,3 46.6 a7 8 - 


aye 


Ai 1 42,3 43,7 44.8 a5 2 46.5 47.6 
413 42,2 43.5 44.8 45.) 46,7 a7 6 
411 42,3 43.5 a4.8 a5 2 46,7 a7 6 
413 42,2 43.6 44,8 45) 465 877 
A111 42,3 43.6 44.6 a5 2 465 a7 7 
41,3 42,2 43.5 a4.g a5 1 466 477 
411 42,3 43,5 44,8 45,2 46,6 47,7 
41,3 42.2 43,7 a4 6 a5) 465 a7 8 
413 42,2 43.6 44,7 451 465 a7 8 
411 42,3 43.7 44.6 45 2 465 87.8 
@1,1 42,2 43.7 44.6 a5 3 465 a7 8 
41,3 42,2 43.5 44.7 451 46 6 87 8 
1.1 42,3 235 44,7 a5 2 46 6 a7 8 
413 42,2 43,5 44.6 45,1 46,7 A7,§ 


+ A, 2 42.3 a3 5 a4 6 a5) 267 a7 8 - Ay 1 223 A355 a4 6 a5 2 AG 7 AZZ 
* Gy 1 422 435 aq.6 45 3 467 a7 gs) . (338) 


Factoring (338) and collecting common terms yields the following 


generic expression for the numerator of Ady: 


num, = G; [Cass 47,6 - 44,6 7,8) (3,5 86,7 - 43,7 ass) 
iJ (as7 a7.8 ~ 44.8 a7 7) 435 26.6 23.6 465 
x lay» (ay 3 45.1 - ay.) 25.3) + a23 (ay 1 a52- a42454)] . (339) 


Substituting appropriate expressions into (339) and using (315) yields 


the following genera/erpression for Asp: 


rr dY (D2) dY3(Do) 
Arn = ie [[(.,00,) Y3(D) dy ze p4(D9) dy ¥i(D2)) 
dY >4(D,) dY,(D,) 


x {ox(0) —F— YD) - xD) Yq) I} 


dY (Do) dY (Dy) 


- p4(Do) Y4(D2) p3(D>) Y,(D2) dy 


dy 


dY>,(D4) dY,(D}) 
x (o3(D1) go Y3(Dy) - (D1) ¥>y(1 ——} 
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dY5,(0) dY (0) 


x LY 34(¥0) (6110) ¥4(0) —— - 02(0) —— Yoa(0) 


dY ,(0) dY,,(0) 


+ Yo,(Vo) p2(0) Gy Yaa(0) - 91(0) Y(0) 





MN 


dy 


dY,(D,) dY3(D,) 


[[(o,00,) Y3(D9) dy p4(D>) Y (D2) 


dy 


dY5,(Dj) dY,(D,) 


x e3(D)) a 


dY (Do) _— dY¥4(Da) 
+ (p4(Dy) ae Vpheps(D>) (Ds) =e 


dY 54 (Dj) d¥3(D,) 
dy Fa(Di) - po(D,) Yo,(D,) dy 


x p3(D,) 
dYon (Vo) dYo,(¥o) 


mee) ge = = age Neal Va) 


dY,,(0) dY,(0) 


x {5,(0) Y (0) dy p2(0) “ay Y2a(0) 
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AY nao) dY 54 (Vo) 


argo Yon(¥o) - YoalVo) 


dy 
on dY,(0) 
x {p,(0) ¥,(0) —— - p9(0) —— —¥2,(0)] 


dY (Dy) 
- [(ps(D2) ¥3(D)) —— reel (0))) pee 7 Y4(D2)) 


dY 5,(Dy) dY,(D,) 


x p3(D;) Y,(D,) = e2(D,) Yop(Dy) dy 


dy 


dY4(D2) — 
+ e4(D2) —{—— dy Y,(D9) - p3(D2) ¥,(D>) sae 
dY54(Dy) dY,(D,) 


x {o3(D,) ay Y3(DD) - p,(D,) Yop(D1) ea 


dY5n (Vo) dY5,(¥o) 


| Katy) maa a aymmesia) 


= ) dY ,(0) 
x Loy(0) ¥j(0) — - p2(0) —¥,(0)] 
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dY9,(¥o) dYon¥o) 


+ (Yale) - Youl¥e) 


dY,,(0) dY ,(0) 


ox) v5 0031]. (340) 


x (p4(0) ¥;(0) rs 


For constant speed of sound and ambient density, (340) becomes 


“ik. 9, D, =k, D 
ee ot a Ufo, -rk, ) 
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“ky D, “ik, D, *iky D, 


«(ogky -pak, De 27ace “sae 3 


+ (0, Ky, + po Ky ) (o; Ky, + p2 k,) 


jk D, +jk D, -jk D 
«(oak spat, De Yoglge. 937 © 30 


; (5, Ky, - p2 i ) leg ky, + po ky) 


-jk OD, -jk. D, +jk. D 
milpak -p3k Je Yo le Wa ore a 
: 


-jk D, +jk D, -jk. D 
(oak, + Pak Je Wy: le %3 2e ue Wy. (341) 


Using the definitions of Ro;, R23, and R34 presented earlier, (341) may 


be reduced to 
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+jK, 





-j2k, (D,- D,) 
lie aeesne een 2 n) 


~jk {D,- D,) -j2k. (D,- D,) -j2k Ci 
[Ros Rox e Me FN aL - Roy Rgge %3 ° - Ray Rage %2 | 


Using conditions (328) and (329), (342) may be reduced to 


ik ney “ik oy 
mean ies ee yo] / 


4nk 
Yo 


Arp, : 


-j2k. D 
| 


[1 ? Ro R3 € 7 (343) 


Equation (343) is equal to equation (256). Since (256) has already been 


verified, we can conclude that our derivation of the solutions for Aap: and 


hence, for Arp. for the four media waveguide problem is correct. 
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For the unknown constant Bo, Mathematica provided the following 


numerator: 


+ 


+ 


+ 


+ 


+ 


ca 


+ 


nums = G; (a).3 42,2 43.7 44.8 a5, 1 46.4 a7,6 
~ 41,2 823 427 448 454 864 276 + Ay 1 223 A374 g a5 0 ag 4 ATG 
~ Ay 1 822 37 gg a5 3 Ag 4 a7 6 - Ay 3 a2 2 AZ 4 agg a5 AG7 A765 


A, 2 42,3 43.4 Ag.g a5 26,7 A7.6 - 
At 42,2 43.4 Aq.g a5 3 G67 a7 6 - 
41,2 42,3 43.6 44,8 a5.) 46.4 a77 - 
At 1 42,2 43.6 Agg 45.3 464 a77 * 
~ 44,2 42.3 43.4 44g a5) 46,6 47,7 * 
Aj .1 42,2 23,4 A4.g 25,3 46.6 A772 - 
A, 2 42.3 43.7 84.6 45 4 AG 4 A789 * 
Ay 2 42,3 43.6 24,7 a5 1 AG 4 A7,8 - 
Ay 1 42,3 436 24,7 45 2 26.4 a7 4 * 
~ Ay 4 42,2 43.6 4,7 25,3 Ag 4 a7 sg - 
@} 2 42,3 23.4 44,7 25 1 466 478 - 
41,1 42,2 23.4 24,7 45,3 46.6 A798 * 
~ 81,2 42,3 43,4 24,6 a5 4 6,7 a7,8 * 
~ A141 82,2 a3 4 a4 6 5,3 67 a7 gl. 


A1,1 42,3 43.4 Ag.g a5 2 G7 a7 6 
81,3 42,2 43,6 44.8 45,1 46,4 87,7 
411 42,3 3.6 a4.g a5 2 AG 4 a77 
41,3 42,2 43,4 a4.g 45.1 46,6 877 
41,1 42,3 43.4 a4.g a5 2 26,6 877 
413 42,2 43.7 a4.6 a5) 46,4 97.8 
413 42,2 43.6 24,7 a5 1 26.4 a7 8 
A114 22,3 83,7 a4 6 45,2 46,4 a7§ 
411 22,2 43.7 44.6 45 3 46 4 A798 
41,3 42,2 83,4 A4.7 a5 1 46,6 87.8 
41 22,3 a3 .4 a4.7 a5 9 46.6 47.8 
Ay 3 82,2 A3,4 A4.6 AS | G7 a7 8 
41,1 42,3 43,4 44,6 45,2 46,7 47,8 


(344) 


Factoring (344) and collecting common terms yields the following 


generic expression for the numerator of Boy: 


hums = Gy lags 47,6 44,6 47,4) a3.4 46,7 - 43,7 a6.4) 
(ago a7.g - 248 a7 7) 43.4 46.6 ~ 43.6 46,4 
x lay» (ay as 3 ~ 413 as 1) + a23 (ay. 4517 411 as2)| 


(345) 


Substituting appropriate expressions into (345) and using (315) yields 


the following genera/ erpression for Boy: 
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- dY,(Dy) dY; (Do) 
Bop = On =o (D5) Y3(Ds) dy p4(D2) dy | Y 4(Do) 


dY5,(D,) dY,(D,) 


x p3(D,) ay Y3(Di) - p2(D,) oD) dy 


dY3(D,) dY (Dz) 


p4(D») Y,(D2) - p3(D2) Y;(D2) or 


dy 


dY 3, (Dy) dY3(Dj) 
De Sy LHD) - enlDy) Y39(D1) 


} 


dy 


dY ,(0) dY3,(0) 
x Ly3(v0) (o9(0) “ay Y2a'0) - p4(0) Y,(0) ewe 


dY,,(0) dY (0) 


+ Yoal¥o) (p\(0) ¥\(0) a) ee ie 


dy dy 


dY ,(D9) dY3 (D9) 


[I(o,00,) v0.) a palDa) ¥i(D)) 


dy 


dY5,(D,) dY,(D,) 


x Co3(D)) —— Y9(D1) - palDy) Yog(D) —— 
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dY3 (Da) dY (Da) 


+ (p4(Do) Y,(D2) - p3(D2) Y3(D2) 


dy dy 
dY5,(Dy) dY3(D,) 


¥3(D,) - pa(D,) ¥3,(0,) ———} 


dy 


AY on(¥o) dY2,(¥0) 


x [(Y3,(v0) 5 


Yon(Vo) 
dY,,(0) dY ,(0) 


Tay 7 P20) Val) 


x {o4(0) ¥,(0) ay 


dY2,(Vo) dY 54 (Vo) 


+ dy Yon(¥o) a Y4,(¥o) dy 


dY,,(0) dY ,(0) 


x Loy(0) ¥4(0) —F— - p20) —— ¥2,(0)5] 


dY4(Do) dY3(D2) 
- [ (03602) ¥3(D2) —— - pa(D2) —— Y4tD2) 


dY,,(D) dY3(Dy) 


x e3(D;) ey = e2(D;) Y op(Dy) dy 
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dY3 (D9) dY, (D2) 
+ (p4(Da) ay Ya(Da) - p3(D2) ¥3(D2) —— 


dY5,(D;) dY,(D,) 


Y3(D,) - po(D,) Yo,(D,) ——— 


x {03(D,) dy 


dy 


dYopl¥o) dYanty,) 


x (y3,(y,) dy ~ dy ava) 


dY,,(0) dY ,(0) 


x {p,(0) ¥,(0) —— - p9(0) ——¥;,(0) 
i dy dy 


Ye (y,) dY 54 (¥o) 
Yop(V¥o) ie Yoq(Vo) ae 





dy 


dY,,(0) dY ,(0) 


« (660) 410) —— - og) ~G— ¥,00} | (346) 


For constant speed of sound and ambient density, (346) becomes 


LS wha 
: 2 4 
Bob. a : (bs ‘ ise ky) 
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jk D, -jk D 
¥, 2 


-jk_ D 
[2 me "4 *T(o, Ky +99 k, ) Gs , - P2 Ky 


+ (o, ky * 0% ky) (5, ky * Pa k,,) 


sje. D, +jk D, -j. D 
«logk, resk, Je yo ewe aieeen os 
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= (a ky, - p2 k, ) lee ky, + p2 k,) 


Using the definitions of Rj, Ro3, and R34 presented earlier, (347) may 


be reduced to 


+jk rex aoe 0) 
es [Rog + Rage %3 | 


4xk 
Yo 

-j2k. sy +jk Oy -j2k. D 
Y5 le y, 0 Y, ! 


/ 


xile Rowe e 
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=j2k. (D.- D,) 
[Ros Rae € ce 


-j2t D, -j2k (D,- D,) 2s D 
- Ro; Rgge 2 We Y3 : ! - Ro; Rage 66 ale (348) 


Using conditions (328) and (329), (348) reduces to 


+jK +jk 
ee ore 


4xk 


y jk, ¥, 7) -i2k, D 
¥2 ° + Roe ¥2 | ¢ 9 / 


-j2k D | 


[1 - Roy R3 Be: (349) 


Equation (349) is equivalent to equation (261). Since (261) has already 
been verified, we can conclude that our derivation of the solutions for B..: 


and hence, for Bob. for the four media waveguide problem is correct. 


For the unknown constant A3, Mathematica provided the following 


numerator: 


num, = G; lay, 42.2 43.5 443 45,1 46,4 47,7 
~ 8,2 42,3 Az 5 A4g As) 464 AD 7 + Ay 223 43.5 agg a5 2 864 A77 
~ Ay 4 42,2 43,5 44.3 53 a6 4 AD - ay 3 a2 2 a3 4 agg a5 AG 5 a7 7 
+ 412 42,3 Az .4 Aq.g AS) 465 877 - Ay) 42,3 43 4 44g AS 2 AGS AD? 
7 A114 42,2 43.4 agg 45 3 AGS a7 7 ~ ay 3 AQ» a3 5 a4 7 a5 1 464 ADs 


172 


+, 2423 435 847 a5 1 26.4 a7 8 - Ay | 423 435 447 AS 2 4G 4 a7 8 
+, 1 422 435 a4 7 45 3 a6 4 a7 gt Ay 3 A292 AZ 4 a4 7 a5) 465 a7 
~ 41 2 42.3 43 4 a4 7 a5 1 465 47.8 t Ay 1 403 a3 4 a4 7 a5 2 AG 5 a7 4 
- Ay | 89.9 434 447 5 3 Ag5 78] (350) 


Factoring (350) and collecting common terms yields the following 


generic expression for the numerator of A3: 


numg = Gy (a4 a7,7 - a4,7 a7,5) eas ag4- 434 aes) 
x lap. (ay 4547 411 as, 3) + a2 3 (ai a5 2 - a1 2 as ;) 651) 


Substituting appropriate expressions into (351) and using (315) yields 


the following genera/ expression for Az: 


dY,(D2) dY ,(D.) 


-k - + 
Ag = 5.5 (pg(Dz) Y9(Dy) ay PalDa) a YG(De) 


dY 54(D1) dY 4 (Dy) 


x p9(D,) dy Yop (Dj) 5 92(D,) YD) dy | 


dY,,(0) dY ,(0) 


x Lys.tv) (5,(0) Y ,(0) anc p2(0) “ay Vaal) Ge, 


dY ,(0) dY,,(0) 


x (5,(0) “ay Y2al0) - p4(0) Y,(0) =) 
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[(o3(0, 


x te} Dye 


i p4(D2) 


x p3(D,) 


x (p4(0) ¥,(0) 


dY , (D2) 


)¥3(D2) —— - 


dY (D2) 


dY54,(D1) 


—orree N) . 


dy 


Y,(D,) - 


me (tes = 


dY3(D2) 


p4(D2) —Go— Yq(Da) 


dY (Dj) 
p2(D1) Y>,(D,) =a 


dY (Do) 
p3(Do) Y,(D2) dv 





p2(D,) Yon (Dj) 


dYon(Vo) AY 2,(¥o) 


Yon(ye) 


dY,,(0) 
dy 


dY,, (yo) 


dy Yon(Vo) 


si 


x (p,(0) ¥;0) —— 
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dy 


—a~ Yalve)) 


dY ,(0) 


- 920) ¥(0)) 


dY op (Yo) 


” Yoal¥o) ee) 


dY ,(0) 
- 620) 4 ¥54(0)}] 


dY (Do) dY3(D) 


7 [(o5(D2) Y,(D2) Pa 04D) 


Y,(D9) 
y 


dy 


; Y,(D,) 
dY5,(D,) d 3%] 


(D,) - >9(D)) ——} 
x p3(Dy) gy ¥3(Dy) e2(D,) Youl 1 dy 


dY (Dy) dY (Do) 


+ (p,(D,) Y, (D2) - p3(D2) ¥3(D,) Sales 


dy 


,(D,) 
dY5,(D,) dY( I 


Y3(D,) - p2(Dy) Ygg(D,) ———} | 


x Lp3(D,) dy 


dy 


dY 54 Vo) dY3,(¥o) 


x [(¥5,(y0) aye aie ano! 


dY,,(0) dY,(0) 


x {p,(0) ¥;(0) ae 


dY5,(¥o) de nlye) 
a 


+ dy Yas) = YoqlVo) dy 


dY,,(0) dY ,(0) 


: a 50) i (352) 
x {o4(0) ¥\(0) —— - o2l0) —¥3,(0)}] 
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For constant speed of sound and ambient density, (352) becomes 


ik. +fk D, -jk D 
hy ote %3 Toe (zoek, )losk, sesh, ) 


-jk OD 
[2 ky, @ erie ky, +02 k, ) (o; ky, = ky} 


ik, DY ly 
«(oak -eaky Je ome 3a ec mw 3 


+ (p, ky, +99 k, ) (5 Ky, + p92 k,,) 


+k, DY tik, D, ary, 


len Ky, + p3 k, Je 2 e 3 e 


7 (0, Ky, po ky ) (5, Ky, * p2 ky) 
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(353) 


Using the definitions of Ra;, Ro3, Rag, and Ty3 presented earlier, (353) 


may be reduced to 


at sik, oy jk yy -jk, D, «jk D 
Ay, - Tas le ¥2 °° + Roe Yo le Yo! eo %3 2 / 
4nk_ 
Be 
-jk, (D,- D,) «jk, (D,- D,) ik. Ds Dy ky 
[ass Ray e ue el ers i, - Roy R3q e : eee! e a 


Ro, R23 € Y3 Cc (354) 


Using conditions (328), (329), (354) may be reduced to 
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+jK, 41k Yo “ik, Yo) “it, D tik, D 
Ns coum o3n1© 2 oye meee le es aa 


-j2k. D 


[1 - Roy Ro3 o (355) 


Equation (355) is equal to equation (274). Since (274) has already been 


verified, we can conclude that our derivation of the solutions for A, and 


hence, for Ay. for the four media waveguide problem is correct. 


For the unknown constant B3, Mathematica provided the following 


numerator: 


Num = Gy [- a1 3 42.2 435 44.8 as | Ag 4 A765 
+12 423 435 agg as) 464 47,6 ~ 41,1 42,3 43,5 agg as 2 AG 4 a7 6 
+ 811 42,2 43.5 44g ag 3 a6 4 AD 6 + Ay 3 a2 2 AZ 4 agg a5) AG 5 ADS 
A1,2 42,3 43,4 44g 251 465 A765 + Ay) 42,3 a3 4 agg A520 aG5 a7» 
~ Aj 1 82,2 a3 4 a4 g As 3 AG 5 A76 + A, 3 a2 2 435 agg as | AG 4 a7 gs 
A} 2 42,3 a3.5 a4.6 a5 26.4 27,9 + Ay A253 435 246 45,2 464 a7 
41,1 42,2 43,5 24.6 a5 3 46 4 a7 - Ay 3 a2 2 434 446 51 AG 5 ADs 
+ 1,2 42,3 43.4 a4.6 a5 1 465 AZ g ~ a1 A223 AZ 4 a4 6 AS 2 AGS A798 


+ Ay 1 422 43.4 a4 6 5.3 465 a7 9] . (356) 


Factoring (356) and collecting common terms yields the following 


generic expression for the aumerator of B3: 
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numy ~ G, (a4¢ a7,5 - a4, 47,6) (43,5 46,4 - 43,4 a6) 
x lay» (ay 3 a5 4 - ay.) A553) + apg (ayy ago - ay a5 1)) (357) 


Substituting appropriate expressions into (357) and using (315) yields 


the following genera/ erpression for Bs: 





; dY (Do) dY;(Dy) 
k, + + 
Bs = 5" (p4(D2) Gy 7 Yale) - p3(D2) ¥3(D2) —— 
dYoy(D)) dY5,(D;) 
x p2(Dy) gy Yon(Dy) - 2(Dy) You(D1) Gy 
dY,,,(0) dY (0) 
x [¥5a(¥0) (61(0) ¥\(0) —— - 020) SG Yoa(0)) + Youle) 
dY (0) dY,,(0) 
x (210) — ¥3.(0) - 01(0) ¥,0) —) 
dY (Do) dY (D3) 
[I(e9(0) ¥500) —— p4(Dz) —— ¥4(Da) 
dY54(Dj) dY,(D,) 
x p3(Di) gy Y3(D)) p2(D1) Y2y(D)) a 
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dY (D2) dY ,(Dp) 
* loq(Da) Y, (D2) - 93(D2) Y3(D2) iw 


dY5,(D,) dY3(Dj) 
x p3(D;) ay Y3D1) = e2(D,) Y5p(Dy) ay 
dYon(¥o) dY2a¥o) 
x Yo,(¥o) dy = dy Yon(¥o) 
dY,,(0) dY ,(0) 
x {p,(0) Y , (0) ce p2(0) ae Y;,(0)} 
dYoa(¥o) —dY 25 (¥0) 
+ dy Yon(¥o) = YoalVo) dy 
dY5,(0) dY ,(0) 
“ {4(0) Y , (0) ae p2(0) Se ll 
dY, (D2) dY3(D,) 


- [(ostD) ¥5D2) —— - 4(D2) —e— Y4(D2) 


dY>,(D,) dY,(D,) 
x p3(D,) ey 3D) B p2(D,) Yop(Dy) ae 
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dY,(D2) dY ; (Dy) 


+ p4(D2) ay Ya(Da) z p3(D>) Y,(D2) Sve e 
dY5,(D,) _—— d¥, (Dj) 
x {03(D,) aay Yau) 5 p2(D,) Yop(Dy) dy 


Volvo) dNsduve) 


x [(¥3,(¥) =a Yale! 


= dY,(0) 
x {p,(0) ¥;(0) —— - p90) —— ——¥p,(0)) 


dY5,(¥o) AY 54 (Vo) 
+ ‘evay Yon(Vo) is Yo,(Vo) dy 


dY5,(0) dY (0) 


x {p,(0) ¥,(0) - p9(0) =—Y5,(0 od]. (358) 


dy dy 


For constant speed of sound and ambient density, (358) becomes 


jk, city Dy ik, D, 
Byopee 73 84 (Zork, ) lock, esky) 
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-jk_ OD 
[2 ky e V4 ae Ky, + p9 k, ) (o, Ky, - 92 k,) 


= (5, Ky - p2 Si ) (0; ay - p2 ky) 
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-jk D, +k D, -ik. D 
x (p4 k, *p3k, ye ete ae es iy (359) 
3 4 


Using the definitions of R24, R23, R34, and T23 presented earlier, (359) 


may be reduced to 


+jk tik) oy oes iat | Mie ae | 
By = Rya Tas le 17 + Ron € ¥, le uP le 73 oof 


4xk 
a3 
-jk, (D,- D,)  +ik, (D,- D, ) Fir e(Ds= Di) (ok aD 
[Ros Reve y, 2 | es aes | Re Ryae Yel ey, | 
sjk (D,- D,) -i2k. D 
- Roy Roz 3 : € a) a, (360) 
Using conditions (328), (329), (360) may be reduced to 
B, -0. (361) 
C 


Equation (361) is exactly what one would expect for a semi-infinite 


medium (i.e., no wave propagation in the negative y direction). Thus, we can 
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conclude that our derivation of the solutions for B,. and hence, for Bs for 


the four media waveguide problem is correct. 


For the unknown constant B,, Mathematica provided the following 


numerator: 


num, = G; [a1 3 a2.5 43.7 44.8 a5 9 ag.4 a7,6 


41.2 42,5 43,7 44.9 a5 3 a6 4 a7.6 - 
+ A, 2 42,4 437 agg a5 3 AG 5 a7 6 - 
A) 3 42,4 43,5 a4.g As 2 467 a76 + 
~ 81.2 42,4 83.5 a4 a5 3 467 A756 - 
+ @1 2 42.5 43.6 a4.g as 3 ag 4 a7 7 + 
~ 41,2 42,4 43,6 44g a5.3 a6.5 a77 * 
~ 41,3 42.4 435 24.8 a5 2 46.6 877 - 
+ 1,242.4 43.5 44g 45 3 46.6 7,7 - 
+ 1,3 42,5 43.6 44,7 25,2 46 4 a7, + 
— 81,2 42.5 43.6 84,7 a5 3 46.4 A798 + 
~ 41,3 42,4 436 24,7 a5 2 465 47.8 - 
+ Ai 2 82,4 43.6 44,7 a5 3 a6 5 a7 8 - 
+ 13 42.4 43,5 44.7 a5 2 46.6 7,8 + 
41,2 42,4 43,5 44,7 a5 3 6.6 a7 9 + 
413 42,4 43.5 a4.6 a5 2 967 a7 8 - 


+ 412 42,4 435 44.6 a5 3 467 a7 8) . 


+ 


41,3 42,4 43,7 a4g a5 2 465 a7 6 
A 3 42,5 43.4 a4.g as 2 467 876 
Aj 2 42,5 43.4 44.8 a5 3 a6 7 a7 6 
41,3 42.5 23.6 44.8 a5 2 46 4 a7 7 
413 42.4 43.6 44.8 a5 2 46,5 a77 
41,3 42,5 43.4 44g 45 2 46.6 47,7 
@1 2 42,5 43.4 a4. a5 3 a6 6 a77 
41,3 42,5 43,7 44.6 45 2 46,4 47.8 
@1,2 42,5 43,7 44.6 5,3 46.4 a7 8 
41,3 42,4 43.7 44.6 45 2 46.5 a7 8 
412 42,4 437 a4 6 a5 3 AG 5 a7 8 
13 42,5 43,4 44,7 a5 2 26.6 a7,.8 
@1 2 42,5 43,4 84.7 a5 3 466 a7 8 
41,3 42,5 43,4 44.6 a5 2 ag 7 a7 8 
41,2 42,5 43.4 44.6 45,3 467 a7 8 


(362) 


Factoring (362) and collecting common terms yields the following 


generic expression for the numerator of B,: 


num, = G, (ay 3 as 2 - a1 2 as 3) lays [fag¢ a7g- 448 ay 6) 


x {a3.4 467 ~ 43.7 a64J + (ass a77 - 


247 27,8) 13,4 26,6 - 23,6 96,4 


- 424 (ase a7g - 448 a> 6) 435 467 - 43.7 a6 5 
5: (agg a7.7 ~ 447 ays) 435 46.6 ~ 43,6 465/11 - 
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(363) 


Substituting appropriate expressions into (363) and using (315) yields 


the following genera/erpression for By: 


dY,,(0) dY,,(0) 


By ~ 5° (020) —— Ya(0) - 02(0) Yzq(0) —— 


dY, (D>) dY (Dp) 


YD) o0s es) 


x Lysa(ve) I(o4(D.) dy 


dy 


dY5,(D,) dY,(D,) 


aye Ya(D)) ‘s p2(D,) yo, (0p 


dy 


dY, (D2) dY; (D9) 


+ (03(Dp) ¥4(Dp) erry 


dY5,(D,) dY (Dy) 


x {p3(D,) Y3(Dj) - p2(D,) Yop (Dy) ae ee 


dy 


dY (Dy) dY ,(D,) 


dY¥5,(D,) dY,(D,) 


“gy YD) i e2(D,) Yo(D1) 





x p3(D,) 
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dY (Do) dY (Do) 


+ (p3(D3) ¥3(D2) —— - p4(D2) 


Y (D2) 


dy 


dY54(Dj) dY,(Dj) 


e(Dy) - p2(D,) Y op(Dy) ——}]] / 


x {3(D,) y 


dy 


dY; (D2) dY, (D9) 


[1 (e502) Y3(D) (D2) ay Y4lDa) 


or 


dY54(Dj) dY,(D;) 


x p3(D,) dy Y,(D,) 2 p2(D,) Yop ——. 


dY, (D2) dY;(D.) 


+ \p4(D9) ay Y4(Da) - 93(D.) Y3(D) dy 


dY 54 (Dj) dY (Dy) 


x {o3(1) —go— 4D) - 621) Y3q(D) ——}] 


IY o4(Vo) IY z4(¥o) 
x (sly, ram mare 


dY,,(0) dY ,(0) 


x {p,(0) 0) Y,(0) ~—- p20) 


i Yy(0) 


dy 
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dYouly) dY 54 (Yo) 


SNe gee ong) = foal Vo) 


dy dy 


dY,,(0) dY ,(0) 


x (p4(0) ¥;(0) ie p2(0) —— ¥3,(0)) 


d¥, (D2) dY3 (D9) 


= {(p3(D2) Y3(D) cia p4(D4) ay Y4(D2) 


AY 99(Dj) ; dY,(D;) 
x p3(D,) dy Y3(Dj) = e2(D,) Yop(Dy) Tay 


dY3(D,) dY (Dy) 


Y,(D2) - p3(D2) Y;(D2) —~—— 


+ \p4(D9) dy 


dy 


dY (D1) dY,(D,) 


Y3(D,) - pg(D,) ¥54(D,) ——}] 


dy 


dvon(yok diay, 


x [(¥3,¢y,) iogee en Yipl¥o)) 


dY,,(0) dY ,(0) 


x (64(0) ¥;(0) — - p2(0) 


: Ypq(0) 


dy 
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+ (~~ Yaplvo) - Yaal¥o) Gy 


~ dY,(0) 
x {p(0) ¥,(0) —— - p20) —— =n (364) 


For constant speed of sound and ambieny density, (364) becomes 


ik jk, D 


ie 2g (20k, de *4 


2 


“ik, Dy tik, y, -iky Dy tik, D 


«Mosk, serky Vest, -eak, De 2 € 2 e€ 3 e€ 3 


ik, dD “ik, ya tik, D, , 


-jk. D 
a (03 ky, — Po k, (0, Sy + Pq ky Je c e . e 3 oe 


*ik, D, “Iky ue “ik, D, , 


+{k D 
Mosk, -ork, )lesk, -psky }e 2 ¢@ ‘2% e @3 “mee 
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189 


ae Dajte (Ductal D 
7 ee en 


x (og k, rpsk, de 2 te e (365) 
3 4 


Using the definitions of Rg;, Roz, Ra4, and To; presented earlier, (365) 


may be reduced to 
+1ik 
a |Kr pe 


le 
4nk 
” 





jk oy -j2k D, -j2k (D,- D,) -j2k_ D 


-jk, y aig De De) 
en? [Ne Reina camer Vl / 


apie (i eed) -j2k (D,- D,) -j2k D 
lene % + 1 - Ro; Rage Y3 : ! e 72 


-j2k dD 
Ro Ro3 e ip) ] ; (366) 


Using conditions (328) and (329), (366) reduces to 
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+jk Sicko se ik 7 y. qlee 
B, 3 : To [Ro Ee 2 re ae ee) / 
c 


4nk 
by, 
-j2k_ D | 


[4 = Roy Ro e — (367) 


Equation (367) is equal to equation (268). Since (268) has already been 


verified, we can conclude that our derivation of the solutions for Bi. and 


hence, for B, for the four media waveguide problem is correct. 
¢ 


For the unknown constant Ay, Mathematica provided the following 


numerator: 


nums = G, ae 2,2 43.5 A472 a5) 464 A765 
— 412 42.3 43.5 a4 7 a5) 26 4 A76 + Ay) 423 a3 5 a4 7 a5 2 AG 4 a7 6 
41.1 42,2 43,5 44,7 453 46.4 476 - A 3 a2 2 A354 Ag 7 a5 AG 5 a7 6 
+ Ay 2 42.3 43 4 44,7 a5 1 46,5 476 ~ Ay 1 42,3 AZ .4 Ag 7 As 2 a6 5 a7 6 
* 11 42,2 43.4 84,7 45.3 465 a7 6 ~ Ay 3 a2 2 43,5 agg as AG 4 aD 7 
+ Ay 2 423 43.5 a4 6 a5) 464 477 - Ay) 42.3 AZ 5 446 a5 2 AG 4 A772 
+ 411 42,2 43.5 44,6 453 46.4 A772 + Ay 3 AD 2 Ag 4 ag 6 a5) AG 5 877 
Ay 2 42.3 43 4 44 6 aS) 65 477 + Ay ADs a3 4 a4 6 a5 2 AG5 A727 


= ay} a2 2 a3 4 a4 6 as 3 a65 a7 7 : (368) 


Factoring (368) and collecting common terms yields the following 


generic expression for the numerator of Ag: 
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numg ~ Gy (a47 a7,.6 - 44,6 47,7) a5.4 46,5 ~ 43,5 are 
x laa (ay 45,3 - 44,3 85,1) - 22,3 Cars 5,2 - ayo a5,1)1-. (369) 


Substituting appropriate expressions into (369) and using (315) yields 


the following genera/erpression for Ag: 


r dY3(D,) d¥ (D2) 


- — (p,(D,) ay ¥3(D2) - p(Da) Y3(D2) —— 


Aq 2 y 


dY 54 (Dj) dY5,(Dy) 


x (patDy) YD) } - p2(D)) —— You (Di) 


dy 


dY,(0) dY,,(0) 


eer | ) (p(0) —— yo taal) - p,(0) Y,(0) = 


dY,(0) dY,,(0) 


- Yql¥o) (p9(0) So ¥q(0) - py (0) ¥4(0) Sl / 


dY (D2) dY3(D) 
[[(c4(0.) Y,(D,) aera p4(D2) mira: Yi(D2)) 


dY 54 (Dj) dY3(D,) 


x e3(D,) mere lO i p2(D,) sekh), ee 
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dY3 (Da) dY{(D,) 


+ p4(D2) dy Y, (D2) = p3(D2) Y,(D9) _ 
dY5,(D,) dY,(D,) 
x p3(D)) gy Y3(D)) - 92(D1) Yo9(D1) Gy 
dY5,(¥,) dYoely,) ; 
| Valve) es YoplVo) 
dY,,(0) d¥,(0) 
x {4(0) Y , (0) ae p2(0) dy Y2al0) 
dY 5 (Yq) dY a4 (Vo) 


2 eerior Yapl¥o) - Yoal¥o) ae aie 


dY,,(0) d¥,(0) 
x {p,(0) Y,(0) F200) Te ON 
dY, (Dp) dY3(D,) 


nL [(os(D.) Y, (Dy) dy = p4(D2) ae Y;(D,)) 


dY (Dy) dY (Dy) 


x e3(D,) Y,(Dj) = p2(D,) Yop(D1) “dy 


dy 


13 


dY, (D9) dY (Dz) 


* \p4(D2) —Go— Y,(D2) - p3(D2) Y,(D2) aan 


dY>4(Dj) dY,(D) 
ay Y3(D) - p2(Dy) Yop (Dy) 


}] 


x 1p3(D,) “ty 
dYon (Vo) dY3,(¥o) 

x onl) dy a tee Yon(Yo) 
dY,,(0) dY (0) 


x Loy(0) ¥4(0) —— - p20) —— Yoa(0) 


d 


dY,,(¥o) dY 54 Vo) 


+ a Youve) = Y5, (yy) = 


dY,(0) dY ,(0) 
x {,(0) Y ,(0) a p2(0) Sell (370) 


For constant speed of sound and ambient density, (370) becomes 
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jk. D 
[2 k,, ea (pi ky, +p k, } (o, k,, - py ky} 
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: (5, Ky, - P2 k, ) (5, Ky, a2 k,,) 


“ky D, tik, D, ey \| 


SM p sk) meet (371) 


Using the definitions of Ry;, Roz, R34, and To, presented earlier, and 
defining the transmission coefficient at the boundary between medium three 


and medium four as follows: 





2 p3 Ky, 
134 =, (372) 
P4 <A + p3 Ky 
(371) may be reduced to 
+jk 

Fe meee lene 

4nxk 

Y9 


-jk Yo ier D, 


rik, y -jk. (D,- D,) +jk, D 
alle 12 Ree ewe See ON 


2] / 


€ 
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-i2x. (D,- D,) eke D.  =f2k 2D, = De) 
2 ae | 
Lao; Rye °3 Sele Ra Rage 2 3) 24 


-i2k, >] 
~RopRoge 2 0a. (373) 


Using conditions (328), (329), and letting 


T34=1, (374) 
and, as aresult, 
Ky = ky (375) 
(373) reduces to 
+ik jk oy ak, ¥,) “tt, D +k, D 
RE = Us gene +*Ry,e 2 Olena ve ia a 
4xk 
Yo 


-j2k D 
| (376) 


[1 -Ra,Ra3e 2 


Equation (376) is equal to both equations (275) and (355). Since (275) 
has already been verified, we can conclude that our derivation of the 
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solutions for A ne and hence, for A ae for the four media waveguide problem is 


correct. 


To summarize the results for the four media waveguide with plane, 
parallel boundaries, for the general case, in which speed of sound and 
density are arbitrary functions of depth, the unknown constants are given 
by (322), (333), (340), (346), (352), (358), (364), and (370). As in the three 
media waveguide case, the denominators in these expressions are exactly the 
same. For constant speed of sound and density, the unknown constants are 


given by (327), (335), (342), (348), (354), (360), (366), and (373). 


Having completed two successful tests of the programming technique, 
we now return our attention to the general waveguide problem. As in the 
simpler examples, the compact vector- matrix system equation (198) applies. 
For the general waveguide problem, the vector matrix quantities involved 
are as follows: 


¢ A is the 28 by 17 matrix of coefficients, 


A- [0 a; > a;3a,;40000000000000] 
[az,;000a25a2,00000000000) 
[0 0 a3 343,000 a,5a,,00000000) 
10000 a5 a46000 a4 19 a4 1,000000) 
[0000000 asg agg 0 0 as 12 as 13.000 0) 
(000000000 ag 49 461100 a6 44 a¢15 0 0) 
(00000000000 ay 42 a7 430 0 an 46 0) 
LO000000000000 ag 44 ag 15 0 ag 17] 
[0 ag 2 ag 3. agg 0000000000000) 
[ayo,1 000 ayos ayo, 00000000000) 
[0 ayy 2 ayy 3 4, 49000000000000] 
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where 


laya 1 0.00 ayo5 a2400000000000) 

[ay3 000 ay35 a43600000000000) 

(0 ay40 444.3 41440000000000000) 

[0 000000 415.8 415.9 00 415,12 215,13 000 0] 
[0 00000000 416,10 216,11 00 416,14 416,15 0 0] 
j0000000 417.8 217.9 00 417,12 417,13 0000] 
{0 00000000 418,10 aigit 00 418,14 418,15 0 0) 
[0 00000000 41910 419.11 00 419.14 419.15 0 0] 
{0 000000 420.8 420.9 00 420,12 490,13 000 0} 
[0 0000000000 421,12 421,13 00 491,16 0] 
(Oo0O00000000000 492,14 422.15 0 20.17) 

[0 0000000000 423,12 223 13 00 293 16 0) 
lo000000000000 424,14 424 15 0 24.17] 

[0 000000000000 495 14 495,15 0 425 171 

[0 0000000000 49612 426,13 00 426,16 0] 

{0 0 427,3 427.4 000 497 § 427.9 00000000) 

[0 000 4228 5 428 6 =10'0 428,10 428.11 00000 0)] 


12 = - pylys) In(k, 0) Y (ys) 
41.3 = P2l¥s) In(k, 0) Yoal¥s) 
41.4 = p2l¥s) Ink, ae Yo(¥s) 
a2.1 =~ pilys) In(k, 0) Yy(ys) 
a2.5 - p2lYs) In(k,0) Y5a(¥s) 
42,6 = p2lys) Jnlk,,0) Yo4(¥s) 


a337 Y5a(Vo) 
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a3.4 = Yo,(¥o) 
a3,5 ~~ Yop(Vo) 
43,9 - - Yop(¥o) 

abe. = Ypa(Vo) 

44.6 = YoalVo) 
4410 ~~ Yop(Vo) 
44.11 = ~ Yop(¥o) 


45.8 ~ 02(¥p,) Ink, o£) Yon(Vp,) 
a5,9 = P2l¥p, ) Jalk,0) Yoo(Vp,) 
45,12 = - e3(¥p,) Inlk, 0) Y3(¥g,) 
45,13 ~~ 3p.) Jn(k,,t) Y3(¥p,) 
46,10 = P2l¥p,) Jalk,,0) ‘GAWE e 
46,11 = P20Vp,) Jnlk, 0) Yon(¥p,) 
4614 7 ~ O3(¥p,) JalK,,0) Y3(¥p,) 
46,15 =~ 03(¥p,) Jn, 0) Y3(Yp,) 


47,12 = P3(¥py) In(k,0) Y3(¥p,) 
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(384) 
(385) 
(386) 
(387) 
(388) 
(389) 
(390) 


(391) 
(392) 
(393) 
(394) 
(395 
(396) 
(397) 
(398) 


(399) 


47,13 = P3(¥g,) In(k,_f) Ya(¥p,) 
47,16 = ~ P4(¥p,) Jn(K,F) a 
ag,14 = 03{¥p,) Jalk,,1) Y3(¥p,) 


48,15 = O3(¥po) In(k,t) Ya(¥p,) 


ag.17 =~ P4 (yp) Salk, P)Y%p,) 
dY (ys) 
ago =- Jak, 1) ao 
dY5,(¥s) 
dga = Jnk,0) dy 
dY,(ys) 
= J,tk. st) Sa 
dY (ys) 
aio = - Ink, 0) Tite 
dY,,(¥s) 
410.5 = Jn{k,°) ae 
dY,,(ys) 
410.6 - Jalk, r) dy 
Aalke, r) 
414.2 =- ke, ——— Yy(ys) 
d(kr, r) 


(400) 


(401) 


(402) 


(403) 


(404) 


(405) 


(406) 


(407) 


(408) 


(409) 


(410) 


(411) 


11,3 = Kr, Y>,(¥s) 
fo r) 
dJalkr, r) 

41147 Kr, Yo,(¥s) 
fo r) 


qi2, tae Jn{ke, r) Y;(ys) 
i255 JalKe, r) Y54(¥s) 


412.6 = Jalkr, r) Y5,(¥s) 





Dalke, r) 
131-7 Ke, Y (ys) 
d(ke, 1) 
dJalkr, r) 
413.5 = Kr, Y2,¥s) 
i) r) 
dake, r) 
13,6 = Kr, Y5,(Ys) 
ry r) 


414.2 = 7 In{kr, r) Y (ys) 


414.3 7 Jntkr, r) Y54(Ys) 
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(412) 


(413) 


(414) 


(415) 


(416) 


(417) 


(418) 


(419) 


(420) 


(421) . 


414.4 = Iq(Ke, 0) Yoal¥s) 








dY 24 (¥p,) 
a5. = Jnk,,0) dy 
dYon(yp,) 
a15.9 = Ink, 0) dy 
dY3(¥g,) 
415,127 ~ In(k,,r) dy 
dY (yp) 
415,13 ~~ In(k, .F) dy 
dY 25(¥p,) 
a = 
16,10 Jalk,,0) dy 
dYon(¥g,) 
a = 
16,11 Ink, f) dy 
dY3(¥p,) 
416.14 =~ Ink,0) dy 
dY3(yp,) 
416,15 == Jnfk,,F) dy 
dJ,{ke, ) 
er ee 
17,8 ~ Ke, Yonl¥p,) 
ee) 
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(422) 


(423) 


(424) 


(425) 


(426) 


(427) 


(428) 


(429) 


(430) 


(431) 


dJq(kr, r) 


417.9 = Kr, ae aa Yoo(¥p,) 
Po r 
alkr, r) 
Ia) ok kp, ——— Yolyq,) 
d(kr, r) 
dJn(Kr, r) 
ais = Kr, eee ap) 
kr, r) 


418.10 7 In{ke, 0) Yoo(¥p,) 
A1g.11 7 In(ke, 7) Yoo(¥p,) 
aig.14 ~~ In(ke, 0) Y3(yp,) 
418,15 ~~ Ialke, t) Ya(yp,) 
dJalkr, r) 
419.10 = Kr, —_—_—~ Yarl¥p,) 
f> r) 
dJalkr, r) 


We etry, 
fo r) 
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(432) 


(433) 


(434) 


(435) 


(436) 


(437) 


(438) 


(439) 


(440) 


419,14 =~ Kr, ) Y3(¥p/) 
d kr, r 
dJalkr, r) 

419,15 =~ Kr, Ya¥p,) 
d(Kkr, ¢) 


420.8 = Jake, 0) Yon(¥p,) 
420,9 = Jake, 0) Yon(¥g,) 
420,12 ~~ Jn{ke, 0) Ya(¥p,) 
420,13 ~~ In(ke, t) ¥3(¥g,) 


dY3(Vp ) 


42112 = In{k,.0) Gy 


dY3(yp,) 
421,13 ~ In(k,,t) Go 


dY4(¥g,) 


421,16 =~ In{k, 6) Gy 


dY3(¥p,) 
422,14 = Inlk, 0) 


dY3(Yp,) 


422,15 = Jak, 0) ay 
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(441) 


(442) 


(443) 


(444) 


(445) 


(446) 


(447) 


(448) 


(449) 


(450) 


(451) 


dY4(¥p,) 


422,17 =i Jatk, 0) ay 


CJalkr, r) 
423,12 = ke, —_— Ya Yg,,) 
s d(ky, t) : 
dJnlkr, r) 
23,13 = kr, ~~ Yal¥,) 
r, 0) 
dalkr, r) 
423,16 =~ Ke, Ya(¥pp) 
d{kr, r) 


424,14 7 Jn{Kr, ©) Y3(¥p,) 
424,15 ~ In(kr, 0) Yal¥p,) 
424,17 ~~ Inlke, ©) Y4(Yp,) 


dja kr, r) 


el een) 
fr, r) 


O]alkr, r) 


22515" key Yala, 
rs r 
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(452) 


(453) 


(454) 


(455) 


(456) 


(457) 


(458) 


(459) 


(460) 


d)q(Kr, r) 


425,17 = - kr, ———— Yal¥g,) 


d(ke, r) 


426,12 ~ Inlke, 2) Y3(¥p,) 
426,13 7 Jntkr, r) Y3(¥p,) 


426,16 ~~ Ja(Ke, r) Y4(¥p,) 


AY 5(¥o) 
427.3 = dy 
dY o4¥o) 
427,4- dy 
dY on(¥o) 
927.8" dy 
dY o4(¥o) 
427.9 =~ dy 
aia) 
4285 a dy 
dY9,(¥o) 
49% 6 = dy 
dY 95 (Yo) 
428,10 =~ ~ dy 
dY op (Vo) 
2841 = 7 ody 


and 
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(461) 


(462) 


(463) 


(464) 


(465) 


(466) 


(467) 


(468) 


(469) 


(470) 


(471) 


(472) 


¢xis the 17 by one vector of unknown constants, 


x = (A; By Apa Bog Coa Dog Gz Arp Bat Cop Dap Az B3C3D3 Ag Ba)’, (473) 


where the superscript T indicates the transpose matrix operator (which 


means that x is a column vector). 


The vector of known constants, b, is simply made up of the right-hand 


sides of (149) through (176), or more specifically 


b-(00000000000000000000000000G,0}', (474) 


where here again the transpose operator is used to indicate that b is a 


column vector. 


One should note that the arbitrary constant G, appears in the known 
constant vector while the arbitrary constant G) appears in the unknown 
constant vector. This occurs because G,; represents the known amount of 
discontinuity required to achieve the free-space Green's function solution 
under the necessary conditions for that solution to exist (i.e., constant speed 
of sound, no boundaries). Go, on the other hand, is really an artifact of the 
method used to derive these boundary condition equations, and, as such, 


should be treated as an unknown quantity in the general case. 


In the general waveguide problem, the form of matrix A leads to some 
complications. First, the fact that the number of equations (i.e., 28) is greater 
than the number of unknowns (i.e., 17) implies that the solution to this ~ 


problem will not be unique. Secondly, the fact that the matrix A is not a 
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square matrix implies that the simple matrix inversion technique of (222) 
cannot be utilized. Thus, an alternate solution technique must be employed 


to set up the problem so that Afz7/(fematica may be used in the solution. 


Gelb (1974) describes this situation (more equations than unknowns) as 
an overdeter mined case. Gelb and Haykin (1986) both suggest the use of a 
pseudoinverse matrix in obtaining a least squares estimate for the vector x. 


The pseudoinverse matrix is defined by Gelb as follows for real matrices A: 


A°=(ATA)! aT. (475) 
where A° is the 17 by 28 pseudoinverse matrix. 


Haykin and Menke (1984, p. 253) both define a similar pseudoinverse 


matrix for complex matrices A as follows: 


A*’=(AH A)! al (476) 


where the superscript H indicates the Hermetian or complex conjugate 


transpose matrix operator. 


In the solution for the general waveguide case, we will use a 
combination of these pseudoinverse matrix techniques. This combination 
incorporates a 28 by 28 weighting matrix, W, which allows us to obtain a 
weighted least squares estimate for the vector x. Therefore, our approach 
will be to use the weighted pseudoinverse formulation suggested by Gelb 


and Menke (1984, p. 54) with the complex conjugate transpose operators 
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suggested by Haykin in lieu of the standard (and less general) transpose 
operator. Thus, the pseudoinverse matrix to be used in this thesis is defined 


as follows: 


A°-(AHWA)! AHW. (477) 


Using this pseudoinverse matrix, a weighted least squares estimate for 


the vector x may be calculated as follows: 


I-A’b. (478) 


The following Afatfhematica code was developed to solve the general 


waveguide problem using this weighted least squares technique: 
b = (0,0, 0, 0, 0, 0,0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0,0, 0, 0, 0, 0, 0, 0, 0, G1, 03; 
w = {({wicl, wlc2, wic3, wlc4, wlc5, wlc6, wic7, wic8, 
wic9, wilcl0, wicl1l, wlcl2, wlcl3, wlcl4, wlcl5, 
wicl6, wicl7, wicl8, wicl9, wic20, wlc21, wlc22, 
wic23, wlc24, w1c25, w1c26, w1c27, w1c28}, 
(w2cl, w2c2, w2c3, w2c4, w2c5, w2c6, w2c7, w2c8, 
w2c9, w2c10, w2cl1, w2c12, w2c13, w2c14, w2cl5, 
w2cl6, w2cl17, w2c18, w2c19, w2c20, w2c21, w2c22, 
w2c23, w2c24, w2c25, w2c26, w2c27, w2c28}), 
(w3cl, w3c2, w3c3, w3c4, w3c5, w3c6, w3c7, w3c8, 
w3c9, w3cl0, w3cl1, w3cl2, w3cl3, w3cl 4, w3cl5, 
w3cl6, w3cl7, w3cl8, w3c19, w3c20, w3c21, w3c2Z, 
w3c23, w3c24, w3c25, w3c26, w3c27, w3c28), 
(w4cl, w4c2, w4c3, w4c4, w4c5, w4c6, w4c7, w4c8, 
w4c9, w4cl0, w4cl 1, w4c12, w4c13, w4c14, w4cl15, 
w4c16, w4cl7, w4c18, w4c19, w4ce20, w4c21, w4ce22, 
w4c23, w4c24, w4c25, w4c26, w4c27, w4c28}, 
(w5cl, w5c2, w5c3, w5c4, w5c5, w5c6, w5c7, w5c8, 
w5c9, w5cl0, wScl1, w5c12, w5cl3, w5cl4, w5cl15, 
w5cl6, w5cl7, w5c18, w5cl9, w5c20, w5c21, w5c22, 
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woc23, w5¢e24, w5cZ5, w5c26, w5c27, w5cZ8}, 

(w6cl, w6c2, w6c3, w6c4, w6c5, w6c6, w6c7, w6c8, 

w6c9, w6cl0, w6cl1, w6c!2, w6c13, w6c14, w6cl5, 

w6cl6, w6cl7, w6cl8, w6cl9, w6c20, w6c21, w6c22, 

w6c23, w6c24, w6c25, w6c26, w6c27, w6c28), 

(wicel, W/eZ, w/C3, W/c4, W705, W/C6, W/c/, w/c8, 

w7c9, w7cl0, w7cll, W7c12, w7c13, W7c14, w7cl5, 

w7cl6, w7cl7, w7c18, w7c19, w7c20, w7c21, w7c22, 

w7c23, w7c24, wW7¢C25, W7c26, W7C27, W7c28), 

(w8cl, w8c2, w8c3, w8c4, w8c5, w8&c6, W8c7, w8c8, 

w8c9, w8&cl0, w8cl1, w8c12, W8c13, W8c14, W8cl15, 

w8cl6, w8c17, w8c18, w8c19, w8c20, W8&c21, W822, 

w8c23, w8c24, w8c25, w8c26, W8c27, w8c28}, 

(w9cl, w9c2, w9c3, W9c4, W9CS, W9C6, W9C7, W9C8, 

w9c9, w9cl10, wW9c1 1, wW9c12, W9c13, W9C14, W9C15, 

w9cl6, w9cl7, w9cl8, W9c19, W9C20, W9C21, W9C22, 

w9c23, W9c24, w9C25, W9C26, W9C27, W9C28}, 

{wl0cl, wl0c2, w10c3, w10c4, W10c5, wW10c6, W10c7, w10c8, 
wl0c9, w10cl0, w10cl1, wl0c12, W10c13, W10c14, w10c15, 
w10cl6, w10c17, w10c18, W10c19, W10c20, W10c21, W10c22, 
w1l0c23, w10c24, w10c25, W10c26, W10c27, w10c28}, 

fw llelewlic2, wl tc3, wil ic4, wl lcs, wl lc6, wl lc7, wl 1c8, 
wllco wliclOewilecl ll wiici2, wiltel3, wlici4, wl lel5. 
wlicl6, W1lc17, wllcl8, wl lel9, wl ic20, w11e21, wl le2z, 
wl1c23,wl1c24, w11c25, w11c26, w11c27, w11c28), 
(wl2cl, w12c2,w1l2c3, w12c4, wl2c5, w12c6, w12c7, w12c8, 
w 269; wizclU, w lizcl i, wizel2, wizcel4.wl2ci4 wizcl5, 
wl2cl6, w12cl7, w12c18, w12c19, w12c20, w12c21, w12c22, 
w i2e23, WiZe24, w12025, wizc26, wl 2¢27, wl 2cZ25), 
(wl3cl,w1l3c2, w13c3, w13c4, W13c5, W13c6, W13c7, w13c8, 
wiscd, wl4acl0,wil3cll, wiselZ, wil 3cl3, wl3el4, wl 3cl5, 
w13cl6, w13c17, w13c18, W13c19, w13c20, w13c21, w13c22, 
w13c23, W13c24, W13c25, w13c26, W13¢c27, w13c28}, 
(wl4cl, wl4c2, w14c3, w14c4, W14c5, W14c6, W14c7, w14c8, 
w14c9,w1l4cl0, wl4cl1, w14c12, w14c13, w14c14, wl4cl5, 
wl4cl6, wl4cl7, w14c18, w14c19, w14c20, w14c21, w14c22, 
w14c23, W14c24, W14c25, W14c26, W14c27, w14c28}), 
(wl5cl, w15c2, w15c3, wW15c4, W15c5, W15c6, W15c7, w15c8, 
wl5c9, w15cl0, wI5cl1, w15c12, w15c13, w15cl14, w15cl5, 


w15cl6, w15c17, W15c18, w15c19, W15¢20, w15¢c21, w15c22, 
w15c23, w15c24, w15c25, w15c26, w15c27, w15c28), 
(wl6cl, w16c2, w16c3, w16c4, w16c5, W16c6, Ww 16€7, w16c8, 
wl6c9, w16cl0, W16cl1, wW16c12, w16c13, W16c14, w16c15, 
w 16c16, w16c17, w16c18, W16c19, w16c20, W16c21, W16c22, 
w 16c23, w16c24, w16c25, w16c26, w16c27, w16c28), 
(wil7cl, wl7c2, w17c3, w17c4 wl7c5, wW17c6, W17c7, w17c8, 
w17c9, w17cl10, W17c11, W17c12, W17c13, W17c14, w17cl15, 
wl7cl6, w17cl17, W17c18, W17cl19, w17c20, W17c21, w17¢22, 
w17c23, w17c24, w17c25, w17c26, w17c27, w17c28), 
(w1l8cl, wl8&c2, w18c3, wW18c4, W18c5, W18c6, W18c7, W18c8, 
wl8c9, w18cl0, W18cl1, W18cel2, W18c13, w18cl4, w18cl5, 
w18c16, wW18c17, W18c18, w18c19, w18c20, W18c21, w18c22, 
w18c23, w18c24, w18c25, w18c26, w18c27, w18c28}, 
{w19cl,w19c2, w19c3, wW19c4, w19c5, W19c6, W19C7, w19C8, 
w19c9, w19c10, W19c11, W19c12, W19c13, W19c14, w19c15, 
w19cl6, W19c17, W19¢18, W19cC19, W19¢c20, W19¢c21, W19¢22, 
w 19c23, w19c24, W19c25, W19c26, W19¢27, w19c28}, 
(w20cl, w20c2, w20c3, w20¢4, w20c5, w20c6, W20c7, w20c8, 
w20c9, w20c10, w20cl 1, w20c12, W20c13, W20c14, w20c15, 
w20cl6, w20c17, w20c18, w20c19, w20¢20, w20c21, W20¢22, 
w20c23, w20c24, w20¢25, w20c26, w20¢27, w20¢28), 
(w21cl, w21c2, w21c3, w21c4, w21¢5, w21c6, w21c7, w21c8, 
w21c9, w21c10, w21¢c1], w21¢12, wW21¢c13, w21¢14, w21cl5, 
w21c16, w21c17, W21c18, W21c19, w21¢20, W21¢21, w21c22, 
w21c23, w21c24, w21¢c25, w21c26, w21¢27, w21¢28), 
(w22cl, w22c2, w22c3, W22c4, W22c5, W22c6, W22c7, W22c8, 
w22c9, w22c10, w22c11, w22c12, W22c13, w22c14, w22c15, 
w22cl6, w22c17, W22c18, W22c19, w22c20, W22c21, w22c22, 
w22c23, w22c24, W22c25, W22c26, W22¢c27, W22c28), 
{w23cl, w23c2, W23c3, W23c4, W23c5, W23c6, W23c7, W23c8, 
w23c9, w23cl0, w23cl1, w23c12, W23c13, W23c14, W23c15, 
w23cl6, w23cl7, W23c18, w23cl9, w23c20, w23ce21, wW23c22, 
w23c23, w23c24, W23c25, W23c26, W23c27, w23c28), 
(w24cl, w24c2, w24c3, w24c4, w24c5, W24c6, W24c7, w24c8, 
w24c9, w24cl0, w24cl 1, w24c12, w24c13, w24c14, w24c15, 
w24cl6, W24c17, w24c1 8, w24c19, w24c20, w24c21, w24c22, 
w24c23, w24c24, w24c25, w24c26, W24c27, W24c28}, 
(w25cl, w25c2, w25c3, w25c4, W25c5, w25c6, W25c7, W25c8, 
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w295c9, W25c10, w25c1 1, w25c12, w25c13, W25c14, w25c]15, 
w25c16, w25c17, w25c18, w25c19, W25c20, w25c21, W25c¢22, 
w25c23, w25c24, W25c25, W25c26, W25c27, W25c28}, 
(w26cl, w26c2, w26c3, w26c4, w26c5, w26c6, W26C7, W26C8, 
w26c9, w26c10, w26c11, w26c12, w26c13, W26c1 4, w26c15, 
w26cl6, w26c17, w26c18, W26c19, w26c20, W26c21, w26c22, 
w26c23, w26c24, w26c25, w26c26, W26c27, w26c28}, 
(w27cl, w27c2, w27c3, W27C4, W27C5, W276, W27C7, W27¢8, 
w27c9, w27c10, w27c11, W27c12, W27C13, W27c14, W27C15, 
w27cl16, W27c17, W27c18, W27c19, W27c20, W27c21, W27¢22, 
w27c23, wW27¢24, W27¢25, W27¢26, W27¢27, W27c28}, 
(w28cl, w28c2, w28c3, w28c4, W28c5, w28c6, w28C7, W28c8, 
w28c9, w28c10, w28c11, w28c12, w28c13, W281 4, w28c15, 
w28cl6, w28c17, w28c18, W28c19, W28c20, W28c21, W28C22Z, 
w28c23, w28c24, W28c25, w28c26, W28C27, W28c28}}; 


a=({({0, alc2, alc3, alc4, 0,0, 0,0, 0,0, 0,0, 0, 0, 0, 0, 0) 
{a2c1, 0, 0, 0, a2c5, a2c6, 0,0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
(0, 0, a3c3, a3c4, 0, 0, 0, a3c8, a3c9, 0, 0, 0, 0, 0, 0, 0, 0}, 
(0, 0,0, 0, a4c5, a4c6, 0, 0, 0, a4c10, a4cl1, 0,0, 0, 0, 0, 0}, 
{0,0,.0,0, 0, 0,0, a5c8,-a5c9, 0, 0; adcl 2, adc) 3, 0, 0, 0, 0}, 
{0/0,0, 0, 0; 0: 0, 0,0) a6c10, a6ci |. 0, 0, 4614, a6cl 5, 0, 0}, 
{0,0..0, 0,0, 0-0, 0, 0, 0, 0, a7cl 2, a7cl3, 0,0, a7cl6, 0}, 
{0,0,0,0,0,0,0, 0,0, 0, 0, 0, 0, a8c1 4, a&c15, 0, a8c17}, 
(0, a9c2, a9c3, a9c4, 0, 0,0, 0,0,0, 0,0, 0,0, 0, 0, 0), 


{al0c1, 0,0, 0, al0c5, al0c6, 
(0, al 1c2, al 1c3, al1c4, 0, 0, 
{al2cl,0,0,0, al2c5, al2c6, 
{al3cl, 0,0, 0, al 3c5, al3c6, 


0,0, 0, 
0,0, 0, 
0,0, 0, 
0,0, 0, 
(0, al4c2, al 4c3, al 4c4, 0, 0,0, 0, 0, 


eocooo: 
oooone, 
ooo oc] . 
a 
a 
ogooococno, 
ooooo-:- 
OO SO Oo 
esses 


(0, 0,0, 0, 0,0, 0, al 5c8, al 5c9, 0, 0, al5cl2, alScl3. 0, 0, 0, 0}, 
(0,0, 0,0, 0,0, 0,0, 0,al6cl10, al6cl11, 0, 0, a16c14, al6cl15, 0, 0}, 
(0-0-0. 0: 0-020 al 7c8. a17¢9, 00-4171 274! 7e1 4-0, 0,0, 0); 
{0,0,0,0,0,0, 0,0, 0, a18cl10, al8cl1, 0, 0, al8cl4, al 8cl15, 0, 0), 
(0,0,0,0, 0,0, 0,0, 0, a19c10,a19cl11, 0,0, a19cl4, a19c15, 0, 0}, 
(0,0, 0,0, 0, 0,0, a20c8, a20c9, 0, 0, a20c12, a20c13, 0, 0, 0, 0}, 
(070, 0,0-070. 00 0.0 Oe aZziciZ.azicl3, 0, 0, aZlicl6.0}, 
(0200000. 0-0, 0 0, 0, 0.0, 00) a22c14,422c15, 0, a22c17). 
{0,0,0,0,0,0,0, 0,0, 0, 0, a23c12, a23c13, 0, 0, a23c16, 0), 
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(0,0, 0, 0,0, 0,0, 0,0, 0, 0, 0, 0, a24c14, a24c15, 0, a24c17}, 

(0, 0,0, 0,0, 0,0, 0,0, 0,0, 0, 0, a25c14, a25c15, 0, a25c17}, 
(0,0, 0, 0, 0,0, 0, 0, 0, 0, 0, a26c12, a26c13, 0, 0, a26c16, 0}, 

(0, 0, a27c3, a27c4, 0, 0, 0, a27c8, a27c9, 0, 0, 0, 0, 0, 0, 0, 0}, 
(0,0, 0,0, a28c5, a28c6, -1, 0, 0, a28c10, a28c11, 0, 0, 0, 0, 0, 0}); 


aherm = Transpose[Conjugatefa]]. 


= (Inversefaherm.w.a]).aherm.w.b 


Again, the lower case w was required to represent the weighting matrix 
W due to the Mathematica variable naming convention mentioned earlier 
(i.e., matrices have variable names which begin with lower case letters since 
Mathematica reserves variable names which begin with capital letters for 
built-in functions). The variable “aherm’ represents the complex conjugate 


transpose of the matrix A (ie., A¥), 


This code ran on a Macintosh I] computer which is equipped with five 
megabytes (MB) of random access memory (RAM) for about 25 minutes 
before halting due to an “out of memory” error. This same code was run on a 
different Macintosh I] computer equipped with similar hardware and a 
software package which allowed access of up to eight MB of hard disk space 
for use as virtual memory. Thus, Af#(fematrca had 13 MB of RAM available 
to it. On this 13 MB machine, the code ran for about 75 minutes before it 
again halted on an ‘out of memory’ error. Before attempting an alternate 
approach, the code was run in steps to see which calculation was causing the 
trouble. The multiplication of A4 Wand A ran successfully as did the 


multiplication of AH W. and b. The problem recurred when an attempt was 
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made to take the inverse of the product of A4\W, and A. We speculate that 
the problem occurs because the program is required to store and operate on 
a large number of string variables, thereby requiring large quantities of 
memory to store intermediate results. To confirm this, we attempted to take 
the inverse of a 17 by 17 symbolic matriz using only simple generic variable 
names (such as j!cl) without success. We must conclude that Mathematica 
requires too much memory to successfully run this type of symbolic problem 
on a personal computer. The methodology should be validated when less 
memory intensive or main frame based symbolic programs become 


available. 


In order to attempt the use of A/z(hematica 's LinearSolve function, we 
must first modify the matrix A so that it is a square matrix. This may be 
accomplished while maintaining the integrity of our weighted least squares 


formulation as follows: 


Recall that 


Axr=b. (198) 


Premultiplying both sides of (198) by the matrix W yields 


WAx-Wb. (479) 


Premultiplying both sides of (479) by the matrix A®™ yields 


At warx-aA4 won. (480) 


If we let 
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j- AU Wa, (481) 
and 
d- Au Wop, (482) 


then we have reformulated the problem as desired to 


jx-d, (483) 


where j is the 17 by 17 matrix defined by (481) and d is the 17 by 1 column 
vector defined by (482). 


The following revised Mathematica code was developed for this 


reformulated problem: 


b = (0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, G1, 0); 
w = (wicl, wlce2, wilc3, wic4, wlc5, wlc6, wlc7, wlc8, 
wlc9, wlcl10, wicl1, wicl2, wicl3, wicl4, wicl5, 
wicl6, wicl7, wicl8, wilcl9, wic20, wic21, wlc22, 
wlc23, wic24, wic25, wic26, w1c27, w1c28}, 
(w2cl, w2c2, w2c3, w2c4, w2c5, w2c6, w2c7, w2c8, 
w2c9, w2cl0, w2cl1, w2c12, w2c13, w2c14, w2c15, 
w2cl16, w2cl7, w2c18, w2c19, w2c20, w2c21, w2c22, 
w2c23, w2c24, w2c25, w2c26, w2c27, w2c28), 
{w3cl, w3c2, w3c3, w3c4, w3c5, w3c6, w3c7, w3c8, 
w3c9, w3cl0, w3cll, w3cl2, w3cl3, w3cl4, w3cl5, 
w3cl6, w3cl7, w3cl8, w3cl9, w3c20, w3c21, w3c22, 
w3c23, w3c24, w3c25, w3c26, W3c27, w3c28), 
(w4cl, w4c2, w4c3, w4c4, w4c5, w4c6, w4c7, w4c8, 
w4c9, w4cl0, w4cl1, w4c12, w4c13, w4c14, w4cl5, 
w4cl6, w4cl7, w4c18, w4c19, w4c20, w4c21, w4c22, 
w4c23, w4c24, w4c25, w4c26, w4c27, w4c28), 
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(w5cl, w5c2, w5c3, w5c4, w5c5, w5c6, w5c7, w5c8, 

w5c9, w5cl10, wW5c11, W5c12, wW5c13, W5Sc14, w5c]5, 

w5c16, w5cl7, w5c18, w5c19, w5c20, w5c21, w5c22, 

w5c23, w5c24, w5c25, w5c26, w5c27, w5c28), 

(w6cl, w6c2, w6c3, w6c4, w6c5, w6c6, w6c7, w6c8, 

w6c9, w6c10, w6c11, w6c12, w6c1 3, w6cl 4, w6cl15, 

w6cl6, w6cl7, w6cl8, w6cl9, w6c20, w6c2!, w6c22, 

w6c23, w6c24, w6c25, w6c26, w6c27, w6c28}, 

(w7cl, w7c2, w7c3, W7c4, w7c5, W7c6, W7C7, W7C8, 

w7c9, w7c10, w7cl 1, w7c12, W7ce13, W7c14, W715, 

w7cl6, w7cl7, w7cl8, W7c19, W7c20, W7c21, w7c22, 

w7c23, w7c24, W7c25, w7c26, W7c27, wW7c28), 

{w8cl, w8&c2, w8c3, w8c4, w8c5, w8c6, w8c7, w8c8, 

w8c9, w8&c10, w8c11, w8c12, W8c13, w8c! 4, w8c15, 

w8c16, w8cl17, w8c18, w8c19, w8c20, w8c21, w8ce22, 

w8c23, w8c24, w8c25, w8c26, w8c27, w8c28}, 

(w9cl, w9c2, w9c3, W9C4, W9C5, W9C6, W9C7, W9C8, 

w9c9, w9cl10, w9cl 1, W9c12, w9c13, W9c1 4, w9cl5, 

w9cl6, W9cl17, w9c18, w9cl9, W9c20, W9c21, W9C22, 

w9c23, w9c24, W9c25, w9c26, W9C27, W9C28}, 

(wl0cl, w10c2, w10c3, w10c4,w10c5, W10c6, W10c7, w10c8, 
wl0c9, w10c10, w10c!1 1, W10c12, W10c13, W10cC14, w10c15, 
wl0cl6, W10c17, W10c18, w10c19, W10c20, W10¢e21, w10c22, 
w10c23, w10¢24, w10c25, w10c26, w10¢27, w10¢28}, 
(wlicl, wllce2, wllce3, wl1lc4, wlicS, wllc6, wl 1lce7, wl 1c8, 
wlic9,wliclO,wlicl!, wllcl2,wlicl3, wilcl4, wlicl5, 
wWlicl6, wiicl7, wllc18, wl1c19, wl1c20, wl1c21, w11c22, 
wllc23,wl1c24.wl1e25, w11c26, w11¢c27, w11¢28}, 
(wl2cl, w1l2c2, w12c3, w12c4, w12c5, w12c6, w12c7, w12c8, 
wi2c9, wi2cl0, wl2cl1,wl2cl2, wl2ce13, w12c14, w12cl15, 
wi2cl6, wl12cl7, wi2cl8&, wl2e19, wl2e20, w12c21, w! 2c22, 
w12c23, w12c24, w12¢25, w12c26, w12¢27, w12c28), 
(w13cl, w13c2, w13c3, w13c4, W13c5, W13c6, W13c7, w13c8, 
w13c9, w13cl0, W13cl1, W13c12, W13c13, W13c14, w13c15, 
w13cl6, w13c17, w13c18, w13cl19, w13c20, W13c21, W13c22, 
w13c23, w13c24, w13c25, wW13c26, W13¢c27, Ww} 3c28), 
(wl4cl, wl4c2, w14c3, w14¢4, w14c5, w1 4c6, w14c7, w1 4c8, 
w14c9, w14c10, w14c11, wW14¢e12, W14c13, W14c14, w14c15, 
wW14cl6, w14c17, w14c18, W14c19, w14c20, W14ce21.w14c22, 
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w14c23, W14¢c24, W14c25, w14c26, w14¢c27, W14c28}, 

(wl 5cl, w15c2, w15c3, w15¢c4, w15c5, w15c6, w15c7, w15c8, 
wl5c9, w15c10, wI5S5cl1, wi5cl2, w15c13, w15cl4, w15cl5, 
w15cl6, w15c17, w15c18, w15cl19, w15c20, w15¢e21, w15c22, 
w15c23,w15c24, w15c25, w15c26, w15c27, w15c28}, 
(wl6cl, w1l6c2, w16c3, w16c4, W16c5, w16c6, w16c7, w16c8, 
wl6c9,w16cl10, wl6cl1, w16c12, w16c13, w16c14, w16c15, 
wW16cl16, w16cl7, w16cl8, w16c19, w16c20, W16c21, w16c22, 
w16c23, w16c24, w16c25, w16c26, wW16c27, w16c28}, 
(wl7cl, w17c2, w17c3, w17c4 w17c5, w17c6, w17c7, w17c8, 
w17c9, w17c10, w17e11, wW17¢e12, w17c1l3, W17cl 4, w17cl15, 
w17c16, W17c17, w17c18, W17c19, W17c20, w17c21, w17c22, 
w17c23,w17c24, W17c25, w17c26, w17¢c27, w17c28}, 
(w18cl, w18c2, w18c3, w18c4, w18c5, w18c6, w18c7, w18c8, 
wl18c9, W18c10, w18c11, w18c12, w18c13, w18c14, w18cl15, 
wW18c16, W18c17, W18c18, W18c19, w18c20, W18c21, w18c22, 
w18c23, w18c24, w18c25, W18¢c26, w18c27, w18c28}, 
(w19cl, w19c2, w19c3, W19c4, W19c5, W19c6, W19Cc7, W198, 
w19c9, w19cl10, W19cl1, W19c12, w19c13, w19c14, w19cl15, 
w19cl6,w19c17, W19c18, w19c19, w19c20, W19c21, W19¢22, 
w19c23, W19¢c24, W19¢c25, w19c26, w19c27, w19c28}, 
{w20cl, w20c2, w20c3, w20c4, w20¢5, w20c6, w20¢7, w20c8, 
w20c9, w20c10, w20c11, w20c12, w20c1 3, w20C14, w20c15, 
w20cl6, W20cC17, W20c18, W20¢19, W20c20, w20c¢21, w20c22, 
w20c23, w20¢24, w20c25, w20¢26, W20¢27, w20¢28}, 
(w21cl, w21c2, w21¢3, w21¢4, w21c5, w21c6, W21c7, w21c8, 
w21c9, w21cl10, w2Icl1, w21c12, w21¢13, w21c14, w2Icl15, 
w21cl16, w21c17, w21¢c18, w21¢c19, w21¢20, w21¢c21, w21¢22, 
w21c23, w21c24, w21¢25, w21c26, w21¢27, w21¢28}, 
(w22cl, w22c2, w22c3, w22c4, W22c5, W22c6, wW22c7, W22c8, 
w22c9, w22c10, w22c11, w22c12, w22c13, w22c14, w22c15, 
w22c16, w22c17, w22c18, w22c19, W22c20, W22c21, wW22¢22, 
w22c23, W22¢24, W22c25, wW22c26, W22¢27, W22¢28}, 
(w23cl, w23c2, W23c3, w23c4, W23c5, W23c6, W23c7, W23c8, 
w23c9, w23cl0, w23cl1, w23c12, W23c13, w23c1 4, w23c15, 
w23cl6, w23cl7, w23c18, w23c19, W23c20, W23c21, W23c22, 
w23c23, w23c24, W23c25, W23c26, W23c27, W23c28), 
(w24cl, w24c2, w24c3, w24c4, W24c5, w24c6, W247, W248, 
w24c9, w24c10, w24cl 1, w24c12, w24c13, w24c14, w24c15, 
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w24cl6, w24c17, w24c18, w24c19, W24c20, W24c21, wW24c22, 
w24c23, W24c24, wW24c25, w24c26, W24c27, W24c28)}, 
(w25cl, w25c2, w25c3, w25c4, W25c5, w25c6, W25c7, w25c8, 
w25c9, w25c10, w25c11, w25e12, W25c13, W25c14, w25c15, 
w25c16, w25c17, w25c18, w25c19, W25¢20, W25c21, W25c22, 
w25c23, w25c24, W25c25, w25c26, W25c27, w25c28), 
(w26cl, w26c2, w26c3, w26c4, W26c5, wW26c6, W26C7, W26c8, 
w26c9, w26c10, w26c1 1, w26c12, W26c13, W26c14, w26c15, 
w26c16, W26cC17, w26c18, W26C19, W26C20, W26c21, W26c22, 
w26c23, w26c24, w26¢c25, w26c26, W26¢27, W26¢28}, 
(w27cl, w27c2, W27c3, W27c4, W27C5, W27C6, W27c7, W278, 
w27c9, w27c10, w27cl1, w27c12, W27cC13, W27c14, w27cl15, 
w27cl6, w27c17, w27c18, w27c19, w27c20, W27c21, W27¢22, 
W2/c2 3. W2/7C24, W2/7C2Z), W2Z/C20, W2/C27, W27C28). 
(w28cel, w28c2, W28c3, w28cC4, W28c5, W286, W28C7, W28C8, 
w28c9, w28c10, w28c1 1, w28c12, W28c13, W28cC14, W28C15, 
w28c16, w28c17, w28c18&, W28c19, w28c20, W28c21, W28C22, 
w28c23, w28c24, W28C25, W28c26, W28C27, W28C28)}}: 


= (0-ale2. alcs, alc4, 0,.0,0. 0,0, 000) 0.0, 0.0, 0, 0} 
{az2cl, 0,0; 0. a2¢5, azc6, 0.0, 0, 0. 0,0; 0, 0, 0, 0,0}, 
(0, 0, a3c3, a3c4, 0,0, 0, a3c8, a3c9, 0, 0, 0, 0, 0, 0, 0, 0}, 
(0,0, 0,0, a4c5, a4c6, 0, 0, 0, a4c10, a4cl1, 0,0, 0, 0, 0, 0}, 
(0,0, 0,0, 0,0, 0, aSc8, a5ce9, 0, 0, a5c12, a5c13, 0, 0, 0, 0), 
(0,0; 0, 0,0; 00.0, 0, a6e10: a6¢1 1.0: 0, a6c14, abel 5, 0, 0}, 
(0,0, 0, 0,0, 0,0, 0, 0, 0, 0, a7c12, a7c13, 0, 0, a7c16, 0}, 
(0, 0,0, 0,0,0,0,0,0, 0,0, 0, 0, a8c14, a8c15, 0, a&c17}, 
(0, a9c2, a9c3, a9c4, 0,0,0, 0,0, 0,0, 0,0, 0, 0, 0, 0}, 


a9 
fal0cl<0; 00.4105, al0cb, 0-0. 0.0.0.0, 0.0/6, 0.0) 
(al te? al ics.al les. 04070-0,0 0, 0.00 0.0.0.0). 
{al2cl, 0, 0,0, al2c5, al2c6, 0,0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
fal sel. 0, 0, 0) al 5c5 -alse6.0, 0.0, 0. 0,0,0.0,0. 0. 0} 
(0, al 4c2, al 4c3, al 4c4, 0,0, 0,0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 


(0, 0,0, 0,0, 0, 0, al 5c8, al 5c9, 0, 0, al 5c12, al 5c13, 0, 0, 0, 0), 
{0)-0;-0, 0,0, 0, 0, 0, G,aloci0, alocll, 0, 0.a16el4, al6ct5, 0, 0}, 
(0,0, 0,0, 0, 0, 0, al7c8, a17c9, 0, 0, a17c12, al7c13, 0, 0, 0, 0}, 
(0,0,0,0,0, 0,0, 0, 0, a18c10, a18c11, 0,0, al8c14, al8c15, 0, 0), 
(0,0,0,0,0,0,0,0,0,a19c10, a19cl1, 0,0, a19c14, al9c15, 0, 0}, 
(0,0, 0,0, 0,0, 0, a20c8, a20c9, 0, 0, a20c1 2, a20c13, 0, 0, 0, 0}, 


Paz lol2, aZ icin 0 sOracc lo, On 
0,0, a22c14, a22c15, 0, a22c17), 
a23c12, a23c13, 0, 0, a23c16, 0}, 
, 0, 
0; 


? ry 


SSeS ote 


SsSoecocoeecs 


0, a24c14, a24c15, 0, a24c17), 

0, a25c14, a25c15, 0, a25c17), 

,a26c12, a26c13, 0, 0, a26c16, 0}, 

‘ Oia a27c9, 0, 0,0, 0, 0, 0, 0, 0}, 
, 0, 0, a28c10, a28c11, 0,0, 0, 0, 0, 03); 
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aherm = Transpose[Conjugate[a]]: 
j = aherm.w.4a; 
d = aherm.w.b; 


LinearSolvelj,d] 


This revised code ran on the five MB Macintosh for about 25 minutes 
before halting on an “out of memory’ error. It also ran on the 13 MB 
Macintosh for about 110 minutes before halting on the memory error. These 
failures have led us to conclude that the solution of the general problem is 
possible with this technique but is not practical with currently available 


hardware/software configurations. 


One additional test was run using the code generated for the general 
case employing Afathematica 's LinearSolve function. In this test case, the 
three media waveguide with plane, parallel boundaries was simulated by 
setting the weighting matrix, W, equal to the identity matrix (using the 
Mathematica command: w = IdentityMatrix[28];), setting appropriate values 
of the A matix equal to zero, and directly assigning values to the Al matrix - 


as follows (in Afztfematica code): 
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aherm = ({0, 0, 0, 0, 0, 
0,0,0,0,0 


00 00070205070) 00 -050.05 0.0, 0, 
iO), 


{alc2Ccc, 0, 0.0, 0, 0,0, 0, age2CC; 0, 0,.0,0, 0,0, 0, 0.0, 


0,0, 0, 0, 0, 0, 0, 0, 0, 0}, 

fale3cc, 0, asc3CC, 0,0, 0,0, 0, a9c3CG, 0, 6,0, 0, 0, 0, 
0,020) 0-070) 0.0.0, 0. 0-a27c5CG, 0). 

{alc4CC, 0, a3c4CC, 0, 0,0, 0, 0, a9c4CC, 0,0, 0, 0, 0, 0, 
0,0,0,0,0, 0,0, 0, 0, 0, 0, a27c4CC, 0}, 

(0,0, 0,0, 0, 0,0, 0,0, 0,0, 0,0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0:0,0;07 0,0}, 

(0,0, 0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0, 0, 0, 0, 0, 0, 0, 
Os050" 0050), 

(a0 0.080 0.0,-0.0-0. 0, 0200-0, 0.050, 6-0, 0: 0, 
OsC 3080 00): 

(Oh Op ascece, 0) ascsCC, 0.0, 0;0,0, 0, 0,0, 0, al ScsCc, 
00. 0-0, 0,0, 070, 0,0, 0-aZ7c6CC, 0), 

(Ge a3c9CC, 0 ascoCG. 0. 0:0, 000, 00,0, 0; al 5c9CC. 
0,0, 0,0, 0, 0, 0, 0, 0, 0, 0, a27c9CC, 0), 
(0.0, 020) 07020) 0,0-0,0. 0, 0,0,.0 0-0, 0.0.0, 0, 0; 
O70; 0;,,0,.0..03. 

(O10) 07020200; 00-0 70,0, 0-0, 0.0, 0070, 02 0,0. 
0,0, 0, 0, 0, 0), 

(000-050. ayel 2Z6C 0-0, 0..0, 0.0.0, 070 ai sci2cc 0,0), 
O2070.070. 00.070. 070) 

(O0F 0502000, 070,00, 0 070-0- 0, 0.0, 0202), 0, 
0.0, 00,0. 0}. 
(0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 0,0, 0, 0, 0, 
070.0, 020,03. 

(0,0, 0,0, 0,0, 0, 0,0; 0, 0, 0, 0, 0; 0; 0, 0,0, 0,0, 0, 0, 
0,0, 0, 0, 0, 0}, 

(070, 0.0. 0.0. 0.0.0. 0-0, 6.0: 0,0-0,0, 0. C0 0.0.0, 
0,0, 0, 0, 0, 03), 


where the notation CC in each of the variable names indicates that the 


element represents the complex conjugate of the appropriate element of the 


A matrix. 


ae) 


This code ran on the five MB Macintosh for ten days without converging 
to a solution. This additional failure confirms our suspicions that the general 


problem is not practically solved using current technology. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The main purpose of this thesis was to obtain the symbolic solution of a 
multilayer (four fluid media) ocean waveguide problem. The waveguide was 
assumed to have depth-dependent ambient density and sound-speed 
profiles in all fluid media, and arbitrarily shaped boundaries between all 
fluid media. A system of 28 equations in 17 unknowns was generated by 
satisfying all of the boundary conditions (including the boundary condition 
at the source) in cylindrical coordinates. A weighted least squares estimation 
technique was employed to formulate a symbolic solution to this 
overdeter mined (more equations than unknowns) case. A computer program 
capable of performing symbolic algebra was sought to minimize the number 
of assumptions required to be made, thereby maximizing the generality of 
the solution obtained. M/a(hematica (version 1.2.1 £33 (enhanced)) for the 
Macintosh computer was selected for this work due to its availability at the 
Naval Postgraduate School and its symbolic algebra capabilities. Ma/hematica 
code was developed which programmed the weighted least squares 
estimation technique for the most general case. Unfortunately, this code was 
unable to provide a solution to the most general case due to software and 


hardware limitations (ie., speed and random access memory problems). 


By relaxing the arbitrary boundary shape assumption, Afathematica 
code was developed which programmed a direct solution to the three media 
waveguide problem for plane, parallel boundaries. This code ran 


successfully, and provided results which could be verified bv direct 
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comparison with the known solution to this classical (i.e., three fluid media 
with plane, parallel! boundaries and constant ambient density and sound 
speed in each medium) waveguide problem. During this process, built-in 
Mathematica functions were used in an attempt to simplify the resulting 
symbolic expressions. This effort revealed that these built-in functions 
lacked sufficient sophistication for applications of this complexity. This lack 
of sophisication resulted in manual reduction of the program output so that 
verification was possible. A/zifematica code was then developed to solve the 
four media waveguide problem for plane, parallel boundaries. This code also 
ran successfully and yielded results which could be verified using Known 
classical waveguide solutions when some judicious assumptions were made 


to mathematically eliminate the fourth medium. 


In addition to validating the symbolic solution technique, this thesis 
provides a series of generic expressions for the unknown constants for each 
of the three and four media waveguide problems with plane, parallel 
boundaries. Each of the generic expressions is a combination of generic 
variables whose definitions are provided in the text. Each of these generic 
variables can be programmed in a high level language (i.e, FORTRAN) as a 
unique subprogram or function. In this manner, the unknown constants can 
be calculated by combinations of calls to appropriate subroutines. This 
modular programming technique is enhanced by the fact that each of the 
generic expressions has a common denominator, which can also be 


programmed in a similar manner. 
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The generic results from the four media waveguide problem for plane, 
parallel boundaries should be programmed in FORTRAN. When this 
programming has been completed, the standard test cases (i.e., three fluid 
media with plane, parallel boundaries and constant ambient density and 
sound speed in each medium) should be run on the new code to verify the 
results. This verification will provide additional credibility to A7athematica 's 


output. 


It is recommended that the most general case (i.e, four fluid media 
with arbitrarily shaped boundaries) be attempted again when one of the 
following conditions are met: 

¢ a Macintosh computer with more than 13MB of RAM becomes available 
at the Naval Postgraduate School, 
¢ an advanced version of Mathematica is released, or 


¢ another computer program with symbolic algebra capabilities becomes 
available for use on a workstation or main frame computer. 
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